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ABSTRACT 


Ship Synthesis Model for Coast Guard Cutters 

by: Michael Jay Goodwin 

Submitted to the Department of Ocean Engineering on 

May 9, 1975, in partial fulfillment of the requirements 

for the degree of Ocean Engineer and the degree of 

Master of Science in Naval Architecture and Marine Engineering 


In order to determine the optimum design of a Coast 
Guard cutter, it is necessary to perform a large number 
of feasibility studies. This is an impossible task if 
done using hand calculations. A mechanized version of 
these hand calculations is developed which will permit 
the calculations to bse performed by a digital computer. 

The final program listing is given together with the 
development of the empirical relations used. 

The program is applicable to search and rescue and 
patrol type cutters with lengths between 150 and 400 feet. 
Twin screw propulsion is assumed. 

A check on the accuracy of the program is made using 
three U.S. Coast Guard cutters now in service and the 
U.S. Navy*s Patrol Frigate design. The program is found to 
be much more accurate tnan existing programs when applied 
to Coast Guard vessels. The program predicts the full load 
displacement within four percent for the Coast Guard designs 
tried. 
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ChAPTsr I 
TivPrODUC TION 


Altnough digital conputers have been used by engineers 
for more than twenty years, their potential as design 
aids is only beginning to be realized. Naval architects 
made use of computers soon after they became available 
for solving specific algorithms such as displacement sheet 
calculations and speed-power estimates. Such algorithms 
are beneficial in the preliminary and later design stages. 
However, witn the growth of systems engineering in ship 
design, more emphasis is being placed on the conceptual or 
earliest phase of a design. At this stage of the design, 
the computer can be used more effectively than at any other 
stage. | 

3ecause of the comolexity of the naval architectural 
Calculations involved in later design stages, the principle 
characteristics of a snip must be chosen early in the design. 
The outout of the conceytual design phase is this set 
of characteristics. Without the use of a computer, the 
best that can be done is to develop a feasible set of 
Characteristics when a optimum set is desired. Systems 
engineering is placing more emphasis on obtaining this 
ootimum,. 

To Choose an optimum or near optimum set of character- 
istics, a great many feasibility studies are required. It 
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ieimoossrbie to do more than a few feasibility studies 
without the aid of computerized models because of the 
time required. At this early stage, the ship must be 
designed as a whole. This precludes dividing the work 
into a number of parallel paths with several people 
working simultaneously on a Single feasibility study. 
gach study is essentially a one man job. 

Witn hand calculations a good engineer can perform 
about one feasibility study per day. As a result, great 
reliance must be placed on the deSigners intuition. In 
the past this has led to series of designs which may or 
may; not have been optimum. Zach new design was a slight 
perturbation of a vrevious successful design. While this 
procedure produces satisfactory designs, bad features 
tend to be duplicated and few new features are added. 
The engineer was not really at fault because no tool 
was available to perform all the calculations necessary 
to obtain an optimum design. 

Comouters permit the mechanization of these hand 
' Calculations. The software required will be referred to 
here as a ship synthesis model. The model aids the designer 
by freeing him from the routine calculations which take 
up most of his time during the early phases of design. 

The algorithms referred to previously also free the 
HeSsronet, mom tedious Calculations. 
what maxes a ship synthesis model different is that 


it goes beyond merely reducing the time and effort required 





to do calculations woreviously done by hand. Using the 
model, hundreds or even thousands of feasibility studies 
can be made. While this still does not guarantee that 
the optimum solution is obtained, the engineer will have 
much more confidence in his design than he would if only 
a few feasidility studies had been made. The designer 
is free to concentrate on the more important asvects of 
the early stages of ship design such as systems level 
tradeoff decisions. These include hull form, sveed, 
endurance? and payload ootions among others. These aspects 
were seldom considered in detail before the advent of 
Ship synthesis models. 

=qually important is the fact that computers permit 
a consistency among feasibility studies. With several 
people doing indesendent feasibility studies or even 
With 2 Single person doins them all, it is difficult to 
insure tnat tne same assumptions will be made for each 
study and that no bias is apvlied. The model guarantees 
consistency. And consistency is important because without 
it two feasibility studies can not be compared and an 
Spernum Solution 1S difficult to select. 

several models have been available for a number of 
years and have proven their value. The U.S. Navy's 
destroyer model, ppo7 ‘1 )*, on which the program developed 
in this thesis is patterned, was started in the 1960's 

* Numbers in varenthesis refer to references at the 


end of this thesis 





10 
anc updating work has been continuing since then. This 


was the first true "systems" model. iiodels also 

have been develoved which attempt to model all surface 
ship types with a single program. For accurate results, 
each different ship type really requires its own model. 

in addition to tne destroyer model, the U.S. Navy 
has developed individual synthesis models for aircraft 
carriers, underway replenishment ships, submarines, 
hydrofoils, and several other ship types. The U.S. Coast 
Gu2rd presently has no model which gives good results for 
cucters. The destroyer model, DDO7, does not apply to 
ships as short as 150 feet in length which is one Coast 
Guard requirement. Some all inclusive models have been 
tried but with voor success, with predictions of displace- 
ment as much as 30 percent off when data on existing 
designs is inputed. 

The prozram, which will be referred to as the cutter 
model, d2veloped in this thesis is an attempt to satisfy 
the Coast Guard’s need. Tne cutter model applies to 
searcn and rescue and patrol type cutters with lengths 
between 150 and 400 feet. It is not applicable to buoy 
tenders, tugs or ice breakers. An effort has been started 
by the Coast Guard to develop a model for buoy tenders 
anc tugs but this work is not yet complete. 

sven in tne early stages of design where the synthesis 
models are of use, a great many calculations are required 
to prove the feasidility of a ship design. This has led 


designers to employ a type of regression analysis where 
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data from existing shivs is scaled up or down in order 


to estimate values for a new design. This estimation 
procedure is employed in ship synthesis models also. 
fore detailed calculations could be made but at the cost 
of sreatly increased computer time since a large number 
of feasidility studies will be run. 

Tne shortcut estiuating techniques employed place 
more empnasis on the relative accuracy among studies 
rather tnan on the absolute accuracy of the result. 
Because of this, all estimating relationships in the 
mocel must show the devendence of the result on all 
Sisnificant inout variables. eee sone this is an im- 
possible task because not enough is known about this 
derendence. As models become more widely used much 
development work will have to be done in this area. An 
estimate of the variation with the most significant 
input variable is the best that can be done at present. 

secause the estimating relationshivs are rather 
crude, ectimates are only made at the ship system level. 
Snio sysven level information consists of such elements 
as total outfit weight, total volume required for machinery 
and required electrical load for ventilation. Although 
total arranzements area iS computed, no attempt is made 
to verform an arrangement of spaces or to compute trim. 

The engineer usually has sufficient flexibility in these 


areas to arrive at a feasible design if sufficient space 


is availaole, 
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One problem in developing a model for Coast Guard 


cutters is that there are few designs available from 
which to obtain data. This causes difficulty in deter- 
Mining rational scaling laws. In many cases the use of 
naval architectural vorinciples indicates the correct 
scaling orocedure. In other cases all that can be done 
m to searcn for the most logical correlation. 

The design produced by a synthesis model should 
Pecenored yO meeG tne conditions required of a feasible 
solution. fFrirst there must be a balance between weight 
and displacement. Second, internal svace available must 
equal internal space required. .This requirement is becoming 
very important as most Coast Guard vessels are volume 
rather tnan weight limited. Third, the energy available 
must equal the energy required. This includes both 
propulsion and electrical power requirements. finally, 
tne distribution of weight and volume must be such as to 
satisty design criteria for transverse stability, girder 
strength and sea keeping. There are also other require- 
ments imolicit in the solution such as hull type, material 
or construction and space requirements. 

In zeneral, the inputs required for a ship synthesis 
model are a description of the “payload“ to be carried, 

a statement of vehicle mobility performance and a selection 
of options to be used. The payload description must be 
more detailed than a simple weight breakdown. The relative 


Vertical centers of gravity, area requirements, electrical 
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power requirements and cost must also be given. Payload 
for tne cutter model also includes people, and thus the 
size of the crew is an input. The description of vehicle 
motility performance includes booth maximum calm water 
speed and endurance speed and range requirements. Options 
include length, Cys Cs; as well as machinery type and 
marzins. ‘Ynere is also an ovtion in the cutter model 

to inout tne maximum sustained horsepower, in which case 
the program will calculate the maximum sustained speed, 

A full description of the vrosram inputs for the cutter 
model is contained in Chapter III (Main Program Description) 
and in Apvendix A. 

The output zenerated by the model must be sufficient 
to demonstrate the feasibility of the design. Wherever 
possible, the output should permit the designer to verify 
the feasibility without doing additional calculations. 

On a Tew of the studies run and in particular on the design 
selected, the design team will have to lay out a rough 
arrangements drawing to verify that sufficient space is 
"available. Suificient information should be given in the 
Survout to vermit this to ove done easily. 

the output generated should also allow the designer 
to rapidly determine the design that best satisfies his 
Optimization criterion. This criterion usually takes the 
Meorm of a total acquisition cost or life cycle cost. 

*’ne definition of coefficients and dimensions used 


in this thesis may be found in Reference (3). 
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Predictions of acquisition cost are often included ina 
ship synthesis model. 

Sinally, the output sSnould provide a basis for beginning 
the preliminary design phase. Generally the preliminary 
pnase is handled by a number or design teams working in 
parallel. ‘Tne output produced must provide information 
to initialize each of these parallel developments. 

sample output listings for the cutter model are 
ziven in Avpendix 3B. 

It is important to emohasize that the cutter model 
oroduces a zreasible solution only and not an optimum 
Semcon, An Optimum solution implies that the best of 
a very large number of feasiole ships has been chosen, 

It is felt that the desizn process works more efficiently 
mmeemie CeGSlSner acts aS a Buide for the computer rather 
than having the computer do the entire optimization on a 
Sinzle run. in this way, many unnecessary deSigns can 
avoided in areas where designs are less attractive than 
others already tried. suck a process takes more real 
tite Dut less computer time. It also gives the designer 
a much better feel for the sensitivity of the optimum 

to changes in characteristics. 

The remainder of tnis thesis is arranged as follows: 

Chaoter iI gives an overview of the cutter model 
and exdlains tne major options. 

CUnapter III goes into detail on each of the sub- 


routines in the model and explains the derivation of all 
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of the empirical formulas used. This chapter may be 
skioped by anyone not interested in the details of the 
prosram. 

wnaoeen L¥ @ives an evaluation of the model including 
test runs made using existing ships. 

Uha.vter ¥ contains recommendations and conclusions. 

Appendices aA and 3 contain input and output formats 
respectively. 

Appendix C is a listing of the program. 

The remaining appendices give background information 


us2=d in this work. 


It is recommended that at least Chapters I, II, IV, V 
and Appendix A be read before attempting to use the program, 
Chapter III is aimed primarily at the person who wishes 
to inodify the program. However, the important aspects 
of each subroutine are covered in the introduction to 


Peen Oo: the subsections of Chapter III. 





CHAPTER II 
GENERAL PROGRAM DESCRIPTION 


This chapter zives an overall description of the 
solution method used bv the cutter model but does not 
discuss the detailed calculations. The detailed calculations 
are explained in Chapter III. This chapter should also 
serve as a guide to understanding how each subroutine 
discussed in Chavter III fits into the total program. 

The organization of the computer program can be 
viewed in two different ways. The first is the program 
control organization which is shown in Figure 1. Sub- 
routine X=CUi2 controls the entire solution of each 
feasibility study case. control is returned to the main 
prozram only when data for a new case is required. 
suoroutine X=CUT=s calls other subroutines as needed and 
graduall: fills all the data storage locations needed 
for output. This subroutine also controls the iteration 
joons for weight and vertical center of gravity. 

The lozical organization of the program is shown 
schematically in Figure 2. This is an easier organization 
scheme to follow and shows more clearly how the program 
solves each invut case. It will be noted that subroutine 
AsvuUTc has not been shown in figure 2. As explained 
above, X:CUT= merely guides the computer around the 
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Figure 2 
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lozical organization chart. Srror messages are also 
vrinted by subroutine X&CUTZ if any constraint of the 
model is violated. 

The program begins with data input in the MAIN 
program. The primary purpose of the mAIN program is to 
read date cards and all input is done in this routine. 
2ome booxkeeping is also done in MAIN such as summing 
payload weights and costs. ‘The final purpose of this 
routine is to estimate initial values of vertical center 
of gravity, KG; displacement, and cubic number. 

Once all data for one feasibility study has been 
read by iwAIN and the initial values of KG, displacement and 
cubic number have been determined, control is transferred 
momeuorousine A2CUTE. This subroutine directs the execution 
of the program until a solution is obtained or until 
pmeonaetralnt 1S violated. When either of these occur 
control reverts to the MAIN program where data for the 
next feasibility study is read. If there is no more 
data the vrogram stops. 

The first subroutine calied by XECUTe is an under- 
waver dimensions subroutine, UWDIM. Subroutine UWDI! 
uses the input values of length, free surface correction, 
Co» C, and Gi;/Seam aS well as the estimated KG and dis- 
placement to calculate the required beam, draft and 
waterplane coefficient, Cyn. To the level of accuracy 
of this program these completely determine the under- 


water Shave of the hull. 
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Two options are available with regard to speed and 


power calculations. Under the first option the maximum 
sustained sveed of the vessel is input. If this option 
is speciried, X:CUTS calls Subroutine HrcALO which, 
together with subroutine CXVAL, estimates the horsepower 
“required to maintain the maximum sustained speed. 

Option two allows tne user to specify the maximum 
sustained horsepower. In this case the program will 
compute the maximum sustained speed. Subroutine HPCALC 
has been written to take speed as an input and return 
avalue of horsevower. Therefore, an iterative procedure 
is required if this second option is used. Control of 
this iteration takes place in subroutine X=CUTE. 

Regardless of which of the above options is used 
the endurance speed is input. Subroutine HFCALO is 
called to determine the power required at endurance speed. 

A Taylor standard Series horsepower estimate is 
used in subroutine HPCALC. Stored values of Cry are 
used. The purpose of subroutine CRVAL is to choose the 
correct C, value from a stored array. Reference (7) 
exdlains this method for estimating horsepower. 

The electrical power routine, EPLANT, is the next 
one called. This routine determines the maximum electrical 
load at sea, the average electrical load and sizes the 
generators. The average load is needed for estimating 
generator fuel consumption. The input electrical loads 


along witn crew size and cubic number are used in this 
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routine. 


The weight of fuel, lubricating oil and potable 
water is next calculated in subroutine MACHLQ. The 
endurance range given as input together with crew size 
and the endurance power calculated earlier are used by 
this routine. lLubd oil requirements are either input, 
in the case wnere horsepower is specified, or are estimated 
based om specified machinery type. Potable water require- 
ments devend entirely on the crew size. If make-up feed 
water is required it should be specified as a cargo load 
item. The specific fuel consumption is either input, in 
the case where horsepower is specified, or is estimated 
based on machinery type and horsepower, 

subroutine WBSIZS is next called if the sustained 
sneed 1s specified. Tnis routine estimates the machinery 
box ieee and minimum depth. If the horsepower input 
ootion is specified these dimensions are input. 

tne next subroutine is one of the most important 
in tne program. Subroutine VOLUIKiz& determines both the 
Sheer line and the deckhouse volume. These complete 
the definition of internal volume and arrangements area 
Since the underwater volume has been determined in sub- 
routine UwWDIix. Because most ships for which this program 
applies are volume rather than weight limited, it is 
important to insure tnat there is sufficient aeeee Or 
all ship functions and equipment which must be installed. 


ane volume is divided into three categories. First, 





ac 
there is the engine room. The dimensions for the engine 


room have been estimated earlier. Second is tankage 
volume and finally arrangements volume. Arrangements 
volume is more useful if converted to an arrangements 
area. 

As a rirst step in determining the space available 
in each of these categories, the sheer line must be deter- 
Mined. subroutine VOLUM= calls subroutine SHEER to 
determine an acceptable sheer line. The sheer line must 
eet several criteria including adequate freeboard forward 
and aft, adequate depth. amidships for the engine room 
ani for adequate strength, and an acceptable sneer line 
eurveture. 

Wwitn tne sheer line fixed, the total volume of the 
hull is estimated. The volume of tne engine room is 
also estimated and then subtracted from the total volume. 
Next, an estimate is made of the fraction of the remaining 
volume wnicn can be used only for tankage because it is 
unsuitable for use as arrangements volume. 

A required tankage volume is then estimated based on 
tne weients oF fuel, luod 011, votable water and aircraft 
fuel. Tnis required tankage volume is compared to the 
volume which must be used for tankage. The larger of the 
(Wo 1S suotracted from the hull volume that remains. The 
remaining; hull area is available for use as arrangements 
volume, This volume is converted into an arrangements area. 


A required arrangements area 1S computed based on 
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empirical data from past designs and on input area values. 
Baca justing the Size of the superstructure or by adding 
a raised deck, the available arrangements area can usually 
be made to equal the required arrangements area. If 
Boo 2arze a deckhouse is required or if a raised deck 
Jonzer than tne length of the ship is needed, an error 
message is generated and a new input case is tried. 

Once tne volume of the hull has been determined, 
the weignt of all the light ship weight groups and of the 
Joad items can be calculated. This is done in subroutine 
WIGHT. Altnough many of the items calculated by this 
routine are later printed as output, the variable of 
immediate importance is the full load weight, WFULLD. 
This value is compared with the initial guess of displace- 
ment. If there is a difference of greater tnan 2 tons 
between the two, a new estimate of displacement is made 
and the orogram returns to subroutine UWDI:i for a new 
eivimate ot dimensions. Once a balance between WFULLD 
and displacement 1s achieved the program proceeds to the 
Seftimation of vertical center of gravity. 

Subroutine VOG maxes this estimate. iost centers 
are vaxen as a fraction of tne average depth of the hull 
or o> the amidships devntn. Tne estimated vcg is compared 
to tne initial guess and a new KG estimate is made if 
tnere 1S no azreement within a tenth of a foot. Once 
peoalancs is obtained, the program proceeds to a cost 


estimate. 
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Tne lead Ship cost is estimated using a Coast Guard 
cost estimating procedure known as "Flanagan's lie thoa" §6) 
This orocedure has been vrosrammed as subroutine COST. 

An estimate or tne sveed the vessel will be able 
to sustain in the North Atvlantic Ocean is made in sub- 
routine scaAssyD. This routine is based on work done at 
the naval onid =nzgineerinsg Center as reported in References 
(4) and (5). 

Sinally, subroutine OUTFUT is called. The purvose 
of this routine is to print the values determined by the 
other routines. The input data is also printed so that 
a compact mecord of the feasibility study is produced. 
After the output is printed, tne program returns to the 
mAlix prozram to read more data and start a new case. 

This concludes the general description of the program. 
A twnucn more detailed description is given in Chapter III. 
ire read2r wno is not interested in the details of the 


Drogram can skiv Chapter IIi with no loss of continuity. 





Sear ioR Lit 


DITAILAZD DESCRIFTION OF THE PROGRAL 


Met Introduction 

Tnis chapter is composed of sixteen sections including 
this introduction. The other fifteen sections each 
exolain in detail one subroutine of the cutter model. 

For the location of these subroutines in the overall 
program refer to Figure 2, Logical Program Organization. 

An attemot has been made to keep each section in the 
same format. A brief description of the routine will 
be given first followed by a list of required inputs. 
Statements will then be exvlained one by one. A list 
of the outputs produced by the subroutine will be given 
next. Finally, a nomenclature list will be given. The 
nomenclature list applies only to the subroutine being 
considered. in some cases the same variable will have 
different names in different subroutines. 

In both the input and output listings a code is 
employed to indicate whether the variable is transferred 
by a labeled common statement, LC(AA), or aS a Sudroutine 
Bai. dumny argument, SCDA. The letters in parenthesis 
after the labeled common code indicate which labeled 


common was used to transfer the argument. 
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3.2 MAIN Program 


Sae.i Introduction 

The MAIN subroutine performs several functions. 
First and most important is the reading of all input 
data used by the model. Each feasibility study case 
starts with the reading of one or more data cards by 
the MAIN program. The actual execution of each case 
is controlled by subroutine XECUTE. This routine is 
discussed in Section 3.3. However, once execution on a 
case is complete, control 1S returned to the MAIN sub- 
routine. A check is then made to see it there is more 
miput*data. If not the MAIN subroutine directs the 
computer to stop. 

Several bookkeeping tasks are performed in this 
Subroutine also. These include such jobs as adding 
together the deck area and electrical load requirements 
of the individual payload items. The final function 
pertormed is the setting of initial values for displace- 
ment, KG and cubic number. These values are used for the 
first iteration and are then replaced by better estimates. 
See.e Inputs 

Block Data subroutines have not been used in this 
model. AS a result there are no variables with assigned 
values at the beginning of the MAIN program. All of the 
variables required are read from data cards or assigned 
values in the routine. 


A flow chart for this routine is given in Figures 
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3a and 3b. A listing of the program is given in Appendix 


C. The nomenclature list at the end of this section 
gives a description of each variable name used. 
3.2.3 Statement Descriptions 

A line by line discussion of the program statements 
follows. The REAL and COMMON statements have not been 
listed but may be founu in Appendix C. 

REAL Statements - An explicit type definition has 
been employed so that the variable names will be more 
descriptive of their actual meanings. 

COMMON Statements - Because a great many variables 
are required to define each case, COMMON statements have 
been used to transfer values between subroutines. This 
also saves computer storage space. All COMMON statements 
used are labeled common. 

Cr array input: 

READ(G,120) CR1 

120 FORMAT(16F5.0) 

READ(8,121) CR2 

mei FORMAT(9F5.0) 

READ(S,122) CR3 

122 FORMAT(6F5.0) 

The three arrays which contain the Taylor’s Standard 
Series C, values are read only once at the beginning of 
the program. 

INDEX input: 


2 READ(8,124) INDEX 
124 FORMAT(I1) 


The variable named INDEX has two functions. First, 


as described below, it tells the program when to stop. 





MAIN SUBROUTINE 
FLOW CHART 








REAL ANDO 


LABELED 
COMPTON 


STATEMENTS 





START MAIN 












READ 
HEAD (TZ), T=1,2¢ 
NOARM 





READ TL=/, NVoARmM 
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READ NOELT 
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Figure 3a 








READ INDEX 





WCARGO =D 
WTG7 +f 
WTAMMO= 
WTAC = g 
ELOAO? g 
LOAD = 






ALEAOH*D 
ALEAHL = GB 
esreg7=¢ 







WG = 
‘CSTGIG eg 






FOR. Lel,MoeLT 
WTGds= ES GLIWTCI) 
ELOAD@ & Expr) 
GLICST ¢ GRICST 7/000, 
CSTGI2e EGELACSTID) 
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NOFF =F LOATMCOFF) READ DdOLHR, READ 













WC Po =Floar(Mecho) GIIND G@SInD ENDDAY 
G2IWD GEIND WACE UL 
VENML 2 FLOR TWOCkEW) G3INO CACFUL 










READ HFo"'S, 
DHA, XLN, WTGZ. 
SFCMHP,SFCHHP, 
VEX D, RGEAD,FRLWAD, 
PCE, Pom, CSTE2 , WILO 







YES 


READ JOoPT 
Ci CompurTe MAX HP, 


2 Comfure SPEED) 





HPMS:@ 






NO Nore: 

















READ VSUS, ATNPE =U , eeee 
VEND, RGEND, BEE Pe) 
HPMS = HPNSH+AESE. Poe , PCM 






CS1G 22 CSTGL BMW. 






READ MTYPE, 
FRSC, LEN, CP, OX 
GmMMVL , DBMAE, 

DELCE 


PREVIOUS 
PAGE 
Alo 
CN= DPTLYx/40. 


CALL XECUTE 


GR2¢esrTre CSTG2 












PBMAL= 1.+08MAL 
CSTG2 = Gk2Cs7T 
CSTGI=G9CST 
CS7673G2CST 

K OTR y= LEW C0IILEN 47,33) 
KGTRY: KGTRY 2. 
DPTey: £,(LEN)> 












YES 


Depry «= f (LEW) 


Figure 3b 
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It is also used to discern between full hull ships and 
fine hull ships for the initial guess ot displacement. 
Stop execution: 
IFC INDEX.EQ.0) STOP 
This is the only place where the program is directed 
to stop execution. However, if a data card is left out, 
the program will stop when there is no more data cards 
to read. | 
Heading input: 


READ(8,123)(HEAD(I),I=1,20) 
123 FORMAT(20A4) 


Variable HEAD contains any alphanumeric information. 
It will be printed without change at the beginning of the 
output for all test cases using the data which follows. 
NOARM input: | 


READ(&,100) NOARM 
100 FORMAT(I2) 


NOARM is the total number of armament, aircraft 

and cargo payload data cards which follow. A maximum 
value of 20 is allowed. 

Armament, aircraft and cargo input: 

Rea 3, 101 CCRT TEn( J,2),J=1,5) ,cARGOW(1L),CARGOC(I), 
GR7WT(1I),GR/7CG(I),AMOWT(I),AMOCG(I),ACWT(I), 
ACCG(I),ELLD(1I),DHAR(I),HLAR(I),GR7CST(1I),I= 
1 ,NOARM) 

101 FORMAT(5A4,12F5.0) 

Up to 20 armament, aircraft and cargo payload descrip- 

tion cards are allowed. The data may be lumped together 


and appear on only one card or it may be broken up in any 


way desired. However, if a weight is given, a vertical 
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center of gravity must also be given on the same card 
immediately following the weight. The first 20 spaces 

on each card are reserved for an alphanumeric description 
of the payload items on the card. This is used as a line 
title in the output listing. Of the four weight categories 
on the card only the group 7 weights are included in the 
light ship weight. From a calculation standpoint it is 
unimportant whether ammunition or aircraft are listed 
separately or considered cargo. However, if they are 
listed as cargo they will be classified incorrectly in the 
output listing. 

All vertical center of gravity values should be 
inputed as a multiplier of the hull depth at amidships. 
This number is easier to estimate than the actual center 
in feet. AS an example, it a payload item has its center 
10 feet above the baseline and the depth at amidships 
is 20 feet the value entered is 0.5. 

The deckhouse area Should include only the area 
that 1S required to be in the deckhouse. If the area 
could be in either the deckhouse or the hull it should 
be included in the hull area. The group 7 cost should 
include only the cost ot materials and should be estimated 
for the time when the ship will be built. Because of 
format limitations the value input should be the true 
cost divided by 1000. All weights are in tons, all areas 
in square teet, all electrical loads in KW and all costs 


in dollars. This applies to all the input items. 





32 
Summing the payload data: 


WCARGO=0. 

WTG7= 0 e 

WTAMMO=0. 

BLOAD=0. 

AREADH=0 e 

AREAHL=0. 

CS TG/7=0 e 

DO Oi =i, NOAM 

WCARGO=WC ARGO+C ARGOW(TI ) 

WIG7=WTG7+GR7WT(1) 

WTAMMO={WTAMMO+AMOWT (TI ) 

WTAC=WTAC+ACWT(I ) 

BLOAD=BLOAD+ELLD(I ) 

GR7CST(I )=GR7CST(I )*1000. 

CSTG7=CSTG7+GR7CST(I ) 

AREADH=AREADH+DHAR (TI) 
10 AREAHL=AREAHIL+HLAR(T ) 


Only the totals or the payload items entered are 
needed tor calculations. These totals are computed in 
the above statements. 

NOELT inputs 


READ(&S,102) NOELT 
102 FORMAT(I2) 


This variable serves the same purpose as NOARM. A 
maximum value of 20 is allowed. 

pilectromres input: 

READ(%,103)((TITLE(J,I1+20),J=1,5),GR4WT(1I),GR4CG(I), 

ELD(1),DHA(I), HLA(I),GR4CST(I),1=1,NOELT) 
103 FORMAT(5A4,6F5.0) 

The comments made about armament input also apply 
here. Group 4 weights are included in the light ship 
weight. All items normally in group 4 should .be included 
Since no allowance is made in the program for navigation 


equipment. 
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Summing electronics input data: 

WTG4=0. 

DO 11 I=1,NOELT 

WTG4=WTG4++GR4WT(T ) 

ELOAD=ELOAD+ELD(T ) 

GR4¥CST(I )=GR4CST(I)*1000. 

CSTG4=CSTG4+GR4CST(T) 

ARE ADH=AREADH+DHA(I) 

11 AREAHL=AREAHL+HLA(I) 

The comments about summing armament data also apply 
here. Note that the area in the hull and deckhouse is 
a combined total ot all the payload inputs, both electronic 
and armament. 

storing initial costs: 


G4CST=CSTG4 
G7CST=CSTG7 


Because the value of variables CSTG4 and CSTG7 are 
modified during the execution of the program, a record 
is required of their initial values so that they will 
be available if another test case using the same payload 
7S Fun. 

Manning input: 

READ(&%,104) NOOFF,NOCPO,NOCREW 

104 FORMAT( 315) 

NOFF=FLOAT(NOOFF ) 

NC PO=FLOAT(NOCPO) 

NENL=FLOAT(NOCREW ) 

The number of officers, chief petty officers and 
crew members are input and then converted to real values 
so that they can be used in non-integer equations. 


Cost indices input: 


READ(8,110)DOLHR,G1LIND,G2IND,G3IND,GSIND,G6IND 
110 FORMAT(6F5.0) 
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The labor rate for the year the ship will be built 


must be estimated as must the inflation corrections for 

material costs. These corrections are for the period 

trom 1959 to the year the ship is to be built. Section 

3.14 gives some guidance as to the values to use. The 

indices can also include special material costs such as 

making the hull out of HNY-80 steel instead ot mild steel. 
Military mission consumables input: 


READ(%,105 )ENDDAY , WACFUL,CACFUL 
fO5 PORMAT( 3F5.0) 


The number of endurance days for dry provisions should 
be given. All other endurance days are taken as a suitable 
fraction of this number. The weight and center of gravity 
ot aircraft fuel must also be given. 

HPMS=0. 

This is done so that there is no difficulty in 
Sarling subroutine XECUTE. 

DOL ei mpuUL t 

READ(%,106) JOPT 

106 FORMAT(I2) 

IF( JOPT.EQ.2)GOTO 1 

This variable denotes the option to be used in entering 
vehicle performance data. If JOPT equals 1, the maximum 
Sustained speed is input. If JOPT equals 2, the maximum 
Sustained horsepower and machinery characteristics are input. 

Vehicle performance input for JOPT equal 1: 

READ(Y,107)VSUS,VEND,RGEND, PCE, PCM 


107 FORMAT(2F5.0,F10.0,2F5.0) 
GOTO 9 





The maximum sustained speed 1S input along with the 
endurance speed and endurance range. The propulsive 
coefficients at the endurance speed and at maximum speed 
must also be given. 

Vehicle performance input for JOPT equal 2: 

1 READ(S,108)HPNS, DHM,XLM,WTG2SFCMHP,SFCHHP,VEND,RGEND, 
PALOAD, PCE, PCh,CSTG2,WTLO 
moe FORMAT(7F5.0,F10.0,5F5.0) 

HPMS=HPMS #1000, 

CSTG2=CSTG2*1 000. 

GReesT=CSTG2 

This option is used primarily when a specific machinery 
plant is desired. Besides the maximum sustained horsepower, 
the length and minimum depth ot the engine room must be 
given. The total group 2 weight and specific fuel con- 
sumptions at maximum and halt power are required. The 
endurance speed and range must be given as before. Pro- 
pulsion auxiliaries electrical load and propulsive coef- 
ficients at full and endurance speed are required as are 
the group 2 material cost and lub oil weight. 

To fit into the format used, the horsepower and 
cost values must be divided by 100U. The cost of group 
2 must be stored for later use as were the costs of group 
4 and group 7. 

Options input: 


Y READ(Y,109 )NLYPE,FRSC,LEN,CP,CX,GMMUL,DBMAR ,DELCF 
moo FORMAT(15,OF5.0,F10.U) 


Machinery type must always be specified even if the 


option of inputing machinery characteristics was used. 


It machinery information was given, any number except 
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gero can be used for MTYPE. A value ofr zero will cause 

the program to assume the next input value is INDEX. 

This can be used to return and stop the program or to 

input a new set of data tor another ship. If the machinery 
particulars have not been specified, MTYPE equals zero 

has the same ettect. However, MTYPE equals 1 will cause 
the program to assume a. diesel plant is to be installed. 

A value ot MTYPE equal to 2 will result in calculations 

tor a CODOG plant. These are the only two machinery 

types tor which estimating relationships have been included 
in the program. If MY'YPE is any other number, the results 
are unpredictable. 

A tree surface correction can be specified. This 
value is interpreted as a rise in the full load KG of the 
ship in an amount equal to FRSC in feet. A negative value 
can also be used. In this case FRSC can be interpreted 
aS a lowering of the average KG value built into the 
program. 

The length between perpendiculars is given next. The 
Bvalues for Ch and Cy, have the standard naval architectural 
eenitivens. -- GliinUL if the fraction GM/Beam, A value 
ot U.1 is commonly used. The next variable, DBMAR, is 
a design and builders margin. This is given as a fraction, 
1.€., a 1U percent margin is input as 0.1. DELCF is the 
SC, value used in the calculation of horsepower. A value 
of 0.0004 is common. 


Pei rieee. 0) GOTO 2 
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One or more option input cards may follow the first. 
The program will assume that each one is a new input 
case with all other data the same. To input an entirely 
new set of data, an option input card with a value of 
ML£YPE equal to zero should be included. No other values 
need be given on the card. This will cause the next 
input card te be interzweted as a new value for INDEX. 

DBV.AR=1.+DBMAR 

CaeG2=Ex2CsT 

CSTG4H=GHCST 

CSTG7=G7CST 

DBMAR is changed to a multiplier for weight and the 
input costs are reset to their initial values. 


Initial KG estimate: 


KGTRY=LEN/(. O465*LEN+4.825) 
KGTRY=KGTRY-2. 


The initial KG value is estimated based on the values 
tor current Coast Guard vessels. These are shown in 
Figure 4. Reducing the initial estimate by 2 feet was 
found to be advantageous. It helps to insure that the 
tirst iteration will not exceed the KG constraints of 
the program. 

imeciay displacement estimate: 


Peiny=(.0005602-. 0011 59*LEN/1 00. )*LEN**3/35. 
ma (INDEX. EQ.2)DPIRY=. 00254 LEN*£39/ 35. 


Two initial values of displacement are available. 
For an INDEX value other than 2 the displacement is taken 
as that of conventional Coast Guard vessels. This 


displacement curve is shown on Figure 4. However, if the 
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vessel has a high speed-length ratio, greater than 1.3, 


it will not be possible to find the horsepower in some 
cases if the conventional ship curve is used. . This is 
due to the limitations on the Cy arrays. By specifying 
an INDEX equal to 2 the program assumes a fine tftorm at 
least for the first iteration. 
| Cubic number estiwate : 

GN=DPTRY*160 

An initial cubic number is assumed for use in estmating 
the electrical load on the first iteration. Values for 
Current Coast Guard vessels range from 140 to 180 times 
the displacement. 

CALL XECUTE(HPMS ) 

Control is transferred to subroutine XECUTE to 
Calculate the vessels characteristics. 


GOTO 9 
e2ND 


The GOTO statement causes the program to loop back 
and try to read another options input data card. The 


END statement stops compilation of the routine. 





3.2.4 List of Outputs of MAIN Subroutine 


MCCG (20) ccescvccves LC (EE) 


AMOCG(20)occceccveeeLC(EE) 
AMOWT (20 )ecsccceee « «LC (EE) 
AREADH. oe cccccccceeeLC( BB) 
000 oLC (3B) 
BORD. ccc ccceceeceeeLC(DD) 
oeeeeeLC(HH) 


AREAHLeccccccccce 


PWOPUL.....-0. 
CARGOC(20)..0000000eLC(EE) 
CARGOW(29)eccccccccs LC (EE) 
MOET. 5 5c oescecsceecLC(BB) 
BU oic cess soceee LOCC) 
BMID152) ee. eeccsee LC(TL) 
R972) ccc ccerceescLC(TT) 
BROL216).cccccccceesLC(T) 
BETO cece c ccc cee oLC CMM) 
DOG +e ccc ccc ccccce ee LC (NM) 
BERG]. oc cccscccce eo LO (he) 
MIs sis 6 0 eStore wee LOOGG) 
DBMAR...0ceccsccees cll (WW) 
EPIC... 0c cceceeeeLC(CC) 
TR ZO) occ ccccecceeelLC(KK) 
DHAR(20)...00. 5s. oeECiGies) 
MN Auer eels nein 2 oC USS) 
OER ccc ccc cece ccc eo LC(HN) 


PUN coc ccc cece s Gl (GG) 


ELD( 20). 50500 
ELLD(20). see. 
ELOAD. occcces 


RG Oistele o 6 siete se 


GLINDisccsece 
GZUNDa « <6 «0.66 


GSN Dis 0:6: s ie: oi 
GOWN s <'6.5: 6 ere 
GUMUL. sc0cese 
GR4CG(20). eee. 
GR4CST(20)...- 
GRUWT (20). 00. 
GR7CG(20)..e. 
GR7CST(20).e. 
GR7WT(20)..eee 
HEAD(20)..00. 
HUA( 20) eccece 
HLAR(20)ecseees 


HPHS ce cccccce 
JOPT . cccccee 
KGTRY.«ccces 
LENe cece ceee 


MIVIPE. <ccscee 
NCPO. cccccece 
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ooee eo LC (KK) 
pie oie ore te Kee 


coco e oe Lo( DOD 
ENDDAY. ccccces 


cece oC LAM) 
coe eeL0(GG) 


cocce esl GN) 
s% eee eel CHRD 
GBINDeccccces 


coeeeeLC(NN) 
ooeeeeLC(NN) 
oceee LC (NN) 
oceeeeLl (GG) 
ooeeeLC (EE) 
oceee eo LC(KK) 
ooeeee LC (EE) 
oceeeeLC (EE) 
oceceeLC (KK) 
ooee + LC (EB) 
cooeecelLC (FF) 
cocee el (KK) 
oe LC (KK) 


coccccve en DA 
coccecveme (CC) 
ccceeeeLC(GG) 
occceeeLC (BB) 


coceeelS(CC) 
se 6a'e LOA ) 
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List of Outputs of MAIN Subroutine (cont) 


Mis ccccccsccescesel (AA) 
RGARMs. ossecesccesesLC(FF) 
MOET... 0 cccceeecLC(FF) 
BE: os 6000005 e000eLC(AA) 
SMUOAD.. ..s0¢ec00eeeL0(DD) 
MME ooo cccecescecseLC(CC) 
BN fis... . «seeeiebo(ce) 
BMETID. ois «oss vecceeeLC(CC) 
MEM so coe sscenee elO(CC) 
BOOMER. oc cccs cece eo LC(CC) 
TITLE(5,40)..0ceeeeeLC(FF) 


* Output only if JOPT = 1 
## Cutput only if JOPT = 2 


VENDE os eeseccee se os LOKCe) 
WSUSUaaies or eerie LG (COD 
WACHULs «+ ccccee sees LO(Ley 
WCURG Oem rrauris cee LC CID) 
WORChnmesmeie tae LC CLE) 
WLAMHO.scccecccces elC(LL) 
WIGZE So ccccccce ces LOCWW) 
WIG. vo ccccccscces oIC(WW) 
WTE7. cc cccecsccecselLO(HW) 
GUI, Ohs, 0 sso eee CLO CEE) 


+t 


Midd os 6 « © 10 eels 6 eb CCoe 
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3.2.5 Nomenclature List 


ACCG(20) 
ACWT(20) 
AMOCG(20) 
AMOWT( 20) 
AREADH 
AREAHL 
AVGKW 

B 

BLOAD 
CACFUL 
CARGOC (20) 
CARGOW( 20) 
CCARGO 
CCST 
CGAC 
CGANMO 
CGCREW 
CGFUEL 
CGLO 
CGPE 

CGPS 

CN 

CP 
CR1(1152) 
CR2(972) 
CR3(216) 


aircraft vceg/Di0, D10 = depth amidships 
aircraft weight, tons 

ammunition vceg/D10 

ammunition weight, tons 

total input deckhouse arrangements area, sq ft 
total inzut hull arrangements area, sq ft 
average electrical load, KW 

maximum beam, feet 

total armament electrical load, KW 
aircraft fuel veg/D10 

cargo veg/D10 

cargo weight, tons 

total veg of cargo inputs, feet 
construction services cost, dollars 
total veg of aircraft, feet 

total vcg of ammunition, feet 

total veg of crew, feet 

total vcg of ships fuel, feet 

total veg of lub oil, feet 

total vcg of personal effects, feet 
total vcg of personnel stores, feet 
cubic number, cubic ft 

prismatic coefficient 

C, arrays for Taylor's Standard Series 


resistance estimate 
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Nomenclature List (cont) 


CSTGi 
CSTG2 
CSTG3 
CSTG4 
CSTG5 
CSTG6 
CSTG?7 
CWP 

CX 
DAVG 
DBMAR 
DOST 
DELCF 
DHA(20) 
DHAR (20) 
DH 
DHV 
DOLHR 


P MPIRY 


ELD(20) 
ELKW 
ELLD(20) 
ELOAD 
ENCVOL 
ENDDAY 
FRSC 


total cost of weight group i, dollars 
total cost of weight group 2, dollars 
total cost of weight group 3, dollars 
total cost of weight group 4, dollars 
total cost of weight group 5, dollars 
total cost of weight group 6, dollars 
total cost of weight group 7, dollars 
waterplane coefficient 

midship section coefficient 

average hull depth, feet 

design and builders margin 

design cost, dollars 

ACp used in resistance calculation 
electronics deckhouse area, sq ft 
armament deckhouse area, sq ft 
minimum machinery box depth, feet 
deckhouse volume, cu ft 

labor rate, dollars/hour 

displacement estimate, tons 
electronics electrical load, KW 
generator rated load, KW 

armament electrical load, KW 

total electronics electrical load, KW 
total volume of hull and deckhouse, cu ft 
number of endurance days for dry provisions 


free surface correction, feet 





Nomenclature List (cont) 


G1IND 
G2IND 
G3IND 
GECST 
G5IND 
G6IND 
G7CST 
GMMUL 
GR2CST 
GR4CG(20) 
GR4CST(20) 
GR4WT(20) 
GR7CG(20) 
GR7CST(20) 
GR7WT(20) 
H 

HEAD(20) 
HLA( 20) 
HLAR (20) 
HPMS 
INDEX 
JOPT 
KGTRY 

LEN 

LRD 

MCST 


cost index for weight group 1 

cost index for weight group 2 

cost index for weight group 3 

total material cost weight group 4, dollars 
cost index for weight group 5 

cost inuex for weight group 6 

total material cost weight group 7, dollars 
GM/Beam ratio 

total material cost weight group 2, dollars 
group 4 vcg/D10 

material cost weight group 4, dollars 

group 4 weight, tons 

group 7 vceg/D10 

material cost weight group 7, dollars 
group 7 weight, tons 

draft, feet 

descriptive heading (alphanumeric) 
electronics hull area, sq ft 

armament hull area, sq ft 

maximum sustained shaft horsepower 

control variable 

input option variable for machinery 
initial KG estimate, feet 

length between perpendiculars, feet 

raised deck length, feet 


miscellaneous costs, dollars 
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Nomenclature Liet (cont) 


MTYPE 
NC PO 
NENL 
NOARM 
NOELT 
NOFF 
NOCPO 
NOCREW 
NOOFF 
PALOAD 
PCE 
PCM 
RGEND 
SFCHHP 


SFCMHP 


SHPE 

SHPM 
TITLE(5,40) 
TOTCST 

VEND 

VSUS 

WACFUL 
WCARGO 
WFULLD 


machinery type option 

number of chief petty officers (real) 
number of enlisted men (real) 

number of armament input cards 

number of electronics input cards 

number or officers (real) 

number of chief petty officers (integer) 
number of enlisted men (integer) 

number of officers (integer) 

propulsion auxiliary electrical load, KW 
propulsive coefficient at endurance spsed 
propulsive coefficient at maximum speed 
endurance range, nautical miles 

specific fuel consumption at half horsepower, 
1bs/SHP-hr 

specific fuel consumption at maximum 
horsepower, lbs/SHP-hr 

shaft horsepower at endurance speed 
Shaft horsepower at maximum speed 

line titles.(alphanumeric) 

total lead ship cost, dollars 

endurance speed, knots 

maximum sustained speed, knots 

weight of aircraft fuel, tons 

total cargo weight, tons 

full load weight, tons 
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Nomenclature List (cont) 


WLSHIP 
WTAC 
WTAMMO 
WICREW 
WTFUEL 
WTGi 
WTGe 
WTG3 
wIch 
WTG5 
WTG6 
WTIG7 
WTLO 
WTPE 
WTPS 
XLbi 


light ship weight, tons 


total 
total 
total 
total 
total 
total 
total 
total 
total 
total 
total 


aircraft weight, tons 
ammunition weight, tons 
crew weight, tons 

ships fuel weight, tons 
grvup 1 weight, tons 
group 2 weight, tons 
group 3 weight, tons 
group 4 weight, tons 
group 5 weight, tons 
group 6 weight, tons 
group 7 weight, tons 


weight of lub oil, tons 


weight of personal effects, tons 


weight of personnel stores, tons 


machinery box length, feet 





4 
3.3 Subroutine XECUTE 


3.3.1 Introduction 

Subroutine XECUTE controls the execution of each 
trial case. All of the remaining subroutines except 
two are called by this routine. Error messages are also 
printed by this routine whenever a constraint of the 
program is violated. 

In three cases, iteration control statements have 
been included in XECUTE. The first set of these statements 
controls the convergence of the spsed-power estimate 
when maximum sustained horsepower is specified. A modified 


(2) is used with the slope ap- 


Newton-Raphson technique 
proximated by a secant to the speed~power curve. 

The second set of iteration control statements 
controls the convergence between estimated displacement 
and full load weight. A similar modified Newton-Raphson 
technique is used. 

The last set deals with the convergence of the 
vertical center of gravity estimate. For this set there 
is a tendency for the values to oscillate and a Newton- 
Raphson aoa could not be used. The average of the 
estimated and calculated values of KG is used instead. 
While this eliminates the oscillations, convergence may 
be slow in some cases. 


A flow chart for this subroutine is shown in Figures 
5a, Sbi, Dc and 5d. 
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XECUTE SUBROUTINE 
FLOW CHART 


REAL AND 


LABELED oe 
START XECUTE COMMON ee 5 ; 
LNDEX:= ! 
STATEMENTS 


YES 
> NO > , CALL UWD/M(p) 
NO if 


ES 
CALL OUTPUT(Z,4,1) 
CALL H PCALC(YSUS, 
Sine) 
ZA 


WRITE 
VCG ERROR 
MESSAGE 


EXCKG =9. 






















WwW RITE 


PROGRAM PkocEgy 
To MELT MFU Tt 


CASE 





YES 
> 









WRITE 
LIMUTS OF 






STOLED 
RESISTANCE DATA 
EXCEEDED 


SHPM=SHPAB 


CAL OUTPUT (YG, 1) 





WRITE 
HP ERROR 
MESSAGE 


Figure 5a 





CALL HeCALCC VMAX, YES [<¢o TO 
SHP,Pem ) = 9908 






SHPI> SHP 
UMAX = .95* 





CALL HPCALC CVMAKL, 
SHP, Poem) 


SUP 22 SNE 
SRPL 
DELHP> Saaidnax! bey 
VMAX = VMAX oes 
eA oe 
; = VM 
eS ie CALL HPCALCC VE, 
4 VMAX I * HePM= HPMS | SHP PCE 
LIS Ss 
No 
VMAK = VMAL/ Ye 
S 
SHPE = SHP e 
NO 


pe ss GOTO 
ALL EPLAUT 1000 


CALL MACHLA 









[Vaarti- vA 









Toe LARGE 
A VOLUME 
REQUIRED 








a MNBSIzEB) YES 
( cain volume ) ( cain volume ) CALL OUTPUT(G,61>) 


Figure 5b 


CALL WEIGHT 


CALL VCG(CFULLD) LNDEX= | 


FACTOR=(0°rrRy- 
WFULLO - OPH wri) / 
(DPTRY~DP1) 
DPl=z DPprevy 
WT12z WFULLD 







YES 





YES No CALL OUTPUTS, } 1) 


LMNoEXx 
LNoEX+ | 





YES 






KGTRY= 
CCPULLD + EGTRY)/2 





JM DEX= INDEX| 












WRITE 
NO K& BALAKME 
Possi@le (NV 

20 ITERATIONS 





EXxCkKG = 
No KGTRY - 
CAXLD 


R< 0.001 


Figure 5c 






OPTRY= DPray~(DPTRY= 
WFULD) /eacrok 
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WTt = WFULLOD 
DP!+ ODPTRY 
DPTRY= WFULLD 










Eo 
Yes 


CALL OUTPUT LE G,1) 










WRITE 
No WEIGHT 

Balance PossiBle— 
IN 20 rTERATION 






Go TO 
1OoooO 





CALL COST(LEN) 


CALL SEASPD(AVSA 


CALL OVTFUT 
CAVSP, ExCKG,2) 


RETVUR»S 


Figure 5d 
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303-2 Input List for Subroutine XECUTE ™ 
DPTRYecceccccces el (GG) 
HPMS cccccccccccces eOUDA Only if JOPT = 2 
JOPT so secvcccese el (CC) 
KGTRY.cccccececeeLO(GG) 
Mis 65s + o0ccee eset ( BB} 
MEUE Ds esse eee ee el: ( CC) 
PCE sc cecccessesevelLC(CC) 
OME ct sc cle + +e sislee lO tec ) 
WE cess seccceecee LC (CC) 
MW Scecccreceses sLCloc) Only if JOPT = 1 
WFULLD. cc cccccceeLC (WW) 


3.3.3 Statement Desoriptions 

REAL and COMMON statements will not be listed with 
individual subroutines but may be found in Appendix C. 

EXCKG=0,. 

JNDEX=1 

LNDEX=s1 

These statements initialize values for the excess 
KG variable, the VCG loop counter and the weight loop 
counter respectively. 

7 CALL UWDIM(R) 

This subroutine is called to calculate the beam 
and draft of the ship. The variable R transfers the 
value of any excess KG. 

IF(LEN.LT.0O.) GOTO 999 

If a beam to draft ratio of greater than 4.0 is 
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required to satisfy the stability requirements, LEN is 


set equal to -1000. in UWDIM. This signals XECUTE to 
print an error message and return control to MAIN. 

IF(JOPT.EQ.2) GOTO 1 

The next group of statements control the speed- 
power estimate. Two procedures are used depending on 
which performance inpuw option is specified. 

CALL HPCALC(VSUS ,SHP, PCM) 

IF(LEN.LT.O.) GOTO 998 

GOTO 2 

These statements are used if the sustained speed 
is specified. Subroutine HPCALC is desiged to return 
avalue of shaft horsepower, SHP, if speed and propulsive 
coefficient are specified. In this case VSUS is the 
speed and PCM the propulsive coefficient used. 

LEN is used here also as a signal that an error 
message is required. A negative length, LEN, will occur 
if the limits of the C, array are exceeded. The ranges 
of the Cy, array ares 

048 < Cy S +70 
2,.<B/7 <4. 
for .001< Cy <.003 
5 < V/JL < 2. 
for .003< C, <.006 
5 < W/L <1.3 

The GOTO 2 statement causes the routine to skip 

the speed estimate used when JOPT = 2. 





1 VMAX = (HPMS*200./LEN)**. 33333 
If the maximum sustained horsepower, HPMS, was 
specified, the following statements are used for estimating 
the maximum sustained speed. The statement given above 
provides a rough first guess at the speed. 
& CALL HPCALC(VMAX,SHP, PCM) 
IF(LEN.LT.0.) GOTO 998 
SHPi=SHP 
VMAX1=.95*VHAX 
CALL HPCALC(VMAX1,SHP, PCM) 
IF(LEN.LT.O) GOTO 998 
SHP2=sSHP 
Two points on the speed-power curve are calculated 
in the above statements, one point at VMAX and one point 
at .95 x VMAX. The constant, .95, was chosen arbitrarily 
to give a secant to the curve which has a slope close to 
the actual slope at VMAX. The same limits of the Cy, array 
apply here. 
DELHP=(SHP1-SHP2)/(VMAX-VMAX1 ) 


This is the approximate value of 


_G{SHR) at VMAX 


dv 

VMAX1=VMAX-(SHP1-.8*HPMS )/DELHP 

By definition the maximum sustained speed is the 
speed that can be made when the shaft horsepower is 80 
percent of the maximum sustained horsepower. This assumes 
the ship has a clean bottom and is operating in calm 
water. Therefore, the function to be solved iss 

£(V) = SHP(V) - .8 x HPMS = 0 

The Newton-Raphson method is an iterative procedure 
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in whichs 


£(Vk) 
Wrsi “"%e ~ F4(¥,) 


where f' (Vy) = a 2G 


ae aly 


in the present case: 

Vy, = VMAX 

Vigd = VMAX1 

f£(V,) = SHP1 - .8 x HPMS 

£*(V,) = DELHP 

DELHP is only an approximation to f'(V;)3; however, 
it produces very good convergence. 

IF (ABS (VMAX1-VMAX).LE.VMAX1#.005) GOTO 3 

The iteration stops when the difference between 
any two successive speed estimates is less than half a 
percent of the best estimate. The above iteration technique 
is very fast and convergence to these limits usually 
occurs in 3 to & iterations. 


VMAX = VMAX1 
GOTO 4 


These statements are used if the two speed estimates 
. are not within the above limit. The new velocity estimate 
is used and another iteration is run. 
3 VSUS = VMAX1 

SHPM = HPMS 

GOTO 5 

The above pair of values are saved for output and for 
use in other subroutines. 


2 SHPM = SHP/,8 


In the case where VSUS is specified, the value SHP 
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is equal to SHPM times .8 as before. 


IF (MTYPE.EQ.1.AND.SHPM.GT.14000.) GOTO 997 

If a diesel plant is specified, horsepowers of 
greater than 14000 are outside the range of validity 
of the model. An error message is printed if this occurs, 

5 CALL HPCALC(VE,SHP, PCE) 

IF(LEN.LT.O.) GOTO 998 

SHPE=SHP 

The endurance horsepower is calculated without regard 
to fouling or sea conditions. This statement is executed 
for both performance options. 

CALL EPLANT 

This subroutine is called to calculate the average 
electrical load at sea and to size the ships generators. 

CALL HMACHLQ 

This subroutine is called to calculate the weight 
of fuel, lub oi11 and potable water required. 


IF(JOPT.EQ.2) GOTO 6 
CALL MBSIZE(B) 


Subroutine MBSIZE calculates the length and mininun 
depth of the machinery box. If JOPT = 2 these values 
are given as input. 


6 CALL VOLUME 
IF(LEN.LT.0O) GOTO 996 


subroutine VOLUME sets the sheer line and balances 
the required area and available area. If insufficient 
area is available, the value of LEN will be negative. 
CALL WEIGHT 


All the weights not previously estimated or inputed 
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are calculated in this routine. The primary value returned 
is WFULLD, the full load weight. 

IF(ABS(WFULLD-DPTRY).LT.2.) GOTO 8 

The calculated full load weight is compared to the 
displacement estimate made earlier. If the two do not 
agree within 2 tons, a new estimate of displacement is 
made and the vrogram rc curns to call UWDIM. If there is 
agreement, the program proceeds to calculate vertical 
center of gravity. 

IF(LNDEX.GE.2) GOTO 9 

WT1=WFULLD 

DP1=DPTRY 

DPTRYeWFULLD 

GOTO 10 

Two values of weight and two values of displacement 
are needed for the Newton-Raphson procedure explained 
below to work. On the first iteration, the above simplified 
procedure is used. 

9 FACTOR=(DPTRY-WFULLD-DPi+WT1 )/( DPTRY-DP1 ) 

DP1i=DPTRY 

WIlL=WFULLD 

DPTRY=DPTRY-(DPTRY-WFULLD)/FACTOR 

In this cases 

f£(DPTRY) = DPTRY - WFULLD(DPTRY) = 0 
Again a secant approximation is made to f°(DPTRY). This 


is the variable FACTOR. 
= . £(pPTRy 
DPTRY;4, ™ DPTRY« - Serpoaay 


10 IF(LNDEX.GT.20) GOTO 995 


In case the iteration does not converge, this statement 
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provides a safety valve to stop the loop. An error 


message is printed. 


LNDEX=LNDEX+1 
GOTO 7 


A new iteration is run beginning with a call to UWDIHN. 
8 CONTINUE 

LNDEX=1 

CALL VCG(CFULLD) 

Once the full load weight is sufficiently close to 
displacement, the vrogram proceeds to calculate KG. The 
weight counter is reset. The value CFULLD is the calculated 
KG. This value is compared with the previous estimate, 
KGTRY. 

IF(ABS(CFULLD-KGTRY).LT..1) GOTO 11 

This iteration stops when the two KG valves are less 
than 0.1 feet apart. 

IF( JNDEX.GT.21) GOTO 994 

Again a safety valve is provided in case the iteration 
does not converge. 

IF(CFULLD.LT.KGTRY) GOTO 12 

13 KGTRY=(CFULLD+KGTRY)/2. 

JNDEX#2JNDEX+1 

GOTO 7 

me Pe(R.LT..OOL) GOTO 13 

EXCKGsaKGTRY-CFULLD 

In some cases the minimum beam to draft ratio may 
provide more stability than is required to satisfy the 
GM criterion. The ship is still feasible if this occurs, 
but may be over designed. When this occurs R will have 
@value greater than 0. In fact, it will be equal to the 


excess KG when KGTRY is used in UWDIM. The input value 
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of KGTRY is modified in UWDIM by adding R. Therefore, 


the value of KGTRY at this point in the program is the 
maximum value of vertical center of gravity for which 
a minimum beam to draft ratio will still be valid. If 
CFULLD is less than this value there is no point in doing 
another iteration. Instead, the program continues and 
eventually prints the cxcess KG available as an output. 
In all other cases a new estimate of KGTRY must be 
made. AS mentioned earlier, this estimate is taken as 
the average of the initial estimate and final calculated 
value. This is done to prevent the two values fron 
oscillating. As an example, in the case of the 378° 
WHEC, a KGTRY of 14 feet returned a CFULLD of 16 feet. 
A KGTRY of 16 feet returned a CFULLD of 14 feet. This 
oscillation may not occur in other cases, but the averaging 
estimate used will converge in every case, although pos- 
Sibly at a slow rate. 


11 CONTINUE 
CALL COST(LEN) 


Subroutine COST calculates any remaining cost items 
not given as input. 
CALL SEASPD(AVSP) 
This subroutine estimates an average speed, AVSP, 
in the North Atlantic Ocean which includes the effects 
of sea state. 
CALL OUTPUT( AVSP,EXCKG,2) 
Subroutine OUTPUT contains the write statements for printing 


all the output and input data. The value 2 causes OUTPUT to 
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print both the input and output data. 


RETURN 
The present input case has been completed and control 
is returned to the MAIN program for more data. 
The error message statements are all similar in format 
and generally self explanatory. Only one typical one 
will be discussed. 
994 CALL OUTPUT(0.,0.,1) 
Subroutine OUTPUT is called with a value of 1. This 
directs OUTPUT to print only the input data. 
WRITE( 5,106) 
106 FORMAT(///' NO BALANCE BETWEEN ASSUMED AND CALC KG 
WITHIN 0.1 FEET COULD BE MADE IN 20 ITERATIONS'///) 
GOTO 1000 
1000 WRITE(5,102) 
102 FORMAT(*PROGRAM PROCEEDING TO NEXT INPUT CASE'///) 
RETURN 
END 
The program prints the appropriate error message 
and returns for more data. This concludes subroutine 


XECUTE. 


3.3.4 Output List for Subroutine XECUTE 
DPTRY ecccescce ell (GG) 
RRCKGs oo eve veces eDODA 
WETRY cco ccceeees LC (GG) 

SHFE «sc scseeeseeLC(CC) 
SHPM es ccccecee elS(CC) 
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3.3.5 Nomenclature List 


All variables have the same definition as given in 


the MAIN program nomenclature list except for the following: 


AVSP 
CFULLD 
DELHP 
DP1 
EXCKG 
FACTOR 
JNDEX 
LNDEX 


SFCH 
SFCM 
So HP 
SHP1 
SHP2 


_ VMAX 
VMAX1 
wrt 


average speed in N.A. Ocean, knots 
vertical center of gravity, KG, feet 
derivative of horsepower function 
used to 2tore DPTRY value, tons 
excess KG, feet 

derivative of weight function 

veg counter 

wejeht counter 

excess KG, feet 

same as RGEND in MAIN 

same as SFCHHP in MAIN 

same as SFCMHP in BAIN 

dummy used for shaft horsepower 
dummy used for shaft horsepower 
dummy used for shaft horsepower 
same as VEND in MAIN 

dummy used for speed, knots 

dummy used for speed, knots 

used to store WFULLD, tons 
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3.4 Subroutine UWDIM 


364.1 Introduction 

Three of the important variables determining the 
underwater hull dimensions are given as input. These are 
length, Cy and Cy,» Two others, the beam and draft, are 
Calculated in UWDIM. The waterplane coefficient, Cwps 
is also calculated, but only as a linear function of 
Coe 

Displacement is an input to this subroutine and 
therefore the value of beam times draft is fixed sinces 

Bx T = Displacement x 35./(Lpp x Cy x Cx) Ref. (3) 

Reference (3) may be referred to for a detailed 
definition of the naval architectural variables used 
in this routine. 

The beam to draft ratio, B/T, is determined by the 
GM/B criterion. The GM available is given by the formulas 

GM = KB + BM - KG - Free surface correction Ref.(3) 
Values for KG, free surface correction and GU/B are 
inputs to this routine. A value for KB is estimated 
uSing Morrish‘s approximate femmes?” 

KB = T - 1/3 x (T/2 + V/A) 

where T is the draft in feet 

Vis the volume of displacement = 


B scerex Lop x C. x Cy Gul son 


p 
A is the waterplane area = B x Lpyp * Cwp» 84 ft 
Substituting for V and A gives: 


KB = T x (5/6 = (Cy x Cy)/(3 x Cyp)) 
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A value for BN is estimated using: 


3 
py = Ebp x BO x Cy 


where C, is a transverse moment of inertia coefficient 
which is assumed to be a linear function of Cywp. 

A number of coefficients are used in this routine 
which are defined by the above equations. These are 


listed below. 


C1 = DPTRY*35./(LEN*CP*CX) =» Bx T 
C2 = .833 - CP*CX/(3*CWP) = KB/T 

C3 = LEN*CALPH/(DPTRY*35.) = BM/B? 
C4 = KG ve FRSC 

C5 = GM/B = GMMULY 


Figure 7 gives a typical plot of the equations 


KB + BM - GMa C2xT+03xCl1° -¢1 x5 
“varae 7 


The curve shown was calculated using data for a 
378° WHEC. 

The value of KG is also plotted in Figure 7. In 
the discussion below, KG refers to the value of KG + FRSC. 

A variable R is now introduced with a value of: 

R = KB + BM - GM - KG... 
When R = 0, the stability criterion that GH = C5 x B is 
just satisfied. This is the desired solution. This 
point is shown as point 3 in Figure 7. 

However, the solution is only valid provided point 3 
lies within satisfactory limits of beam to draft ratio. 


These limits are specified by the speed-power estimation 
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It can be seen from Figure 7 that if the value of 
KG is lower than the value of KB + BH - GH at point 1 
there would be no draft for which R #2 0. The minimun 
GM that could be provided would still be greater than 
o> x B. 

The solution method used in UWDIM proceeds as follows: 

(1) The value of B/T at point 1 -is calculated. 
Then a check is made to see if point 1 lies between the 
acceptable limits of B/T. 

(2) If point 1 lies within the limits of B/T, the 
value of R is calculated. If R is positive, no solution 
exists to R = 0 but the ship may still = tral In 


fe 
/¥ mony )4 


this case the minimum GM is greater than us x Be When -—< | 
this occurs the KG estimate is raised until R = 0 at 
point 1 and the program returns to subroutine XECUTS. 
If it is found later that the actual KG of the ship is .v: 
lower than the new KG estimate, ths difference is printed 
as an excess KG. 

If the value of R at point 1 is negative, the program 
executes the procedure described in step 4. 

(3) If point 1 does not lie within ecceptable 
limits of B/T, the value of R at point 2 is calculated. 


If R is positive at this point the procedure is the same 
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as that outlined in step 2. If R is negative, execution 


proceeds to step 4, 

(4+) If the value of KG is greater than the value 
of KB + BM - GM at point 1 or point 2, the value of R 
at point 4 is calculated. If the value of R is negative 
at this point, no draft can be specified within the 
limits of B/T that will satisfy the stability criterion 
and the program will have to go on to another input case. 

If R is positive at point 4, the solution must lie 
between point 4 and point 1 or point 2, whichever is 
lower. It is assumed that if point 1 is within the 
limits on B/T that the solution does not lie between 
point i and point 2. 

(5) To solve for the draft at which R = 0, a synthetic 
division technique is used together with the Newton- 
Raphson iteration technique. This is sometimes known as 
the Birge-Vieta eae The equation for R must be 
modified slightly to use this technique. 

£(T) = TR/C, = T - Cyt /Cy -CyCeT*/Co + C5C42/Cp = 0 

In the Birge-Vieta method three columns are used as 


shown belows 


The values in the first column are the coefficients 


in the equation for f(T), namely: 


ay = -Cy/Co, Aap = -CyCc/Co, a, = 0, ay = 0304°7/C, 
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The values in columns two and three are calculated 


using the equations: 
by zag+zxdb,, 
Cc; = by + ZX C3_y 
where z2 is the current estimate of the draft, T. 
The new estimate of draft is then: 
Zaow Z 2 - dy/e- 
New values are computed for columns two and three using 
Znewe The procedure starts with the draft at point 4 


and continues until the difference between z w and z 


ne 

js sufficiently small. This iteration procedure is very 

fast and usually only three or four iterations are required. 
A flow chart for this subroutine is shown in FPigures 


6a and 6D. 


3.4.2 Input List for Subroutine UWDIN 
Che eveccceseeelLC(GG) 
Gees ceosce ee ele (GG) 
CP ineuewre «000 oo « LC (GG) 
Ciel. « < vieleiele « « « LOGG) 
Wms cc cess sel (GG) 
FRSC eeccees eee LC(GG) 
[omMels cs cece oe ok (Ee) 
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UWDIM SUBROUTINE 


FLOW CHART 


DPTRY* DISPL/ 1.0/4 

CWP: .806S ¥CP4+. 2645 

| CAL PH of 2#CWP ~0 04/4 
Cl OPTeY#35/LEwwce Hc) 
C2Z= +833 -CRECX L340 CWP) 


REAL, DIMENS/OA 


AND LABELED 
START UWDIM eG 


STATEMENTS 








HMAX= SQRT(C!Z.005) C42 C4 +FRSC 














Bmins= Saer((cs + 
<2, 
Bmin= Ci/#max BTR < 2.09 SQRT (CS#42 4) 2NCIMCZACD 
Re C2#HMAX + OR 7 (€6%C3)) 


CZRBMINENZ ~ CF BTL > 4.00 Hains C1/Bmiv 
-—C5#*#GmiV BTR: Bm iW /Hot 








R= CU* Hm) 
+ C3" BMINAKS 


R< A ~C4-CS¥8miN 


NO 


C42 C4-FESC HR 


C4= C4-FRSCt+R 
Bl = BMI) 


NO | Brs amin 
H s+ HMIAN 


Hz HMAX 





| eS 


HMins = SarrdAss 
Bmaxea Clam in) 
Re: Cok Hm 
+€3 » BMAX KHZ 
—-C4- CS BMAX 


RETURN RETUBZN 


ACi=/ 


£0 AWG 
Al2): = 


A(a: - Cie 


Al4y=d 
LEN = -/666. 


Ats)= Cot 
c2 
B(i= | 
Ca)=] 
¢>Hmian 





RETURKAI 


END 


Figure 6a 












FoR 12,5 
Bie AZ) t2*G(r-1) 
C(p2a) t24c(r-N 


—_— — ce — 


BCS. 












[2NEW-z2] 


Se ZNEW < .Ofk 2VEW 


H=SNBW 
Biz CI/H 


C42C9 - FRSC. 
R=, 





Figure 


6b 


RETURN 
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3.4.3 Statement Descriptions 


DPTRY = DISPL/1.014 

DISPL has the value assigned to DPTRY in other 
subroutines. It is assumed that the total displacement 
is 1.4 percent greater than the molded displacement. 
DPTRY has the value of molded displacement in this sub- 
routine. 

CWP=.8065*C P+. 2645 

This relationshin is plotted in Figure 8. 

CALPH=.12*CWP-. 0414 

This relationship is plotted as Figure 8-8(b) in 
Reference (1). This relationship applies to destroyer 
ship types but was not checked against Coast Guard designs. 

CisDPTRY*35./(LEN*®CP#CX) 

C2=.833-CP*CX/( 3. *CWP) 

C3aLEN*CALPH/(DPTRY*35. ) 

CuxC4+FRSC 

The derivation of these coefficients is described 
in the introduction to this section, 

BMIN=SQRT( (C5+SQRT(C5#*2+12. *C1*C2*C3) )/(6.#C3) ) 

HMIN=C1/BMIN 

BTR<BMIN/HMIN 

IF(BTR.LT.2..0R.BIR.GT.4%.) GOTO ! 

The ratio of B/T at point 1 in figure 7 is calculated 
using these statements. The formula for BMIN is derived 
by setting dR/dE = 0 and solving for B. 
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R=aC 2*HMIN+C 3* BMIN**3-C4-C 5*BMIN 
HF(R.LE.O) GOTO 2 
C4eC4-FRSC+R 
BI=BiMIN 
HsHMIN 
RETURN 


These statements compute the value of R if point 
1 lies within the acceptable B/T limits. If R is positive, 
the program returns the value of beam and draft at point 
1 and increases the estimate of KG by R. 


1 HMAX=sSORT(C1/2.005) 
BMIN=Ci/HMAX 
R=C 2*#HMAX+C 3*BMIN*®#3-C4-C5*BMIN 
IF(R.LE.O) GOTO 2 
CheCh-FRSC+R 
Bi=BMIN 
H=HMAX 
RETURN 


These values are calculated for point 2 in Figure 7. 
Actually they are slightly off of point 2 to insure an 
acceptable value in the speed-power routine. 


2 HMIN=SQRT(C1/3.995) 
BMAX=C01/HMIN 
Re 2*HMIN+C 3*BMAX**3-C4-C 5*BHAX 
IP(R.DT.O.) GOTO 6 


These values are calculated for point 4 in Figure 7. 
If R is negative at this point, an error message must 
be printed. 


A(1)=1. 
A(2)=-C4/C2 
A(3)=-C1*#C5/C2 
A(4)=0. 
A(5)=C3/C2*C1**3 
B(1)=1. 
G(1.)dh.. 
SHIN 
3 DO Me i= 250 
B(I )wA(I )+Z*B(I-1 
4 Mer meer kan tat 
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The columns of the Birge-Vieta method are calculated 
in these equations. 
ZNEW=Z-B(5)/C (4) 
IF(ABS(ZNEW-Z).LE..OL*ZNEW) GOTO 5 
2=ZNEW 
GOTO 3 
A revised estimate of draft is calculated and the 
iteration is repeated unless the difference between the 
old and new values is less than 1 percent of the new 
estimate. 
5 H=ZNEW 
Bi=C1/H 
C4=C4-FRSC 
RQ e 
RETURN 
The results are returned for point 3 in Figure 7. 
6 LEN = -1000. 
RETURN 
END 
These statements signal subroutine XECUTE to print 


an error message and proceed with a new data case. 


3.4.4 Output List for Subroutine UWDIM 
Blacccccccsvees el (GG) 
BW ss ve ve <aitiewebe(Ge) 
He esccvesccvece es LC(GG) 


3.4.5 Nomenclature List 
All variables have the same definition as given 
in the MAIN program nomenclature except for the followings 
A(5) celumn 1 of the Birge-Vieta method 
B(5) column 2 of the Birge-Vieta method 





ae 
oe 
a 
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Nomenclature List (cont) 


Bi 
BMAX 
BMIN 
BIR 
C(5) 
C1 

C2 

C3 

cl 

> 
CALPH 
DISPL 
HMAX 
HMIN 


ZNEW 


Beam, same as B in MAIN 
dummy value of beam 
dummy value of beam 
beam to draft ratio 
column 3 in of the Birge-Vieta method 
coefficisnt, see text 
coefficient, see text 
coefficient, see text 

KG or KG + FRSC, feet 
GMMUL in MAIN 
coefficient, see text 
DPTRY in MAIN 

dummy value of draft 
dummy value of draft 
index variable 

KB + BHM - GM = KG - FRSC 
dummy value of draft 
dummy value of draft 
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3.5 Subroutine HPCALC 


3.51 Introduction 
Subroutine HPCALC performs a Taylor's Standard 
(7) 


Series power estimation. The procedure used is very 
similar to the hand calculating procedure. However, two 
problems are of interest with regard to doing the calculations 
on a computer. First, the C, curves plotted in Reference 
(7) must be converted to array form and stored in the 
program. Subroutine CRVAL, described in Section 3.6, 
chooses the correct value of C, from this stored array. 
Second, the surface area of the hull must be estimated. 

One approximate formula for wetted surface is: 

S = 1.7°L‘T + v/a!) 

where L is the ship length 

T is ths draft 
VY is the volume of displacement 

This form was used as a starting point and an empirical 
relationship was derived to bring the wetted surface 
more closely in line with that specified on the curves 
' given in Reference (7). To do this a number of points 
were chosen on the curves given. These points were 
compared to the above formula and some correcting factors 
for B/T, Cy and C, were added to reduce the difference. 
The formula given below gives values within 5 percent 
of the values calculated from the curves of Reference (7) 
over the entire range. 


Sm (1.7+L*T+¥/T)(.0053(B/T)*~.02(B/T)+3*Cy+. 08°C p+. 926) 





This formula is probably not the best method for 
estimating the surface area. One method would be to 
convert the curves to arrays and interpolate to find the 
correct area. Such a method requires much more time and 
compute storage to produce a result only slightly better 
than this formula. Most of the error between a power 
Calculated by this routine and by hand is probably due 
to error in the surfece area prediction. 


A flow chart for this routine is given in Figure 9. 


3.5.2 Input List for HPCALC Subroutine 
Beccccvsccscee ell (GG) 
Biles ese caeliee se LO UGG) 
DELCF. csc ceseeclLC(CC) 
DPTRY.eseeeeee eLC(GG) 
He wssccccevcce el (GG) 
Eels occ cee el BB) 
Po ccecccccccesee eDCDA 
VK awe cccccevscee eDCDA 


3-5-3 Statement Descriptions 

DISPL=DPTRY/1.014 

DISPL is the molded displacement. It is assumed 
that the total displacement is 1.4 percent greater than 
the molded displacement, 

VOL=DISPL*35. 

BDR=B/H 

RNs VK*LEN #1 32050, 

CFn,075/(ALOG10(RN)-2. )*#2+DELCF 


SLRAT=VK/(SQRT( LEN) ) 
CV#VOL/( LEN) ##3 
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HPCALC SUBROUTINE 
FLOW CHART 





DISPL= DPTRY/L014 
VOL = DISPLISS 
DR: B/H 
Bo . VKeLENY 32050: 
x Sieg : P 
slrare /SakVUEY) 
CV= VOolL/tEw 






REAL AND 


LABELED 
PCALC 
START HPCAL Boe 


STATEMENTS 

























v e TS 5S i D <, 
CV>,003 ie CV <.00! ne SLRA No cP< 48 ae BDR<Z 
AND OR OR a oO on 
SLRAT>/.3 CV?.Cob6 SLRAT ? 2, CP >.7 BOR7 FY 
NO 
YES 






CT: CErCQRw.00/ 
S=CL7H ENE 4VOL/A) © 
C0053 &oe*- .02 BOR+3CV 
+ .O8 CPs ,926) 

Az, 0CO4IGHS.YIOSHS 

EXPBar= Aaya CT n 
CHPRPP = (IME 7-2 SRR ORE 
EXP= o WE HP APP 







CALL CRLYAL 
C BOR, CV, CP SLRAT, CR) 







ExP= (1.94 - 6+ SURAT) 
x EUPAPP 


YES 
EHP = 1.27 SYPAPP S 


SHBS CEHP/PC 








Figure 9 





| 78 
These formulas calculate, respectively, volume of 


displacement, beam to draft ratio, Reynold‘’s number, 
Cr, speed length ratio, and C,. Cp, is the coefficient 
of frictional resistance calculated based on the ITTC line. 

IF(BDR.LT.2..0R.BDR.GT.4.) GOTO 1000 

IF(CP.L7T..48.0R.CP.GT..7) GOTO 1000 

IF(SLRAT.LT..5.0R.SLRAT.GT.2.) GOTO 1000 

IF(CV.LT..001.0R.CV.GT..006) GOTO 1000 

IF(CV.GT. .003-ANDeSLRAT.GT.1.3) GOTO 1000 

These statements check to make sure none of the 
constraints of the C,. array are violated. If any is 
violated, the program will write an error message and 
go to the next input case, 

CALL CRVAL(BDR,CV,CP,SLRAT,CR) 

This subroutine is called to determine the value 
CR, the residuary resistance coefficient. The routine 
is described in section 3.6. 

CTaCP+CR*, 001 

Sa(1.7*LEN*H+VOL?7H) #*(.0053*BDR*BDR-. 02*BDR+3. *“CV+ 

o 08*CP+.926) 

Az. 004381 e 990535 

EHPBHeA*VK**3*CT 

The total resistance coefficient, wetted surface 
- and finally the bare hull effective horsepower are cal- 
Culated in these statements. 

EHPAPPe(1.467-.2*SLRAT) *EHPBH 

The effective horsepower with appendages has been 
based on data for the 378° WHEC. This data is shown 
in Figure 10 together with the above approximation. 

EHP=.9*EHPAPP 

IF(SLRAT.LT.1 e+. AND.SLRAT.GT. e 9 /EHPe(1i o74-.6*SLRAT) * 


EHPAPP 
IF(SLRAT.LE..9) EHP21.2*EHPAPP 
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These equations approximate the EHP/EHPrayljor curve 


for the 378° WHEC as shown in Figure 11. 


SHP=EHP/PC 
RETURN 


Finally, the shaft horsepower is calculated using 
the input propulsive coefficient. 
1000 LENs -1C00. 

RETURN 

END 

If a constraint of the C, array is violated the 
value of LEN is made negative to signal subroutine XECUTE 


to write an error message. 


3.54% Output for Subroutine HPCALC 
The only output is SHFS, which is returned as a 
call statement dummy argument. 


3.5.5 Nomenclature List 
All variables have the same definition as given in 


the MAIN program nomenclature except for the following: 


A dummy coefficient 
BDR beam to draft ratio 
CP frictional resistance coefficient 
CR residuary resistance coefficient 
CT total resistance coefficient 
CY Taylor's volumetric coefficient 
DISPL molded displacement 
EHP effective horsepower 
EHPAPP effective horsepower with appendages 


EHPBH effective horsepower bare hull 
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Nomenclature List (cont) 


PC 

RE 

RN 

S 
SFCH 
SFCM © 
SHPS 
SLRAT 


VOL 


propulsive coefficient 
same as RGEND in MAIN 
Reynold’s number 
wetted surface, sq ft 
same as SFCHHP in HAIN 
game as SFCMHP in HAIN 
shaft horsepower 

speed length ratio 
same as VEND in MAIN 
speed in knots 


volume of displacement 
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3.6 Subroutine CRVAL 
3.6.1 Introduction 

Although this routine is physically one of the 
longest in the cutter model, the job it performs is one 
of the simplest. Its only function is to choose from 
stored data arrays the values of the residuary resistance 
coefficient, Ce which straddle the input point and 
then to interpolate to find the value of C, at the input 
point. The input point is specified by four variables. 
The value of each of these four-beam to draft ratio, 

Cpr C,, and speed length ratio-is specified as input 
to the routine. 

The constants stored in the three arrays which 
contain the C, data were determined using the curves 
given in Reference (7). Although the three arrays - 
CRi, CR2, and CR3 - are each one dimensional, the data 
they contain is four dimensional, varying with each of 
the four variables which specify the input point. Each 
of the arrays contains values for beam to draft ratios 
of 2.25, 3.0, and 3.75. For each beam to draft ratio 
Cy was varied from 0.48 to 0.70 in steps of 0.02. For 
each of these values of B/T and Cp, all the possible 
combinations of Cy and speed length ratio given below 
were stored, 

In array CR1 is stored data for C, values of 0.001 
and 0.002 and speed length ratios of 0.5 to 2.0 in steps 


of 0.1. Array CR2 corresponds to Cy values of 0.003, 





0.004, and 0.005 plus speed length ratios from 0.5 to 
1.3 in steps of 0.1. Data for a C, value of 0.006 and 
spsed length ratios of 0.5 to 1.0 in steps of 0.1 is 
stored in CR3. 

In each of the three arrays, the first element 
corresponds to the lowest value of each of the four input 
variables. Speed lengcn ratio is increased most rapidly 
followed by Cpe Cs end finally beam to draft ratio. 

The limits of the C,. arrays given above are extended 
somewhat by extrapolation. Beam to draft ratios between 
2.0 and 4.0 are accepted. The limiting value on C, has 
been extended by 0.001 for speed length ratios between 
1.0 end 2.0. The mezimum ranges that are accepted are 


given below. 


2.0 < B/T < 4.0 
48 < ° 
8 Sn < .70 
ao61 <C. < .006 for 0.5 < V//L < 1.3 
NOOINEG) =000gm ferm0.5 <V/JL < 2.0 


Sixteen values of C,, are required to straddle the 
input values of B/T, Cp, Cy, and V//L. after those , 
sixteen numbers have been chosen, the value of Cc. at the 
input point is found by interpolation. A total of fifteen 
interpolations are required, eight to select the correct 


speed length ratio, four to select the correct C two 


p’ 
to select the correct beam to draft ratio, and a final 
one to select the correct C,. 


Even with the extensions listed earlier, the constraints 
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on the combined C, array cause the major limitations 
of the cutter model. The C.. array data is listed in 
Appendix D. A flow chart for this routine is shown in 


Figures 12a and 12b. 


3.6.2 Input List for Subroutine CRVAL 
BIReccccevcscseee enc DA 
CPrccccccesseceeeeQCDA 
CR1(1152) cece ee lC(TT) 
CR2(972)ccccccs ell (TT) 
M216 ics cee oo lrd ) 
CVeccseccccccccveseeDCDA 
Vic ccceccccesseee eC DA 


3.6.3 Statement Descriptions 
I1(J,K,L,H)2(K-1 ) *#384+ ( J-1 ) #1924 (L-1 ) #16+H 
I2(J,K,L,M)=( Kel )#324+ ( J-3) #108+ (L-1 ) *944 
13(J,K,L,M)=(K-1 ) #72+(L-1 ) 64M 
These are statement functions which convert the 

four dimansions - J, K, L, and H - into a single dimension 

for use with arrays CR1, CR2, and CR3 respectively. 
IF(BIR.GE.3.) Ki=2 
EF( BTR. LT. Be ) Kizl 
K2=zKi+i 
An index is chosen in these statements .for the 
beam to draft ratio below the input point and above the 
input point. 
CPiz.68 
3 IF(CP-CP1+.0001) 1,2,2 
1 CPlxCP1-.02 
GOTO 3 


2 LieIFIX((CPi-.46)*50.1) 
L2aL14+1 






a 


a 
a ee 
; 


Jit 
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CRVAL SUBROUTINE 
FLOW CHART 














BTATEAIENT ? S 
IVF ALM) = (K- «384 + 
CJ-1)* 192 + (L-i) eve tm 
205, K, LM) = (K-19 5244 
(J-3) a fo8 + Ci-1)e9 eM 
J36T, K,L,a)= (K-72 4 
Z-1)* 64m 


LABELED 
Common 


STATE - 
MENT 













START CRVAL 






Liz tFIX((CPi-. 4b) % 
So. 


t) 






£2= itl] 
CVI < -005 


No 
cpiecpi- (WE eo) YES ot gre] 
JizJ-! 


K2= K/+l 
GVi= CV(— OO [ 













YES 


CV4.0000! No 
. < Cvi e 








Cviz CVI 
200! Mi= LeIX((WLI-- 4) * 
10./) JZ= Ti+ 
M2: Mi+tl 
INOEX< |] 
cont. 
NEXT 
PAGE 





Figure 12a 
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Aa C7 
J~JI2 


TéMOEx <2 





aJo 


NES 


Dp2c7 
Cee (D-A*(cv- CUI} 
/000.tA 














DT=LICJ), K2, cise) 
CeE® cailf) 

Cail), Kz, Li, M2) 
Carscel 5 
eat Seo ee) 


DeLils,Ki,c4,me) 
CHAZ CRI CL) 
T22i(F,xX1,41, M2) 
C£62CF1(T) 

l=) Coie ena 











cac> Ceict) CRO= CRIT) 
Is Ll, Ky, KUM2) To2i(3,K2,62,"2) 
CKo=cR1(1) CRH=- cer(r) 


T2t20T, es Li, mM 
CRAACRZ(IVN 
TeITeTes, ler zt, 2) 
CREW CRuE ) 
LeLl2ly, Ki 42, Al) 
Cro: CRUZ) 
Bate 2 £24.Mnt 
CRO= CAT) 


Te I2ly,Kr,b1, mM) 
c@ke- CRUTY 
LeTUCJ, Kr, Ul, mY 
C22 C€2422Ct) 

CRar CH 

Lertls «2,L%,me) 
CRHeCRCT 
















if »FULI,M1) T-I303,K2,LIM?) 
CRA CR3(L) CRED CRRCT) 
Bert Cy Ki, L,m2) Le Z3CI RL, Lt M2) 
CRE = CRICI\ CRE: CR3Z(Z) 


Ler (5 KF, cum) 
CRe CA3CL) 
T=ef3l,KHl pLU,N2) 


CLO» CL3CF) 


I-F3 CJ ,K2, L2,m!) 
CRG > CRBC) 

E-l2 LJ KU, L% MAU 
Q@RH=cRICr 











V2 (VE-VLI}E JO. 
Cz (CP-CP1) #50. 

BIR FLOAT ( (KINI) 
ae ( are-ete SA 
Ciz (CRO-CRA) Het CRA 


CU~ (CAO CRO) HV CRO 


C3 (CRE-CREJAV 4 CRE 
C4 -(cRW-cRG)tV+CKO 
CEA C-CI)ACHC! 
CG 2664-63) MHC TCS 
C2 = 6C6-C5)ME 4 CS 







END 


RETIVURDO 


Figure 12b 
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In these statements, an array index is chosen for the 


Cy below and above the input point. The extra 0.0001 
in statement 3 and the 50.1 instead of 50 in statement 
2 are used to account for small errors in the computer 
calculations. 


CVis,. 005 

& JF(CV-CV1+.00001) 5,6,6 

5 CVieCV1-.001 
GOTO 4 

13 JisJlel 
CViseCV1i-.001 
GOTO 14 

6 JisIFIX(CVi#1001. ) 
IF(J1.EQ.5-AND.VL.GT.1.) GOTO 13 
PC J1.BQ.20AND. VireG?.1.3) GOTO 13 

14 J2Z=J1i+1 


Array indices for C, are chosen below and above 
the input point. The statements connected with statenent 
13 extend the limits of the array by extrapolation. 
VL1i=1 e 9 
7 IF(VL-VL1+.0001) 8,9,9 
8 VLisVL1 = 1 

GOTO 7 
9 MIisIPIX((VL1-.4)#10.1) 

M2=hi1+1 

Again. erray indices are chosen, this time for 
Speed length ratio values. 


ec ait 
INDEX=1 


The value of C, is held fixed each time through 
the following statements. When C.. values at each of the 
two C, values have been calculated, an interpolation 
is done to find the C, value at the correct C, value. 
11 IF(J.GE.3) GOTO 10 


IeIi(J,K1,L1,éM1i) 
CRA=CR1i(I) 
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I=I1(J,K2,L2,M2) 
CRHeCR1 (I) 
GOTO 50 


Variables CRA to CRH are used for the eight C,. values 
that remain when C, is held constant. In this set of 
statements array CRi is used. 


10 IF(J.EQ.6) GOTO 12 
IsI2(J,K1i oui 1) 
CRAsCR2(I) 


I=12(J,K2,L2,H2) 
CRH=CR2(1) 
GOTO 50 


This group of statements is identical to the proceeding 
group except that arrey CR2 is used. 


12 I=I3(J,K1,L1,}1) 
CRA=CR3(I1) 


InI3(j,K2,L2,H2) 
CRH2CR3(T) 


Again these statements are identical with array 
CR3 being used. 


50 Vs3(VL-VL1 )*10. 
C=s(CP-CPi )*50. 
BTR1=FLOAT( (Kiel ) )*.754+2.25 
B=(BTR-BTR1)/.75 
Ci=(CRB-CRA)#V+CRA 
C22(CRD-CRC )*V+CRC 
C 3=(CRF-CRE )#V+CRE 
C4=(CRH-CRG)*#V+CRG 
C5=(C2-C1 )#C+C1 

Ry eos eee 

C7=(C6-C5)*B+C 5 
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These are the interpolation statements required 


to first calculate the value of C,, at the input speed 
length ratio, then at the input Cy and finally at the 
input beam to draft ratio. 

IF(INDEX.GE.2) GOTO 51 

A=C7 

J=J2 

INDEX=2 

GOTO 11 

When INDEX = 1 the lower value of C, is used in the 
interpolation above. The next step is to do the interpolation 
at the upper value of C,,. The value of C, at the lower 
C, is saved as variable A. 

51 DeC7 

CR=(DeA)*(CV=-CV1 )*1000.+A 

RETURN 

END 

One final interpolation is performed to calculate 
the desired value of CR and the routine returns control 


to subroutine HPCALC. 


3.6.4 Output for Subroutine CRVAL 
The only output is the calculated value of CR which 
is returned as a subroutine call dummy argument. 


3.6.5 Nomenclature List 
All variables have the same definition as given in 
the MAIN program nomenclature except for the followings 


A Cy. value at lower Cy, value 
B fraction of BTR interval 
BTR input value of beam to draft ratio 


BTR1 lower value of beam to draft ratic 





Ci 
C2 
C3 
C4 
C5 
C6 
C7 
CP1 
CR 
CRA 
CRB 
CRC 
CRD 
CRE 
CRF 
CRG 
CRH 
CV 
CVi 


INDEX 


Ji 
J2 
K1 
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fraction of Cy, interval 


at correct V//L, lower BIR, lower Cy 


at correct V/JL. lower BTR, higher Cy 


at correct V//L, higher BTR, lower C., 


Cy 
Cy. 
Cy 
C. at correct V//I, higher BIR, higher C., 
Cy. 
Cy 
Cr 


at correct V//L and Cy» lower BIR 


at correct V//L and C_, higher BTR 


p 
at correct V//L, Cp and BIR 

lower Cy, value 

Cy. value at input point 

~ at lower V/JL, C,, and BTR 

C. at higher V/J/L, lower Cy 


C. at lower V//L and BTR, higher C 


and BTR 
p 
C, at higher V//L and C» lower BIR 
C. at lower V//L ana C, higher BTR 
C. at higher V//L and BIR, lower C., 
C. at lower V//L, higher BTR and C,, 
C,. at higher V//L, BIR, and C, 
input value of C, 

lower value of C,, 

Cy value at upper C,, value 

array index 

Indicates whether Cy, is at lower or upper value 


used for either J1 or J2 


index for lower value of om 
index for higher value of C, 
index for lower value of BTR 





Li 
L2 


< 


VL1 
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index for higher value of BTR 


index for lower value of Cy, 
index for higher value of C 


p 
index for lower value of V/J/L 
index for higher value of V//L 
fraction of V//L interval 
input speed length ratio 


lower value of speed length ratio 
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3.7 Subroutine EPLANT 


3.7.1 Introduction 

This routine performs two functions. ‘The first is 
to estimate an average plant electrical load for use 
in calculating fuel requirements. The second is to 
size the generators. Generator size is used primarily 
for estimating weights, but it is also a useful output 
which provides information needed for beginning the 
preliminary design phase. 

Three ships were used in developing electrical 
load relations. These were the 378° WHEC, the 210° WHEC, 
and the 150° WPB. The electrical loads for these vessels 
were categorized as shown in Table 1. This table also 
gives the loads in KW that were used. 

The first two categories, electronics load and arn~ 
ament load, are inputs. Propulsion auxiliary load includes 
only equipment needed for the direct operation of the 
main propulsion machinery. For example, lub oil] service 
pumps are included in this category but lub oil transfer 
pumps are included under other auxiliaries. Air conditicn- 
ing and ventilation also includes electric heaters. 

Hotel] is primarily galley and lighting. 

Estimating relationships were developed for each 
of the last four categories in Table 1. If the propulsion 
auxiliary load is not input, a value of i KW is assumed 
for diesel plants and 71 KW for CODOG plants. Air 
conditioning and ventilation is calculated as a function 





Table i 
ELECTRICAL LOAD REQUIREMENTS 








Electrical Load 150° WPB 210° WHEC 378° WHEC 
Electronics 4.8 21.99 117.84 
Armament 42.7 0.0 46.64 
Propulsion Aux. 1.02 71 71.03 
Air Cond. & Vent. Zhe 49.13 308. 34 
Hotel 20.9 3504 170.51 

Steering gear ° 54 4.48 37.38 

Aux. Machinery 18.26 40.08 70-72 

Shops 1.0 22 4,98 
Other Aux. (total ) 19.8 46.76 113.08 
Total 110, 32 153.99 827.44 


All loads are in KW 


of cubic number times the number of accommodations. 
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Hotel 


load is a function of the 
Other auxiliaries load is 
Cubic number only. These 
Figures 14, 15 and 16. 
The total electrical 


Sum of the six categories 


number of accommodations only. 
assumed to be a function of 


relationships are plotted in 


load required is equal to the 
above with a two percent growth 


Margin on all equipment plus a 100 percent margin on 


electronics. An additional 25 percent is included above 


the growth margin to allow for detericration and starting 


loads. 
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NAC x CN x 107° 
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It is assumed that three generators will be installed. 


Each generator is sized so as to be able to carry the full 
electrical load. This represents the current Coast Guard 
practice for new designs. Less capacity has been installed 
on past designs. 

The smallest generator that will be installed is 
assumed to be 100 KW. Other possible plant sizes lie 
between 100 KW and 250 KW in steps of 50 KW, between 
250 KW and 1000 KW in steps of 250 KW, and above 1000 KW 
in steps of 500 KW. The smaliest generator that is able 
to supply the total load with the margins stated above is 
chosen. 


A flow chart for this routine is shown in Figure 17. 


3.7-2 Input List for Subroutine EPLANT 
BIOAD. oe 0 ccceccieeeeeb&(DD) 
Cilivn » 0 cususc cece cece sols (BB) 
ELOAD. eee eeccev veces lc (DD) 
JOPT. ccccccvescceeselK(CC) 
TPE weer cc sence veel (CC) 
NCPO. cevcvcccvcecvevecell (AA) 
NENLeccceveveccveee eLO(AA) 
NOPF.ccccccccsccvceelLC(AA) 
PALOAD. oeccvcecccece cl (DD) 





START EPLANT 









KiU/= ELoADr BLOAD tT 
PAL ORD + CALODAD +HLOAD 
+ AVYLOADS 

AvVGKW= KkW-8ro0ad 
iwi 2 bh UIS*KKW+ELDAD 


ELKW:/00. 








ELkw = 
ELku 500 
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EPLANT SUBROUTINE 
FLOW 


CHART 


REAL AND 
LABELED NO 


COMMON = =]. 


STATEMENTS 





OALORD = .000264S 4CA) 
NAC2NOFEFANCPOtNENL 


AVLOAD2/5.% *S2beN ACKCA) 
Oo00- 


Eckw tZ50 


Figure 17 
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3.7.3 Statement Descriptions 


IF(JOPT.EQ.2) GOTO 1 
The estimate of PALOAD is skipped if a value is 
given as input. 


PALOAD=1. 
IF(MTYPE.EQ.2) PALOAD=71. 


Propulsion auxiliaries electrical load is calculated 
for diesel and CODOG plants respectively. 


1 OALOAD=.0002695*CN 
NAC=NOFF+NC PO+NENL 
HLOAD=.917*NAC 
AVLOAD=15.+.376#NAC*CN/100000. 
KWaELOAD+BLOAD+ PALOAD+OALOAD+HLOAD+AVLOAD 
AVGKW=KW-BLOAD 
KW1i=1.25#(1.02*KW+ELOAD) 


The individual load items are calculated and sunmed 
to give an estimate of the total load. The average load 
is taken as the total load minus the armament load. 
Margins are then added to the total load for sizing the 
generators. 


ELKW=100 ® 

2 IF(ELKW.GE.KW1) GOTO 3 
IF(ELKW.LT.240.) GOT 4 
IF(ELKW.LT.900.) GOTO 5 
ELKWxELKW+500. 
GOTO 2 

& ELKW=ELKW+50. 
GOTO 2 

5 ELKW=ELKW+250. 
GOTO 2 

3 RETURN 
END 


These statements choose the generator rating, ELKW, 
which is the next size larger than the total load with 


Margins 2 
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3.7.4 Output List for Subroutine EPLANT 


MM Mies sisiclsisibiciols e+ ee eo LC (DD) 
EMTs s . ouiile cess cee oe LO DD) 


3-7-5 Nomenclature List 


All variables have the same definition as given 


in the MAIN program nomenclature except for the followings 


AVLOAD 
HLOAD 
KW 

KW1 
NAC 
OALOAD 


SFCH 
SFCM 


air conditioning and ventilation load, KW 
hotel load, KW 

total load without margins, KW 

total load with margins, KW 

number of accommodations 

other auxiliaries load, KW 

same as RGEND in MAIN 

same as SFCHHP in MAIN 

same as SFCMHP in MAIN 

same as VEND in MAIN 
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3.8 Subroutine MACHLQ 


3.8.1 Introduction 

At this point in the program, the shaft horsepower 
and size of the electrical plant have been determined and 
the weight of fuel, lub oil and potable water can be 
calculated. Subroutine MACHL performs these calculations, ' 
In conjunction with the fuel weight calculations, the 
specific fuel consumption, SFC, at endurance power is 
estimated. 

For diesel plants it is assumed that the SFC at 
rated power is 0.42 lbs/shp-hr and at one half rated 
power is 0.46 lbs/shp-hr. A linear relation is used 
for other power ratios. The SFC used for CODOG plants 
is 0.42 if the endurance power is less than 7000 SHP. 
This assumes that the diesels will be sized to be at 
maximum power at endurance speed. If the endurance 
power requirements are greater than 7000 SHP, the gas 
turbines must be used at endurance speed. Three maximum 
power levels for the gas turbines are used - 25000 HP, 
35000 HP and 50000 HP. At each of these horsepower 
levels the SFC is taken to be 0.50 1bs/shp-hr. At half 
power, the SFC is .595. A linear relation is used for 
other power levels. 

If a specific machinery plant is input, the specific 
fuel consumptions at half and full power are given. A 


linear relation is used for other power levels. 
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The specific fuel consumption of the ship service 


generators is assumed to be about 0.48 lbs/shp-hr or 
0.65 lbs/KWehr. The generator fuel consumption is computed 
using the average electrical load. 

Fuel consumption for hotel services is taken to be 
0.32 lbs/man-hr. This assumes some waste heat recovery. 

The lub oil requirements for diesels is assumed to 
be 5 tons plus 1 ton per thousand horsepower. For CODOG 
plants this same relationship is used for the diesel part 
of the plant (endurance power is used in computing the 
requirements). However, additional lub oil for the gas 
turbines is included. For the gas turbines, the oil 
requirements are 3.5 tons plus 1 ton per ten thousand 
horsepower. 

Potable water requirements are taken to be 50 gal/man 
or 0.186 tons/man. A flow chart for this routine is 


shown in Figure 18. 


3.8.2 Input List for Subroutine MACHLQ 


AVGKW. sc cccscceeeeeLC(DD) RE wccccccccccccces cll (CC) 
OPT sss cse seve cews LG(CC) SFCHeccccccccccese ell (CC) 
MPYPE ccc cccc cece LC(CC) SPCH. wcccccccccseeeLC(CC) 
NCPO. occ ccccccceeeLC(AA) SHES ais oe ciciclccies see LO UCe) 
NENLececcccccee ee eLC(AA) SHPH. ccc ccccccccee el (CC) 


MME Te. <4, 6000 0 ¢ 0 os LC ( AW) Mee relueee ne a-clere seb (CO) 
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MACHLQ SUBROUTINE 
FLOW CHART 


REAL AND 


LABELED | NAC = NOFFeNChot 
START MACHLQ aa NENL 
STATEMENTS WTPS= -/8S*® NAC 















WTLOs 


No ,OOlASHEM ES. 


YES 


WTLO = 
001M SHPM 
f.ov1*SHPE 
tE3.s 


SHPE < 7000 











SFC = 
sFome+(i—- SHee\e2# 


(sFcH-SeC™) 


















SFC=.42 






FR MAW « SFCHSHPE 
FRHOTL 2 3 2H#NACT 


FRG EN) = 17.778 -, 000076 & 
BG Kal ™4 ABHAVG Ke 
FR=FRMAIW TFL both + FRGEN 


WTEU EL? FReRE 
idpeeve 


EwWD 








fo] 
SFC = S ~08¥ pred 
ShPH: SHPM 
2 






SHPM> 7000 





b 
HPe< SHeH 






SFC = 
(1.387 38 


She) a 
25000 





SHPM > 25000 





Figure 18 
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3.8.3 Statement Descriptions 


NAC=NOFF+NC PO+NENL 

The total number of accommodations is calculated, 

WTPS=.186*NAC 

Additional items will be added to the weight of 
personnel stores later but for now it includes only the 
weight of potable water. 

IF(JOPT.EQ.2) GOTO 1 

This skips the lub oil and standard SFC estimates 
Bince these values are inputs. 


WILO=.001*SHPM+5. 
IP(MTYPE.EQ.2) WTILO=.0001*SHPi.001*SHPE+8.5 


The lub oil weight for diesel and CODOG plants 
respectively are calculated here. 


IF(MTYPE.EQ.2) GOTO 2 

SFC=. 5-. 08*SHFZ/SHPM 

SHPH=SHPM/2. 

IF (SHPM.GT.7000..AND.SHPE.LT.SHPH)SFC=. 5=-.16*SHPE/SHPM 
GOTO 3 


These statements calculate the SFC for diesel plants. 
In case the shaft horsepower is greater than 7000 it is 
assumed that four diesels will be installed. If the 
endurance power is less than half of the maximum, only 
two diesels need be in operaticn at endurance speed. 
2 IF(SHPE.LT.7000.) GOTO 4 
SFC=(1.38-. 38*SHPE/50000. )*.5 
IF(SHPM.LE.25000. )SFC=2(1.38-. 38*SHPE/25000. )*.5 
IF(SHPM.GT.25000..AND.SHPM. LE. 35000.) SFCa(1. 38- 
¢ 38*SHFE/35000. )*. 5 
GOTO ? 
4 SPC =. Z 
GOTO 3 


For CODOG plants the SFC is calculated using the 
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statements above. 


1 SFC=SFCM+(1.-SHPE/SHPHM) #2. #(SFCH-SFCM) 

If SFCH and SFCM are input, the above linear relation 

is used for calculating the SFC at endurance speed. 

3 FRMAIN=SFC*SHPE 

FRHOTL=. 32*NAC 

FRGENs. 65*AVGKW 

FReFRMAIN+FRHOTL+FRGEN 

All of these fuel rates are in lbs/hr. 

WIFUEL=FR*RE/(1788. #VE) 

RETURN 

END 


The factor 1788. is derived from: 


Fn*Faeg”"se 
where s 


Fi z .05 is a tailpipe allowance 


Z =z 1,03 is an extra margin 

Paeg™ 1.05 is a degradation allowance 

Feo = 1.10 is an allowance for sea conditions 
and fouling 


The value 2240 converts pounds to tons. 


3.8.4 Output List for Subroutine MACHLQ 
WRRUEL 6 06 cece ere ee LC(LL) 
Wee sc cccc cc ec ecole (LL) 
WTPS oo cccvcccccccveelLC (LL) 


3.8.5 Nomenclature List 
All variables have the sams definition as given in 


the MAIN program nomenclature except for the followings 
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Nomenclature List (cont) 


FR 
FRGEN 
FRHOTL 
FRMAIN 
NAC 

RE 

SFC 
SFCH 
SFCM 
SHPH 


total fuel rate in lbs/hr 

generator fuel rate in lbs/hr 

hotel fuel rate in 1lbs/hr 

main engine fuel rate in 1bs/hr 

total number of accommodations 

Bame as RGEND in MAIN 

specific fuel consumption at endurance power 
same as SFCHHP in MAIN 

same as SFCMHP in MAIN 

one half maximum shaft horsepower 


same as VEND in BAIN 
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3.9 Subroutine PBSIZE 


3.9.1 Introduction 

Very little data is available concerning the minimum 
allowable depth and length of the machinery compartment. 
As a result, only very simple calculations are made in 
this subroutine to estimate these two dimensions. 

The product of length of machinery box times the 
maximum beam of the ship was calculated for two diesel 
designs, the 210° WMEC and the 150° WPB, and for one 
CODOG plant, the 378° WHEC. A value of about 1050 sq. 
ft. was obtained for both diesel designs and a value of 
2300 sq. ft. was obtained using the CODOG plant. These 
values are divided by the beam of the ship to obtain 
the required engine room length in the progran. 

A check is also made to insure that a minimum length 
is not exceeded. This minimum length was taken to be 
28 feet or 56 feet for diesel plants below or above 
7000 BHP respectively and 45 feet for a CODOG plant. 

A rough deck height approximation was derived using 
data from the same three vessels mentioned above. The 
formula used gives a deck height of about 17 feet for 
the 150° WPB and about 19 feet for the 210° WMEC. The 
deck height for all plants with horsepowers greater than 
7000 is taken as 22 feet, 

A flow chart for this subroutine is shown in Figure 
19. 
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MBSIZE SUBROUTINE 
FLOW CHART 


REAL AND 


BELED Y 
START MBS!IZE eon 
STATEMENTS 


XLM= 
SHPMS> 7000. 
XLM=56. be SES Manin 7 1050/8, {LM = 2300./8 
XLMA<SG 


XLM=23 













DHM ~Il+ .001S7x 


SHPM™ 


ES 
DHM=22. eon) 
No 





Figure 19 
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3.9.2 Input List for MBSIZE Subroutine 


Pius 6 sugges 0 00 066006 oD UDA 
WATE PEisisis «oss cee oLl(CC) 
UREA reve: ¢.0,08.01010 @BOKCEC ) 


3.9.3 Statement Descriptions 

IF(MTYPE.EQ.2) GOTO 1 

XLM=1050./B 

IF(SHFU.GT.7000.. AND. XLM.LT. 56. ) XLM=x 56. 

IP(SHPM.LE.7000. AND.XLM.LT.28.) XLi=x28, 

GOTO 2 

These statements set the engine room length for 
diesel plants. (MTYPE =» 1 for diesels) 


1 XLM=2300./B 
IF(XLM.LT.45.) XLM=45. 


These statements set the engine room length for 
CODOG plants. 
2 DHM=11.+.00157*SHPM 
IF(DHM.GT.22.) DHM=22, 
RETURN 
END 
The minimum deck height, which is the same for both 


Plant types, is set by these statements. 


3-9.4% Output List for Subroutine MBSIZE 
All variables have the same definition as given in 
the MAIN program nomenclature except for the following: 


RE same as RGEND in MAIN 
SFCH Same as SFCHHP in MAIN 
SFCM same as SFCMHP in MAIN 


VE same as VEND in MAIN 
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3.10 Subroutine VOLUME 


3.10.1 Introduction 

Many of the ships designed today are volume rather 
than weight limited. Because of this it is important 
that a ship synthesis model determine a balance between 
required and available volume as well as between weight 
and displacement. Subroutine VOLUME performs this function 
in the cutter model. 

The dimensions of the hull below the load waterline 
are determined in subroutine UWDIH. Only the sheer line, 
deckhouse volume and the length of raised deck, if any, 
remain to be determined to define the volume available 
in the ship. The enclosed volume of the hull and deck- 
house will be completely determined by these variables 
since the waterplane coefficient is calculated in sub- 
routine UWDIM and since the general shape of the hull 
is implied by the estimating relationships used in the 
model. 

Deckhouse volume is constrained primarily by stability 
considerations and by external deck area requirements. 

If the deckhouse becomes too high, it will not be possible 

to make the vessel stable. Also a certain amount of 

free deck area is required in nearly all designs, restricting 
the horizontal spreading of the deckhouse. At the other 
extreme, the ship may have excessive stability if too 

Small a deckhouse is installed. 


These deckhouse volume restrictions have been included 
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in the cutter model by restricting deckhouse size to be 


within the limits of past practice. Both an upper and 
a lower limit are placed on the size of the deckhouse 
as shown in Figure 20. The deckhouse volume of Coast 
Guard vessels has been generally greater than that of 
Navy designs. Upper and lower limits of 0.00150 x L? and 
0.001 x L?, respectively, are used in the Navy‘’s destroyer 
model. This difference is probably due to the fact that 
the Navy requires more external deck area for armament. 
The deckhouse volumes shown in Figure 20 do not 
include uptekes. It is assumed that sufficient additional 
area is available for uptakes since uptake area was 
excluded from the data points also. However, it is 
{jmportant to check a design be sketching a layout for the 
external decks and the superstructure before acceptance. 
The volume in the hull is also eatimated in VOLUHE. 
The underwater volume is already available. fhe first 
step in determining the above water hull volume is the 
development of an acceptable sheer line. This step is 
performed by subroutine SHEER and is discussed in con- 
junction with that routine. 
Once the sheer line and with 1+ the average freeboard 
have been estimated, the above water volune is found 
by multiplying the area of the waterplane by the average 
freeboard and by a factor to account for flare. The 


total volume in the hull is the sum of the above and below 


water volumes, 
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This total hull volume includes the volume to the 


main deck only. The hull volume cannot be reduced below 
this size but a raised deck can be added to increase the 
volume of the hull. For vessels of the size of Coast 
Guard cutters it is impractical to add more than one 
raised deck. Also the deckhouse size can be varied some- 
what to change the available volume. - 

Required hull volume is broken devwn into three cat- 
egories. These are the machinery box, tankage volume 
and arrangements volume. The object of subroutine VOLUME 
is to alter the available volume so that there is sufficient 
space for each of the three categories of required volume. 
The model essumes that all tankage volume and the machinery 
box volume are in the hull. The deckhouse contains only 
arrangements volume. 

The machinery box size is the first category estimated. 
This volume is found by multiplying the midships area 
by the length of the machinery box and by a machinery 
box prismatic coefficient. Machinery box volume is 
subtracted from the available hull volume to the main 
deck leaving a volume which can be used for tankage and 
arrangements, 

some of the remaining volume cannot be used for 
errangements because of hull shape. The next step is to 
estimate the fraction cf the remaining volume that can 
be used only for tankage. This is the minimum tankage 
volume for the ship. If less tankage is required for 
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the liquids that must be carried, the excess volume must 


pe used for voids or cofferdams since it cannot be con- 
verted into machinery box space or arrangements volume. 
If more tankage is required than this minimum, some of 
the arrangements volume must be used for tankage. 

The weights of ail the liquids which must be carried 
are calculated in subrontine HMACHLQ. In VOLUME, these 
weights are converted to volumes and added together. 

A twenty percent allowance is added to this figure to 
account for required cofferdams and voids. The volume 
required for liquids is then compared to the minimum 
tankage volume to determine the amount of arrangements 
volume, if any, that is required for tankage. 

All of the volume remaining in the hull after tankage 
volume and machinery box volume have been subtracted,can be 
used for arrangements. The entire superstructure volume 
is also available for arrangements. Arrangements space 
requirements are more commonly given as areas rather than 
as volumes, so the next step is to convert the available 
arrangements volume to an available arrangements area. 
For the deckhouse this is a simple task since the deck- 
house is nearly wall sided. The volume is divided by 
the deckheight of 8.5 feet to find the area. For the 
hull no such simple procedure is possible because of 
hull curvature. 

The Navy performed extensive work on this problem 
while developing their destroyer model, DDO7. ‘This work 
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was made use of in the cutter model not only for the 


arrangements area formulation but also for the flare 
factor, machinery box prismatic coefficient and tankage 
volume fraction estimates discussed earlier. The relation- 
ships used in the cutter model were taken directly fron 

the original notes (8) of Mr. Jim Hills who was largely 
responsible for the development work done by the Navy. 

The major deviation between the cutter model and the 
destroyer model lies in the engine room location. The 
destroyer model assumes that the engine room is centered 
at amidships but Coast Guard practice usually places the 
engine room aft of amidships. No correction is made for 
this difference and the output produced for test runs 
with existing ships indicates that there is little effect. 

By using the Navy's formula, the arrangements area 
in the hull can be calculated. This area plus the area 
of the deckhouse must be compared to the area required 
for all the ship’s functions plus the area specified as 
input. The area required for ship's functions such as 
stores area, living area and auxiliary machinery area, 
are estimated in the cutter model based on past practice. 

Data was obtained by measuring the areas of five 
past designs and assigning the areas to various categories. 
This data is summarized in Table 2. Linear relationships 
were developed from this data as shown in Figures 21 
through 25. Berthing and sanitary areas are based on 


Minimum Coast Guard standards. fhe area required for 
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PaTeGORY 378 Bey 255 210 150 
WHEC WHEC WHEC WIIEC WPB 


Office Spaces 493 465 347 130 65 
essing Facilities 2713 2348 1379 1452 339 
Crew Special 25 609 7A? 324 32 
Officer staterooms mies? 1338 966 567 187 
Officer Sanitary 243 133 120 193 30 
CPO Staterooms 1040 364 336 283 110 
CPO Sanitary 250 100 - 75 27 
Crew Berthing 3641 2655 2568 1658 375 
Crew Sanitary Zio Lily 492 384 90 
CO Stateroom & Cabin 697 Sie 273 406 - 
CO Pantry 58 128 22 = z 
Commissary Stores 600 72! S31 480 35 
Other stores 3015 2354 2397 1294 435 
Workshops 1161 818 Sf 544 - 
Passages 2567 2588 2138 739 229 
Repair Lockers 148 18 97 110 = 
Steering Gear & Windlass 736 548 645 270 294 
Chain Locker 47 30 36 42 16 
A/c & Fan Spaces 369 324 39 165 - 
fc. «= Gyro room Biz 143 232 SO = 
Aux. Machinery Spaces 1860 1874 4414 820 543 
Uptakes (Hull) 462 1216 64 ~ - 
Pilotnouse, Chartroom & 
Gelewc. cs 856 eeu 422 300 
Sud total 24135 20588 15006 10535 3107 
me TRAS 
Crew Recreation Rooms 1445 - 682 = 71 
Cargo & Flume Tanks 1219 249 - - 241 
Electronics & Stores 1618 VCS 1206 319 157 
Arinament & Stores 844 822 7 30 280 ~ 
| Balloon Snelters 512 156 120 - = 
Massion Cffices & Labs 746 336 482 - = 
Mission Workshops = =~ 70 15 = 
Bow Thruster 344 ~ - = = 
rotal 30863 22926 18376 11242 3576 


Table 2 
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repair lockers is taken as a constant 100 square feet. 


The area for chain lockers as 40 square feet and the area 
for uptakes in the hull as either 65 or 460 square feet 
depending on whether diesels or gas turbines, respectively, 
ere installed. A constant value of nine percent of the 
total area is assumed for passages. Passages ranged 

from 6,84 to 13.17 percent of the total area on the 

five designs used, 

All of the required areas are summed to determine 
the total required area. This area is then compared to 
the available arrangements area in the hull and super- 
structure. 

Comparisons are made in the following order. First, 
the arrangements area required to be located in the 
deckhouse i8 compared to the area of the largest allowable 
deckhouse. If the required area exceeds that available, 
the design is infeasible and an error message must be 
printed. Next, the total required arrangements area 
is compared to the total available area. If the required 
area is the smaller, the deckhouse size is decreased. 
The deckhouse size used is either the size which balances 
the area requirements or the smallest allowable deckhouse 
size or the area required to be in the deckhouse, whichever 
is greatest. Should the total required area exceed 
that available a raised deck is added. 

The length of raised deck is determined by an iterative 


procedure. On each iteration the area required to be 
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included in the raised deck is computed and the curve 


shown in Figure 26 is used to convert the area toa 


raised deck length. This curve is based on the main 


deck of a 378' WHEC. The iteration continues until two 


successive iterations fall within 5 percent of the length 


of the ship or until twenty iterations is reached. In 


the later case the design is ruled infeasible. A check 


js also made to insure that the raised deck length does- 


not exceed the length of the ship. 


If the raised deck length lies between 0.4 and 0.6 


times the length of the ship, the raised deck is extended 


to 0.6 times the length of the ship. 


A flow chart for this routine is shown in Figures 


27a and 270-6 


3.10.2 Input List for Subroutine VOLUME 


MEP MTHs .. oo elise os LOEB) 
RMIDATIL,. osc cs cee cele (Be) 
MMs soc ce ce cssees so LONGG) 
SIaMe\ 6 01610) 01 6.6 elsheners’s LONIGG ) 
GMP. crcsccvcccces sLC(GG) 
MEM 6.6.axe's.s.0 ese HOGG) 
RM ss 60:0 a's boos LOGAN) 


DMBecccsvccssccceelt (JJ) - 


BMANs6 6.0 6 0s 06 oe stleniGG) 
TORT irre seeeeeeeceLC(CC) 


WBN. aisers's so 0eekes ee LC (BBY 
MEM BM ioc: sieeve’ e sieleverererere LC (J J} 
PRY. + ese viene LO(CC) 
NOPOs se ecccsesce oo eLC( Aad 
SEM Iie -aiavsie:c sleraie ave ote LC (RAD 
NOMMieis cles vieevelere cea LO WAR) 
WAGFUL « « cicies see's oe LC (LL) 
WICUEL sess eccc sees L0 (LL) 
WILO. covcccsccveecslt (LL) 
WErSe.sceccce cece esl (LL) 





t 
Ted 
oe 
° 
. 
>. 
. 
1 
4 
ee 
















































e > 
~miwe ‘ aa | Fee < 
ee on Es 
aie ee - a 
-_ #@eee aaa? 1 
\- . ; see eh pe te te he fe eee 
. e+e . b t- — emf em fee fee eo 
sete $ -4 eo - dipit-e f- + mee ee ay 
° a jee = —_— s+ —- - -_—~--Ft . 
‘ Go = a ° enone = 4% ep ee ee ee oR eee + ee - —_— ed aa ho ar ~—- oo a 
ie ee ee ee | mee el et stot : | 
-}°- sends 7 aie eax one —_— ee -+- i —- +2 -—p-« 
e eas = 6 © § « } ~ $= © owe +t Se 0500 « ~ eS eee cane 
_— - i a oe -¢+ 
es ee tna > o fo eaten eons dh whe: 
oe -a* & cee ee et - 
oeee . ee t --<o. 1. oe -—- 
ot be tebe : + Op baw Ot Oap ee e 
eb es o- —- -—— —., wae ee wef 
‘- - e 4 o- mah ee st eee i. —- - 
ee -e eg ° 3 eee b= ee hd ee ———— —— ++ oo — | -- e- 1 26 fe 
ete eet «<e*s2 @ eeomehteune: +. me ¢ ewe te wee fae er ——- 3 —e eee —e f- 2 — > 
- oe ; ~s s+ fe —-—4—: -= oan. (eee —_—_ ort meted Meet lal wat 8 Pee ee 
—— oe ae oe > = ow —-_ 
. . '"@e e ’ om 8 coon eg aa we & eee ote “c-- 
a | o- 6 a > : i ——- + ——+ + Hee oy © 
| ‘. t ¢ . ——+ ee 5 wens tae =e one. oe mtee Gore 
ag . . ° 3 » : - ~-~e = ee ne a — ban bate eee 
- ° -_ +e & Gaon th een eee ep eee 
ee. = 5 ’- a are > (a be meee fe rr be eee tw pw meee pe ob, 2 
mse - a : a a ee Rae | ~-7 * -we- ¢ - ome “ga: eds eee 
i cE . ' 3 - eatin esos Sel ate on 8 enlaa h 6 en © = ee 
x = , e te ar tt oe ae es a ee ee ene 
- ee — te 6. 6 ee ———- pe 
eee re - =) ene = ee be et ot = ¢t eer peer 
oo §-: oa. ° stm wg maw ape eee 
1 - e . o- eo - os o-% « e-a- a 
a <6. - es ee oe + - 
_ ee oe — —— _ —_ Sad 
« @ ee - -* 
ah . SSA =< 
ee . - - . -~et- 


7 
+4 


‘ 


t 
2 ww tne a GP ee te ee 
e 
- 
' 











ooe 
ee 
° 


ja 
‘ 
t 
aS 
-| 
1° 


. 
. 
=e 





° 
-_— > 









ee Se 
te - - 
Si: 2 : 
> ' = 5 
. - - 
° . e ¢. 
. eee 
es, lene 
ot > 
oe > 
- ie > 
_ epee . + - 
- ~—pPp-zece 6 
e 
® 
° 
ae 
ae 
be 
e 
ee 
at 
re 
fe 
fe oT c 
a 5; rt Sey + e 
me 1 © ae amek 
: oe 
"4 
be-e 
= ° © 
He eee ase eg poe 
be ee y owe hoe- = 
ee . - 
bes 





MO3SG GsasiVY HLONSA 


12 


Figure 26 


RAISED DECK AREA 
BEAM x LENGTH 


X= 





126 


VOLUME SUBROUTINE 
FLOW CHART 
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Figure 27a 
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3.1003 Statement Descriptions 


KNDEX=0 

NACONOFF+NC PO+NENL 

DECKHT=0. 

LRD=0. 

LRD1=0. 

KNDEX is an ;teration counter used for the raised 
deck length iteration. The variable DECKHT is explained 
in the SHEER subroutine. The two values of raised deck 
length start at zero. 

13 CALL SHEER (DMB,FO,F10,P20,D0,D10,D20,LEN,H,DECKH?) 

The sheer line of the vessel. is defined by the 
freeboard at stations 0, 10, and 20. These values are 
returned as FO, F10, and F20 respectively. The associated 
hull depths are aiso returned. 

AREA=LEN/6. © (PO+F20+4*F10) 

FAVG#AREA/LEN 

DAVGeFAVG+H 

The average depth i8 calculated using the projected 
lateral area. A better estimation could be made by using 
the local beam to weight the freeboard values but this 
was not done,.because the value of F7 which corrects for 
flare is based on the above average depth. 

p7a~ 4018828 (DAVG/H) ##2+ 6 18098*DAVG/HH « 71599 

This is the flare correction factor. The formula 
used was taken from the Navy's work in developing their 
destroyer model. 

HVAWaLEN*B#CWP*PAVG*F 7 


HV BWaLEN*B*H#C P#CX 
THV=sHVAW+ HV BW 


The total hull volume +o the main deck {g calculated 
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by these statements. 


AbiaB*H*CX+B*F10 

Xe LMB/LEN 

CPMa-, 525*X*X+.055*X+.99675 
CPN=C Pit+ 3 ’ 2eX#X*(CP-.6 
VMB=LiiB*AM#C Pil 

AHV OM=THV -ViB 


The volume of the machinery box and the remaining 
available hull volume are calculated in these statements. 
The formula for CPM has been taken from the Navy's work 
and applies to machinery boxes centered about amidships. 


fri 90371154. © 2139727 #D1 0-1. 38263E-2*D1 OF #24 

008058E-4*D10##3-5, 489481E- 6*#D1 OF #44 
2.892153E-8"D10##5 )#CP/, 5B*CX/, 814% 
(1.11 *LMB/LEN+.667) 

ATVOM=Q*AHVOM 

VFUEL=WTFUEL*43. 

VLO=WTLO#39. 

VFW=.1862NAC #36. 

VACFUL=WACFUL#44, 1 

VRE Q= ( VFUEL+ VLO+ VF W+VACFUL) #1.2 

RTAA=0,. 

IF(VREQ.GT.ATVOM) RTAA=(VREQ-ATVOM)/8.5 

AAV CMs AHY OM-ATVOH 

TANKVL=ATVOM+RTAA*8. 5 


The available volume which can be used only for 
tankage and the required tankage volume are calculated 
and compared in these statements. The formula for Q 
was again taken from the destroyer model work. Ths 
“required tankage arrangements area is calculated if the 
required tankage volume is in excess of that available. 

Px (6.51041 E-7#D1 0743-7, 85148 ~5*D1 0%#24+3.438375E-3% 

D1 0+.0530399 )#( 191 *CP+, 6893)7(. ea Loe sd 
(=. 105*LB/LEH+1 0315) . 
AVAA@AAVOM*P-RTAA 
The available arrangements volume is converted to 


an available arrangements area again using a Navy developed 


formula. 
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DHVs, 00223*LEN*#*#3 
DHAA=DHV/8. 5 
CN=( DAVG+LRD/LEN*8, 5) *LEN#B 
RDHA= ( AREADH+24.4.00225*CN+2037#LEN )#1.09 
IF(RDHA.GT.DHAA) GOTO 1000 
TAAA=DHAA+AVAA 


The area in the largest acceptable deckhouse is 
compared to the area required to be in the deckhouse.' 
If the required area is smaller, the program continues. 
If not, an error messazce is printed. The areas required 
to be in the deckhouse include the input areas plus the 
commanding officers cabin and the piiothouse. An allowance 
is also made for passageways. 


ODA=2. 12*LEN-232. 


PHDA=100.4+.00225*CN 
The statements above and those included between 
them are used to calculate the required deck areas. 
TRAA=1. 09*( AREADH+AREAHL+ ODA+XiF DA+CSDA+ OSRDA+OHDA+ 
CPODA+C POHDA+C BDA+C HDA+C ODA+C5S TDA+ CS DA+WDA+ 
RLDA+C LDA+UTDA+S GDA+ACDA+XICDA+AHDA+PHDA ) 
IF(TRAA.GT.TAAA) GOTO 10 
The total required arrangements area is computed 
and compared with the area available in the hull to the 
Main deck plus the largest accaptable deckhouse. A 
decision is then made to either make the deckhouse smaller 
or to add a raised deck. 
DHS V2. 00151 *#LEN#*#3 
DHAMI N=ANAX1 (DHS AA,RDHA) 
TAAA=DHAMIN+AVAA 
IF(TRAA.GT.TAAA) GOTO 11 


12 DHV=DHSV 
DHAA=DHV/8. 5 
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ENCVOL=DHV+THY 
RETURN 
11 DHAR=TRAA~AVAA 
DHS V=DHAR*8 e 5 
GOTO 1l2 


These statements compute the deckhouse size required 
if the decknouse is to be made smaller. First, the 
larger value of the minimum deckhouse size or the required 
deckhouse area is chosen and added to the hull voulume. 
This new available area is then compared to the total 
required area. If the available area is still the greater, 
the difference is excess area that cannot be reduced 
since the minimum size hull and deckhouse are already 
being used. If more than the minimum deckhouse size is 
required a calculation is made to balance the total 
required area and the total available area. 

10 RHAA=TRAA-DHAA 

RDAA=RHAA~AVAA 

X1=RDAA/(B*LEN ) 

IF(X1.GT..7) GOTO 1 

IF(X1.GT..2) GOTO 2 

LRD=-3.22*X1#X14+2.031¥*X1-. 0184 

GOTO 3 

1 LRD=.463*X1#X1+.446#X1+.22 

GOTO 3 

2 LRD#X1+.059 
3 LRD=LRD*LEN 

DECKHT#8, 5 

IF(ABS(LRD-LRD1).LT..05*LEN) GOTO 15 

IF(KNDEX.GT.20) GOTO 1000 

KNDEX#KNDEX+1 

LRDiI=LRD 

GOTO 13 

If more area is required than is available in the 
largest deckhouse plus the hull to the main deck, a raised 
deck is added. The assumption is made that a ship with 


the maximum deckhouse size will result in the least 
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displacement. This is probably valid so long as an 


aluminum deckhouse is used. 

The length of raised deck required is computed in 
the above statements. The new raised deck length is then 
compared to the former estimate, LRD1. If the two do 
not agree within five percent of the length of the ship 
anew iteration is made- A default counter has been 
included to keep the computer from being caught ina 
nonconverging loop. 

15 Cis oS*LEN 

C2= a 6*LEN 

pe( ORD. LT.C1.OR.LRD.GT.C2) GOTO 4 

AVAAz=AVAA+(C2=-LRD)*B 

RDAA&RDAA+(C2-LRD)*B 

LRD=C2 

These statements check to see if the raised deck 
length falls within 0.4 to 0.6 times the ship length. 
If so the calculated values are corrected so that the 
raised deck length is equal to 0.6 times the ship length. 

4 IF(LRD.GT.LEN) GOTO 1000 

THV=THV+RDAA#8 e 5 

ENC VOL#2DHV+THV 

RETURN 

A check is first made to insure that a raised deck 

length longer than the length of the ship is not required. 

A corrected total volume is then calculated and the program 
returns. 
1000 LEN=-1000. 

RETURN 

END 

statement 1000 is a default option warning subroutine 


XECUTE to print an error message. 





133 


3.10.% Output List for Subroutine VOLUME 


3.10. 


WMD. + 0 oe slocieierese «LPF ) 
WOON. 0 bc eeetes 00 SOU PP) 
BUA, « Gite + < « vseCGPP) 
CEMA, ccccsite sess Lourr) 
DMs «4s oelmisln neko urP) 
CUD iAens «os cclees oe eUCerP) 
iiss «0 0 so 6 etehble ofS (DE) 
GODAY . « ccleseecsee eLOUPRP) 
ERODAs «ce cceesees o LPP) 
OROHDE.ccsceeesv evel OPE) 
BoB Miccsccceeces o MC( FE) 
SODA es « ad iieevtlels'e LO VRP) 
ID 6's: + ave Gace 6 oo CHOU ED ) 
BOO. 0200s settee ss ee (11) 
IANO, cca lsxsiermettreres ¢ « «OOK La) 
IN Ge. oo oo cleetete ae’e BOND) 
WHA Mrs. 6 « cuteness LOUPP) 


5 Nomenclature List 


Dis eleielel ele eleleie cielo GCs ) 
ENCVOL. o+eeeeee+ eLC(BB) 
WED ines sins sees LC CSG) 
GERM ccs c cece so lLC(PP) 
OHDA Mis ss +e ee eevee LC (PP) 
GSDA MN. «oss sees 0LC(PP) 
OSRDA. ee eeeceeesLC(PP) 
BHOUM. Ge cose sess sLC(PP) 
RGPDAOL.s s oe sass se LCC PP) 
SGUAM + ocoeeeeeslCCPP) 
MAK Ves cles os eee lc (Pr) 
UTDA.scveecccces LC(PP) 
MEN . 60040 se sesebOCEr) 
MDA%is soc ee0c0vessLC(PP) 
FANE DN cle servic o ele GRE) 
XMFDAsseccccceee oLC(PP) 


The definition of all variables is the same as 


given in the MAIN program nomenclature except for the 


follo 
AAVOM 
ACDA 
AHVOM 
AM 
AMDA 
AREA 


wing: 


available arrangements volume in hull, cu ft 


A/C & ventilation deck area, sq ft 


available hull volume, cu ft 


cross sectional area amidships, sq ft 


aux. machinery spaces deck area, sq ft 


projected lateral area above water, sq ft 








aa) 
oe 
i 






ce | 
meeeiiaeeh ti 


Peete ee ates: ,e4% 


peeei iii iii! 


ATVOM 
AVAA 
C1 

C2 
CBDA 
CHDA 
CLDA 
CODA 
CPh 
CPODA 
CPOHDA 
CSDA 
CSTDA 


DO 

Di 0 
D20 
DECKHT 
DHAA 
DHAMIN 
DHAR 
DHS AA 
DHSV 


FO 
F? 
F10 
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available tankage volume, cu ft 


available arrangements area in hull, sq ft 
0.4 x length, ft 

0.6 x length, ft 

crews berthing deck area, sq ft 

crews heads deck area, sq ft 

chain locker deck area, sq ft 

C.O. S.R., Cabin & pantry deck area, sq ft 
machinery box prismatic coefficient 

CPO S.R. deck area, sq ft 

CPO heads deck area, sq ft 

crew special deck area, sq ft 

commissary stores deck area, sq ft 

same as ENDDAY in MAIN 

depth at F.P., ft 

depth at amidships, ft 

depth at A.P., ft 

height of raised deck, ft 

deckhouse arrangements area, sq ft 
minimum deckhouse area, sq ft 

deckhouse area required, sq ft 

emallest deckhouse arrangements area, sq ft 
Smallest deckhouse volume, cu ft 

same as DHA in MAIN 

freeboard at F.P., ft 

flare factor 


freeboard at amidships, ft 





F20 

FAVG 
HVAW 
HV BW 
KNDEX 


LRD1 
NAC 
ODA 
OHDA 
OSDA 
OSRDA 


PHDA 


RDAA 
RDHA 


RHAA 
- RLDA 
RTAA 
SFCH 
SFCH 
SGDA 
TAAA 
TANKVL 
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freeboard at A.P.,ft 


average freeboard, ft 

hull volume above waterline, cu ft 

hull volume below waterline, cu ft 
counter 

same as XLM in MAIN 

used to store LRD value 

number of accommodations. 

office deck area, sq ft 

officer heads deck area, sq ft 

other stores deck area, sq ft 

officer staterooms deck area, sq ft 
coefficient 

pilothouse, chartroom & CIC deck area, sq ft 
coefficient 

raised deck arrangements area, sq ft 
area required to be in deckhouse, sq ft 
same as RGEND in MAIN 

required hull arrangements area, sq ft 
repair locker deck area, sq ft 

required tankage arrangements area, sq ft 
same as SFCHHP in MAIN 

seme as SFCMHP in MAIN 

steering gear & windlass deck area, sq ft 
total available arrangements area, 8q ft 


tankage volume, cu ft 





THV 
TRAA 
UTDA 
VACFUL 
VE 
VFUEL 
VFW 
VLO 


VREQ 
WDA 


Xi 
XICDA 
XMFDA 
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total hull volume, cu ft 


total required arrangements area, sq ft 
uptake deck area in hull, sq ft 
volume of aircraft fuel, cu ft 

same as VEND in MAIN 

volume of ships fuel, cu ft 

volume potable water, cu ft 

volume of lub oil, cu ft 

volumes of machinery box, cu ft 
required tankage volume, cu ft 
workshops deck area, sq ft 

LMB/LEN 

RDAA/(B*LEN ) 

I.-C. & gyro room deck area, sq ft 
messing facilities deck area, sq ft 
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3.11 Subdroutine SHEER 


3.11.1 Introduction 

Subroutine SHEER selects an appropriate sheer line 
by specifying the freeboard at the forward perpendicular, 
amidships, and at the after perpendicular. This sheer 
line must satisfy a number of criteria. 

First, there must be sufficient depth amidships 
to accommodate the engine room. The depth amidships 
must also be sufficient to insure adequate structural 
strength. A minimum value of L/16 is used. Because 
Coast Guard vessels normally have engine rooms aft of 
amidships, a slightly greater depth amidships is used 
rather than the machinery box depth. This accounts for 
the drop in the sheer line aft of amidships. 

The next criterion which must be satisfied is adequate 
freeboard at the forward and after perpendiculars. 

Navy derived formulas were used. These formulas give 
the minimum acceptable freeboard based on a deck wetness 
criterion. 

Once independent determinations of freeboard at the 
forward and after perpendiculars and at amidships are 
made, a check is required to insure that a reasonable 
Sheer line can be drawn between the three. Data was 
obtained from past designs and from the standard merchant 
Ship sheer curves. This data was monaimenal cnaiived by 
Subtracting the freeboard amidships from the other free- 
board values and then dividing by the length of the ship. 
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The resulting curves are plotted in Figure 28. 


The sheer line is assumed to be practical if the 
forward sheer fraction lies between 0.01 and 0,03 and the 
after sheer fraction lies between 0.001 and 0.0075. 

If this is not the case, one or more of the freeboard 
values is adjusted until the criterion is satisfied. 

If a raised deck is added to the hull, a new sheer 
line is calculated. It is assumed that the required 
freeboard of the main deck at the forward psrpendicular 
can be one deck height or 8.5 feet below that required 
by the deck wetness criterion. This will result in a 
flattened sheer line as is common for raised deck vessels. 


A flow chart for this routine is shown in Figure 29. 


3.11.2 Input List for Subroutine SHEER 
TUS ccc meee te eee sn liste che sre eletete@ siete ec LA 
DEN 666s ce eee cee OCA BECHHT cc ccccsveeDC DR 


3.11.3 Statement Descriptions 


Di 0=.0025*LEN+DMB 
IF(D10.LT.LEN/1i6) DiO=LEN/16. 


The minimum depth amidships is selected by these 
statements. 

FOe1.011827#(100. *H/LEN )-.000636215*LEN+2. 780649 

FO=LEN*FO/100. 

FO2F0-DECKHT 

DO=F0+H 

The minimum depth at the forward perpendicular is 


selected here. 
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Figure 28 
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SHEER SUBROUTINE 
FLOW CHART 













INPUTS 


2-78 
LEN 
c H REAL. DIO = .0oze* LEN 
START SHEER ene ae 7 Ee 








FO=1 O18 22XKUIOML EY) — 
s 00OL3EUSHLEV 42, 780649 
FO> LENXFO// SO. 


F20¢ ole LEN 
C2125 +,002S#LEy)) 





Mo 
- FB-DECKHT 
D20 = F2Zo tH c oo; rH 
Yes 
OHIO. LET. 

$$ = DS-O/g 

Lo-= S$/LEN 

$Z¢ = D2g-0/¢ 







R2g- S2¢/LEN 





ZEZL R2d <,001 
oO 
? 
END R2O > ,007S 





YES 


Figure 29 
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F202. O1*LEN#*(2.125+. 00125*LEN ) 
D202aF20+H 


Finally the minimum depth at the after perpendicular 
is selected. | 
7 $0=D0-D10 
RO«S0/LEN 
IF(RO.LT..01) GOTO 1 
IF(RO.GT..03) GOTO 2 
The sheer fraction forward is calculated and checked 
against the acceptable limits. 
6 $20=D20-D10 
R20=S20/LEN 
IF(R20.LT..001) GOTO 3 
IF(R20.GT..0075) GOTO 4 
GOTO 5 
The sheer fraction at the after perpendicular is 
calculated and checked against acceptable limits. 
1 SO=>.01*LEN 


DO#S 0+D10 
GOTO 6 


If the freeboard at the bow is too low the depth 
at the bow is increased. 
2 S0=,.03*LEN 
D10=D0-S0 
GOTO 6 
If the freeboard at the bow is excessive, the depth 
amidships is increased. 
3 S202. 001 *LEN 
D202D10+S20 
GOTO 5 
If the freeboard aft is too low the depth aft is 
increased, 
4 S$20=.0075*LEN 


D10=D20-S20 
GOTO 7 
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If the freeboard aft is excessive the depth amidships 


$s increased and the sheer fraction at the bow is rechecked. 


5 FO=D0-H 


F20=D20-H 
Fi0=D1 0-H 
RETURN 
END 


Values of both freeboard and depth are returned 


for use in subroutine volume. 


3.11.4 Output List for Subroutine SHEER 


Milmtsictews’e @is)s e eles s «oC DA i) Oita etele cieieiele es eee s OULU 
Putlteses ee 6s ¢ susele eOC LA DmOeae eeewies se eeeeoUla 
Utes beece eee eee COULDA D200 ccceesasveses cOCDA 


3.11.5 Nomenclature List 


DO depth at forward perpendicular, ft 
D1i0 depth at amidships, ft 

D20 depth at after perpendicular, ft 
DECKHT height of raised deck, ft 

DMB machinery box depth, ft 

FO freeboard at bow, ft 

Fi0 freeboard amidships, ft 

F20 freeboard at stern, ft 

H draft, ft 

LEN length between perpendiculars, ft 
RO (DO-D10)/LEN 

R20 (D20-D10)/LEN 

SO DO-D10 

S20 D20-Di0 
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3.12 Subroutine WEIGHT 


3.12.1 Introduction 

Subroutine VOLUME balanced the required and available 
volumes in the null. This subroutine estimates weights 
so that a balance can be made between full load weight 
and displacement. 

Several of the weights required have already been 
Calculated or were given as input. Calculated weights 
include the fuel weight and lub oil weight. All com- 
munications and control, group 4, weights and all ermament, 
group 7, weights are given as input. The weight groups 
noted refer to the NAVSHIPS Hull Group Weight Clessification 
of 1965. The categories in this classification scheme 
are given in Appendix E. This is the classification 
system now used by the Coast Guard. 

Also given as input are the weights of cargo, air- 
craft, ammunition and aircraft fuel. All other weights 
which are a part of the full load weight are estimated 
in subroutine WEIGHT. A flow chart for this routine is 
.8hown in Figure 30. 

There are few patrol cutters of recent vintage 
which can be used for collecting weight data. Three 
have been used throughout this model. These are the 
378° WHEC, the 210° WMEC, and the 150° WPB. Wherever 
possible additional data has been shown. Light ship 
weight data was taken from the weight breakdown sheets 


prepared for each design. Load weights were taken from 
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WEIGHT SUBROUTINE 
FLOW CHART 











REAL AND WIM © 0O169IS#CN 


WTttle JO.% 0013751 ¥ 



















COMM ON WT203s. Sharlene 
STATEMENTS Consens 


START WEIGHT 





WT2M = 256.3 


(7000- SHPE) 


WTzamM: , 
ES : 
256.3 SHPE > 7000. 


WTHZs of/ARwWTIM 


WT300 = OSFS#ELLW W203 : 





WT30/=, 01335 % ELKW 
WTIK =~ 0109 BLEW 
Vol.» ENCVOL //d0000, 
WTS 2-2.804 e¥DLH Es EZOL 
WT6L 2 ,00c14lIZENCVOL 
NAC NOFFIVE Po tHVEML 
WTéEI =, Do lSEmplachs. Sen D 
QeAe /) 5 ti. 


YES 
RAPTS = 
AiwyT lay 


BonTtwT: 1.4 





WTCIL= BOATHT + 2e26F 
WITEIV= Do 2g ACT 259650 
WIGS = VTE pwTEr + 

ATED wre 







AJO Con AJo 
S e Yes 2 


WIPER, JOSWANOFP ¢ 0737 
NCPOF  O2[QVHN ELE 

woPS: (,022%2+,08202 ed) WAC 
+.001aS wr OR WAC + [GOW AC 

APAALS 2 MESHIP PUT ChEWt 
WTPE euTAS THT WT EVR 
+ LTCALGO pLWTAMM OP EWwThe 4+ 
WAC EV 


NES 












WTQZs wTtMewTtO3 

WT POST TG) ew Ed pWTGS 
+wWTG7 

WTHIS O6LS mw DS 

WIG WTI NewTiv eet Ctwrs 


WLS MPs CLUTE4 WTS 2 eS 
WTO HGS 0 wh wer) wb Ad E 
WTCREWY* OFJ27 7 NAC 










RETURN 


EWD 


Figure 30 
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either load weight calculations or from inclining experiment 


datae | 

The first hull weight group is the hull structure, 
group 1. In order to estimate this weight group more 
accurately, three of the sub-groups are estimated sep- 
arately. The remainder of the sub-groups are lumped 
together as a single group. In the program thie group 
is referred to as WT1M. The three groups subtracted 
are group 111, superstructure; group 112, foundations 
for propulsion plant machinery; and group 113, foundations 
for auxiliaries and other equipment. These are identified 
in the program as WTi11, WT112, and WT113, respectively. 

Figure 31 is a plot of WT1M versus cubic number, 

The correlation here is very good. Also included on the 
figure is the estimating line used in the Navy's destroyer 
model, DDO7. 

Superstructure weight, WT111, is shown in Figure 32. 
plotted against deckhouse volume. The line labeled no 
uptakes is the one used in the program. The line used 
in DDO7 is also shown. This line assumes that the deck- 
house ayy eemtinum. 

The weight of propulsion plant foundations, WT1l1i2, 
is shown in Figure 33 plotted against the weight of 
propulsion machinery minus shafting and propellers. Again, 
& line used in DDO? is shown. 

WT113, other foundations, is shown in Figure 34 
versus the total weight of groups 3, 4, 5, and 7. The 
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WT2M = WTG2-WT203 


Figure 33 
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150 
correlation here appears reasonable; however, there is 


a difference between the data and the DDO? estimate, 

Propulsion machinery is classified under group 2. 

In this program, an estimate of the weight of sub-group 

203 (shafting, bearings, and propellers) was made separately. 
The remaining sub-groups were estimated as one group, WT2H. 
A plot of WT2M versus shaft horsepower for diesel plants 

is shown in Figure 35. Here more data is available fron 
non patrol types. The DDO7 estimate for diesel plants 

is poor. For CODOG plants it is assumed that the diegel 
plant will be 7000 horsepower or under and that only 

the endurance power will cause a difference in weight 

from that of a 378° WHEC. A lower endurance power will 
reduce the weight somewhat. 

Shafting and propeller weight, WT203, is calculated 
as weicht per foot of shaft length versus shaft horsepower, 
Twin shafts are assumed withatotal shaft length 0.56 
times the length of the ship. A plot of this relationship 
is shown in Figure 36. 

Group 3 includes the electrical plant weights. The 
Coast Guard intends to install electric plants of relative- 
ly larger capacity in future ships so a simple linear 
estimation using current designs is unrealistic. Instead, 
the weights of a number of diesel generating sets of 
various capacities have been plotted as shown in Figure 37. 
The varying RPM of these plants resulted in three weight 


ranges * 
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Figure 37 





154 
The weight of the generating sets is only a part 


of the group 3 weight. The weight of switchboards, WT301, 
is shown in Figure 38 as is the remaining group 3 weight, 
WT3R. The correlation for these weight groups is quite 
bad. Fortunately, these weights make up a small percentage 
of the total weight of the ship. 

Since group 4 weignts are given as input, group 5 
weights, auxiliary systems, is estimated next. The 
destroyer model uses the estimating relation shown in 
Figure 39. The three Coast Guard vessels agree very 
well with this relationship . 

Group 6, outfit and furnishings is the last light 
ship group estimated since group 7 weights are given as 
input. Group 6 is broken into four sub-groups for estimating 
purposes. The first, WI6I, includes groups 600, 602, 

603, and 605. These are hull fittings, rigging and canvas, 
ladders and gratings, and painting, respectively. The 

second group, WI6IIJ, includes groups 604, nonstructural 
bulkheads and doors; 606, deck covering; 607, hull insulation; 
and 608, storerooms, stowages and lockers. Group three, 
WT6III, includes only group 601, boats, boat stowage 

and handling. The final group, Wf6IV, includes all the 
remaining sub-groups of weight group 6. 

WT6I plotted against enclosed volume is show in 
Figure 40. wWT6II and WI6IV are plotted against number of 
accommodations as shown in Figure 41. WT6III is estimated 
in three parts. The first is a boat weight, the second a 
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liferaft weight and the third a life jacket weight. 


All of the lightship weights are summed and then 
multiplied by a design and builders margin. The design 
and builders margin is an input. 

Crew weight is taken to be 0.0737 tons per man or 
165 pounds per man. Personal effects are calculated 
as 0.105 tons per man for officers, .0737 for chiefs 
and .029 for enlisted men. 

The weight of personnel stores includes provisions, 
general stores, and potable water. The weight of provisions 
js based on Table 3. The data in this table was reduced 
to a formula for which only the endurance days of dry 
provisions need be given. This formula is: 

(.0222+.00202 x D) x NAC 

where D is the endurance days for dry provisions 


NAC is the number of accommodations 


Table 3 
days endurance pounds/person 
dry prov. frozen chilled total 
30 | 97.0 36.8 51.8 185.6 
4s 145.5 5502 52.8 25365 
60 194.0 736 53.8 321.4 
90 291.0 110.4 55.8 457.2 


The weight of general stores is taken to be 
000135 x D x NAC 

and that of potable water to be 
0186 x NAC or 50 gallons/man 
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The full load weight can then be calculated by adding 


Ms os alealoe os oe LOGRD) 
TD ei eseisisleisieuelc «+ ss eo LC GARY) 
DBMAR..oeceeceee eLC( WW) 
MENistie’s oc oo ces e eLOUBE) 
BU Whie ssc occ eee LC(DDD 
MICWOL. osceececeolv( Bp) 
DM sis sco oe o waieCeD 
Ue sgeis os + 0 00-00 LC (BED 
MTYPE ccc scseces el (CC) 
Onis nieces sees ol (hh) 
ME Tisisie oc cece ee oe LOC AM) 
WOPF. « cwisccccces LC CAR) 


the load weights to the light ship weight. 


3.12.2 Input List for Subroutine WEIGHT 


SHEES Mies cites 76 LO (CC) 
SHPMecccccccesseselO(CC) 
MAGEUL. o + ccs se oes LO (LL) 
WOERGO.¢2 ccs oe 6 « eLCieiiy) 
WITCH. ic 055-c6 sicioteiere oo Cua) 
WEAMMO. 00 0sc0ccees Le (Lb) 
WIFUEL. occccscccesLC(LL) 
WEG2s coe ccccececeeLO(WW) 
WEG. ceccccccececeLC(WW) 
WEG is csccsccse veo LCiW) 
WHEO. ccc cccsccce solu) 


*only if JOPT = 2 


3212.3 Statement Descriptions 
A detailed description of each statement will not 
be given for this routine wince the flow of the program 
follows the description given in the introduction and 
the relationships are those plotted. 
One area that will be discussed is WT6III. Here 
the number of 15 man liferafts required is first calculated. 


Q=eNAC/15.+1 e 
RAFTS=AINT(Q) 


where AINT is a function which drops the decimal fraction 
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from its argument. Boat weight is assumed to be 1.4 tons 
up to 25 men, 5 tons between 25 and 50 men and 10.1 tons 
above 50 men. This weight includes launching equipment. 


Life jacket weight is taken to be .002 tons per man. 


3.12.4 Output List for Subroutine WEIGHT 


WiaUID oniices cola eee « LCN) 
MISHEPs cles oecis ce « LONE) 
WGI 1 occ cc cienleiess. sLO CVV) 
WEE, cies cece ce oe LCUu) 
WMMN 3 a oieieloxave'e «eles caf CGNAV) 
Mls cosccccossecseelt Vy) 
WE208 «sss ciceicee ves ICV) 
ELEM Mais sso Sioleteisae es LOG) 
WABOO ss cilcccle ce ccieele vy) 
WRI OL viele s sleielelerses oe LOGY ) 
MAUI oo slooaeveneiere cierereeCl ave) 
MIO oes secede celeee LU (VY) 


WEGIT ceiieeecscs es oll OW) 
MTGLIT oc cceeccescelO(VV) 
WIOE View elelee cleo woiels LOM) 
WICREW. oc ccccseee sIC (LL) 
WEGIUr asic e sie so cle oo LO OND 
MIGZs « os cleee + ccc oc liClWN) 
WIGGiss oo cece ccs e es Le WM) 
MEGS. « seisisisisisie's bo SLCK WM) 
W BGC ets ciciciels ciceisic s LONE 
WEP Wittee cise c cl citeeriLe (LI) 
WES rccrersreiete 6 «ie sis leo Cel) 


* only if JOPT = 1 


3.12.5 Nomenclature List 
The definition of all variables is the same as 


given in the MAIN program nomenclature except for the 


followings 

BOATWT boat weight, tongs 

D same as ENDDAY in MAIN 
NAC number of accommodations 


Q dummy argument 





RAFTS 


SFCH 
SFCHI 
VE 
VOL 
WT111 
WT112 
WT113 
WTiM 
WT203 
WT2M 
WT300 
WT301 
WT3R 
WT6I 
WT6II 
WT6III 
WT61V 
WTFDS 
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number of liferafts 


same as RGEND in MAIN 


same as SFCHHP in HAIN 


same as SFCMHP in MAIN 


same as VEND in MAIN 


enclosed volume x 107? 


weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 


of 
of 
of 
of 
of 
of 
of 
of 
of 


see text 


gee text 


see text 


see text 


group 111, tons 
group 112, tons 
group 113, tons 
remainder of group i, tons 
group 203, tons 
remainder of group 2, tons 
group 300, tons 
group 301, tons 


remainder of group 3, tons 


sum of weight groups 3, 4, 5, and 7 
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3.13 Subroutine YVCG 


3.13.1 Introduction 

The vertical center of gravity of the ship is estimated 
in this routine. The designer has some latitude with 
regard to the location of the center of gravity; there- 
fore, the value calculated by this routine should be viewed 
as a feasible rather tnan as the best value. 

- -—fhe effect that the veg has on the design can be 
investigated by using different values of free surface 
correction. The free surface correction is added to the 
KG. Therefore, either a positive or a negative value 
can be used to represent a shift in the KG location. 

In order to determine feasible locations for the 
various centers, data for the 378° WHEC, the 210° WHEC, 
and the 150° WPB was used. This data is shown in Tables 
4 and 5. On these tables the vcgs of the various weight 
groups are shown. These centers are related to various 
hull depths by a factor K. The hull depth used is given 
in the right hand column of these tables. An average 
_ coefficient, K, was chosen and is given in the column 
labeled Coef. K. 

The variables to which the vegs are referred were 
chosen arbitrarily by observing which gave the most 
consistent set of coefficients among the three ships. 

some of the weight groups shown in the tables are 
really input items. These weight groups are included 


to assist the designer in choosing realistic values to 





Weight 


ATL 


Wr1ii 


dT112 


i113 


WT2ki 


WT203 


WT300 


W7301 


WTSII 


WT61II 


WIOLV 


WIGY 


WIG? 


V2RTICAL CENTERS OF GRAVITY 


373" 
WHEC 


15.78 
© 553 


37 01 
1.4 


D773 
20 


21.05 
© 794 


15.14 
oak 


4.02 
«52 


20.35 
-768 


ele 
» 804 


Aa) 
oe o 


19.43 
© 733 


ZOnCo 
78 


26.48 
1.0 


41.75 
1.575 


24.88 
© 939 


Bio 
1.191 


Doece 
ES) 


210° 
WhieC 


14.72 
oe 


Dee 
1.26 


uu 
0157 


Ono 
wo) 


10.40 
eo 


6.0 
0214 


14.7 
s2) 


11.8 
e421 


19.41 
693 


16.638 


© 596 


20.43 
°73 


21.93 
783 


33.7 
1.738 


2057 
0749 


350 
1.25 


24.45 
873 


Table 4 


0. 
WPB 


955 
© 584 


23.6 
1.47 


25 
0156 


15.1 
S44 


8.50 
0531 


2.6 

0163 
19.0 
1.19 


21.0 
1.31 


17.79 
oe 


13.5 
, B44 


15.09 
943 


19.23 
ac 


2509 
oe 


13.46 
o 841 


22.6 
1.41 


21.3 
1.33 


Coef. 
K 


055 


1.4 


a2 


s57 


0 54 


016 


i. 3 


1.3 
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times 


DAVG 


D1 0 


D1 0 


Center groups 


Stat tay 


D1 OMN 


D10 


D1iO 


D10 


D1 0 


D10 


D10 


D10 


D1 OMN 


D1 0 


D1i0O 


D10 





Weight 


Crew & 


veg 


Fersonal K 


sifects 


Provisions vcg 


General 
Stores 


Ammo 


Potable 


Water 


L.U. 


ruel 


Aarcrait 


K 


vcg 
K 


veg 
K 


vcg 
K 


veg 
K 


veg 
K 


VCe 


A/c Fuel veg 


DAYG 
Dido 
D1 Ow:N 


VERTICAL CANTERS Or GRAVITY 


B40) 
WHEC 


26426 
«964 


15.0 
oe 


24.05 
° 883 


15.99 
587 


7.06 
«26 


3.58 
shod 


7.67 
eon 


39.01 
1.43 


9.56 
0351 


28.54 
27225 
27625 


210° 
Wik C 


20.43 

1,048 

By fe 
2 536 


16216 
«829 


8.4 
0431 


8.5 
2 30 


Zao 
o12é 


5.17 
«292 


26.91 
Z2oa0 
19.5 


Table 5 


50 
WPB 


VS a0 
0938 

11.0 
«688 


11.0 
688 


19.0 
1.19 


4.2 
© 263 


12.0 
0°75 


6.6 
e412 


16.36 
16.0 
16.0 


Coef. 
K 


ait 


1.4 


oon 
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times 


D1 OMN 


Di 0 


D1 OMN 


D1 OMN 


D10 


D1 OMN 


D1 0 


D1 0 


D10 





START NCG 
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Figure 42a 
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WIT203)/W TEV 











CGCREW= .G7KDIOMN 
CGPE = CGCREw 
CGPW2 .26*+DdD/0 
CGG3= .GSxd/O MD 
CGPV = ,SS*DIO 











CGP S52C, OF22 +,002% Dd) 
C GPY + 0012 CD¥Y CGGS + 
ou err °2982+.0933 
Db 







+OGhl xy wraAc + 
CACEL# WAC FUL) 


Cane 2COLSHI be WLSwiP + 
LOMOMN)/WEULLD 
CFULD = CFVLLO 
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input. 


A flow chart for this routine is shown in Figures 
42a and 42b. 


3.13.2 Input List for Subroutine VCG 


Oe Crm. « » epencaiee lc (EE) 
AOvUlhs + os es 6 ove LC (EE) 
MNOS e's clelereieis eee LOer..) 
BOW ices. sisisvecraO(EE ) 
PMORUIG s+ sisiere ¢ au0.0 LO UHH) 
RG OC cisisrsis « eicueiellC (EE) 
CARGOW, 6. .00005eLC(ES) 
Dieriiela¥eia sis’ siayara LOK) 
DUMONT, oo 0000s 00 cLe (IL) 
ON ica terete ss aieie B ) 
BOWARs 00000000 0LC( WR) 
GR4CG..ccceceeee LC (EE) 
GREW, cc ceseceske (BE) 
PEO ccs bees e sLOCES) 
CRYO. occ cece de LOtEe) 
TO iels ccc cie sees eG (5B) 
EMM eisis cee s cs oso SUCKS ) 
WOAEON . 0a cccc ves LUFF) 
PROENUD ss so ee oe eHOUET) 
MBG AU Ios si eheltie sores Olay) 
POARGO.. «cco eBOCED) 
WFULLD.. 12 eves ell (WH) 


MIMS corel ies ¢ oes EC (VN) 
WISN? os since ss 0s LC(VV) 
WT203.cccvcccceeelC(VV) 
WE2H. occcccccceceLC(V¥) 
BT IO0 <ccieialoleseo.s. 0 LC CVV) 
WT301..ccccecees LC(VV) 
WEDR s sane wysee LCV) 
Wile e cee ce cLOCUN D 
MECUT is. stl otk CCVND) 
WL6IIIeececeeeseeLC(VV) 
WI6IV. ccccccceeeeLC(VV) 
WEAKMO. secceceee eLC(LL) 
WICREW. oc ccceceeeLC(LL) 
WEUEL.a sees cess LCCLL) 
WTGL wes cecceeee LOCWW) 
WIG2. cece cc cove LOEW) 
WHGS. .. sce s aes. LOE) 
WiGiecccess ccee c elem) 
WIGS. .deee. «scowl (WW) 
WIG6 se eccesece ce oLC(WW) 
WIGC7 seve ccccesce LE(WE) 
Mawes. ccs le (LL) 
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ISTHE ¢ occ cle os see's eo LCCRW) Ween et eras anedtereie.eie LC Llu) 


WE111ecccccccccvceeLC(VV¥) W2PScecccccceeecLC(LL) 
MIM 2 osc 0 oe eee ee LO) 


3.13.3 Statement Descriptions 

Detailed statement descriptions will not be given for 
this routine. In general the statements consist of vcg 
estimates using the relationships mentioned in the intro- 
duction or of calculations to combine weight groups by 
summing the vertical moments of individual groups and 
dividing by the total weight. 

Note that DiO is different from D1OHMN only if the 
length of raised deck is greater than one half the ship 
length. It must also be remembered that input moment 
arms are given as fractions of the main deck depth at 


amidships and not as a distance in feet. 


3.13.4 Output List for Subroutine VCG 


CCARGO. se eeeeeeeeLC(00) CGGsdes ssc LC( 00) 
OEULDs.....eeeueee SODA CGGH. ..000e00018(00) 
SAUL Dy). /-/e\cle ss)Jeeta (OO) WESigoas soqoonverioes 
BGiCis. .. seve ss SOLO CHH) CCEOMe. + eee LC (00) 
CGAMMO.....00e0eeLC(HH) Caan dees). sEC(O0) 
CGCREW. ...00050+eLC(HH) CGWOMs...>.++eLC (HH) 
MOP UE Lise vesscee s ee lot Hi) CGPies 6s sle oe se LO (HH) 
EGC... +s seers Le (OO) GGPSW s+. 0. +++ -LCC(HH) 


BEGOh... >» eeeeerel(00) CLSHIP...++++eLC(00) 





Bel3e5 
ACMOM 
AMOMOM 
CACFL 
CFULD 
CFULLD 
CG111 
CG112 
CG113 
CGiM 
CG203 
CG2M 
CG300 
CG301 
CG3R 
CGé6I 
CG6II 
CG6III 
CG6IV 
CGG1 
CGG2 
CGG3 
CGG4 
CGG5 
CGG6 
CGG7 


CGGS 


Nomenclature List 
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total vertical aircraft moment, tons 


total vertical ammunition moment, tons 


vertical 


center 


same as CFULLD 


vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 
vertical 


vertical 


center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 
center 


center 


aircraft fuel, ft 


of 
of 
of 
of 
on, 
of 
of 
of 
of 
of 
of 
of 
of 
of 
of 


of 
of 
of 
of 
of 
of 


ship at full load, ft 


group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 


111, ft 
iterate 
Get 
iM, ft 
203, ft 
2M, ft 
300, ft 
301, ft 
3R, ft 
61, ft 
6II, ft 
ie emer 
6IV, £t 
1, ft 

2, ft 

3, ft 

4, ft 

5, ft 

6, ft 
7, ft 


general stores, ft 





CGOMOM 
CGPV 
CGPW 
CLSHIP 
CPO 
D 

DO 
Di 0 
Di OMN 
D20 
ENL 
G4bi Oni 
G7MODI 
HLEN 
I 

OFF 
WTiM 
WT111 
Wi11i2 
WT113 
WT203 
WT2M 
WT300 
WT301 
WT3R 
WT6I 
WT6IT 
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vertical moment of cargo, tons 


vertical center of provisions, ft 
vertical center of potable water, ft 
vertical center of light ship, ft 
same as NCPO in MAIN 

same as ENDDAY in HAIN 

depth at F.P., ft 

depth amidships to main or raised deck, ft 
depth amidships to main deck, ft 
d@épth at &.P., ft 

same as NENL in MAIN 

total group 4 moment, tons 

total group 7 moment, tons 

half the length of the ship, ft 
index 


same as NOFF in MAIN 


weight 
weight 
weight 
welgnt 
weight 
weight 
weight 
weight 
weight 
weight 
weight 


group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 


iM, tons 
111, tons 
112, tons 
113, tons 
203, tons 
2M, tons 
300, tons 
301, tons 
3R, tons 
61, tone 
6II, tons 





172 
WT61II weight group 6III, tons 


WI6IV weight group 6IV, tons 
XLDMOM total load moment, ton-ft 





173 
3.14 Subroutine COST 


3.14.1 Introduction 

This subroutine is based exclusively on the pre- 
liminary cost estimating procedure currently used by the 
Coast Guard. To use the procedure, the weight of each 
of the seven light ship weight groups must be known. 
Linear graphs for man hours and material cost are then 
entered with these values to determine the cost and labor 
required for each weight group. 

Various relationships are also given for determining 
design costs, construction services cost and miscellaneous 
costs. Miscellaneous costs include initial outfit equip- 
ment such as life jackets as well as retrofit costs, spare 
parts and administrative or resident inspector's office 
costs. 

The curves used in the procedure were converted 
to formulas and the procedure was programned as subroutine 
COST. A flow chart for this routine is given in Figure 
43, 

All of the material cost curves are based on 1959 
prices. Indices are used to correct the prices to the 
year the ship will be built. These indices are based 
on Bureau of Labor Statistics data for the particular 
weight group. 

The labor rate for all of the weight groups is assumed 
to be the same. The labor rate is inputed as an average 


dollars per hour for the year the ship will be built. 
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Material cost indices and labor rate data are shown 


in Figures 44 and 45. There is a large fluctuation in 
costs between shipyards so these numbers can only be 
viewed as average values. 

In addition to the corrections for inflation, eighty 
percent is added to labor cost for overhead and ten 
percent is added to ali costs for profit. 

Initial outfit costs are assumed to be two percent 
of the total cost. Spare parts cost is the sum of the 
cost of spares for propulsion machinery, electrical 
machinery, electronics and outfit. The cost of spares 
as a percentage of the weight group material cost is 
given below: 

Group 2 --- 9% 

Group 3 --- 8% 

Group 4 --- 25% 

Group 5 --- 10% 

Retrofit cost is assumed to be 4 percent of the 
total cost. Administrative socts are assumed to be 3.5 


percent of the total. 


3.14.2 Input List for Subroutine COST 


CSTC2. ec ccees ee LC(MB) GUIND. «0 0c00eas eLOCNN) 
CSTG4 oe cece reese e LO(MM) G2IND..cccocceeslO(NN) 
CONG] soc ccc cee LOCH) G3IND...0.000066LC(NN) 
DBMAR..scccceeceLC(WW) G5INDE.. 22 se+ os sLC(NN) 


DOLHResescoeececLC(NN) G6INDs oecsccee ee LC(NN) 





INDEX 


1.7 


1,6 


15 
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1,2 
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BOPT. a. 00000. LC(CC) WIG3. cccceccces LC (WW) 
TEI) ocoicle tele 6 «.0, 2/0 6 atHXESCDA WIG. oecccccceeoLC (WW) 
WRB OREE EE ois aie,s0ss ss LC (CC) WIGS. ccc ceseeeeoLC (WW) 
BIG ooo bese eee os LOGHW) WLG6occcccccccecLC( WW) 
MEG?» occ cc secs oe sLOCWH) WIG7 sc cccccsseeoLC (WW) 


fonly if JOPT = 2 


3.14.3 Statement Descriptions 
Wi=WTG1 *DBMAR 


W7=WTG? *DBMAR 
The design and builders margin is evenly divided 
among the weight groups for costing purposes. 


GIMH=10.#*(ALOG10(W1 )*.75217+2.97276) 


G7MH=10.**(ALOG10(¥W7)*.75217+2.97276) 

These formulas represent the relationships between 
man-hours and the weight of individual weight groups. 
These relationships are assumed linear when plotted on 
log-log paper. 


GiMC=10.**(ALOG10(W1 )*.9542434+2. 790484 ) 
G2HMC=CSTG2 
IF( JOPT.EQ.2) GOTO 1 
G2HC=10.**(ALOG10(W2 )#1.018885+3.634328) 
IF(MTYPE.EQ.2) G2hiCeG2MC+1500000. 

1 G3MCH10.*#*(ALOG10(W3)*1.0725514+3.667812) 
G4MC=sCSTG4 
G5MC 210. **( ALOG10(W5)*.870906+3.974192) 
G6MC=10. **(ALOG10(W6) #1. 068265+3. 603833) 
G7MC2CSTG7 


These formulas represent the material cost curves 


used by the Coast Guard. Some of these costs, CSTG2, 
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cSTG4Y and CSTG7, are input. Group 2 cost is an input 


only if the machinery plant is specified. Otherwise, 

the machinery cost is calculated. The formula used 

represents the cost of diesel plants. An extra $1.5 

million is added if a CODOG plant is to be installed. 
G1LC=G1MH*DOLHR*1.98 


G7LC=G7MH*DOLHR#1.98 

The labor cost by weight group is calculated here. 
The value, 1.98, represents an overhead factor of 1.8 
times a profit factor of 1.1. 

GiNC=sGiIMC#GLIND#i.1 

IF( JOPT.EQ.2) GOTO 2 

G2hiC=G2hic*G2IND 

2 G3MC=G3MC*#G3IND*1.1 

G5NC=#G5NC FGSIND#1.1 

G6MC=G6MC #GSIND#1.1 

G2hiC=G2iic *1.1 

GUMC2#G4NC*1.1 

G7MC=G7iMC#1.1 

In these formulas the material cost is corrected 
for inflation and a factor of 1.1 is included for profit. 
Note that the costs given as input are not corrected 
for inflation. Therefore, the cost of the equipment 
must be corrected for inflation before being inputed, 


LABCST=G1LC+G2LC0+G3LC+G4LC+G5LC+G61LC+G7LC 
MATCST=G1MC+G2iC+G 331C+G41iC+G 56C+G6MC+G7MC 


The total labor and material costs are calculated 


in these statements. 
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LOA=1.08*LEN 

FACTOR=(.114*LOA-1,6)/100. 

DLO=sFACTOR*LABCST 

FACTOR=(-.01*LOA+7.5)/100. 

DMC=FACTOR*MATCST 

These statements calculate the design labor, DLC, 
and design materials, DMC, costs for the first ship of 
a class from the beginning of the conceptual design stage. 
The two values of FACTOR are based on Coast Guard curves. 
Since the length overall is required in these curves, its 
value is estimated as eight percent greater than the 
length between perpendiculars. 

PACTOR=(.014*LOA+7. )/100. 

CLCeFACTOR*LABCST 

FACTOR=(-.015*L0A+6.5)/100. 

ChCa2FACTOR*HATCST 

The construction services labor and material costs 
are calculated here. Again the formulas for the two values 
of FACTOR are based on Coast Guard curves. 

TOTCST=LABCS T+MATCST+DLC+DiIC+CLC+CHC 

ORICST=.02*TOTCST 

SPARESa#(.09*G2MC+.08*G3MC+.25*G4MC+.1¥*G5MC )/1.1 

RETRO=.04*TOTCST 

ADMIN=.035*TOTCST 

The four items that are included in the miscellaneous 
cost are calculated here. Since profit is included in 
the individual weight group material costs it is divided 
Out when calculating the cost of spares. 


CSTGizG1LC+G1biC 


CSTG7=G7LC+G7MC 
Total weight group costs are calculated here 
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DOS T=DLC+ DMC 
COSTsCL0+Chic 
MCSTsORICST+S PARES+RETRO+ADNIN 
TOTCST=TOTCST+HCST 
RETURN 
END 


Finally, the total design cost, construction services 
cost, miscellaneous cost and the total first ship cost are 


calculated. 


3.14.4 Output List for Subroutine COST 


BCC To. oo os ste ALCOR) CSTC6. ccc ceec cece o LC (ME) 
CSTGIG. .. « clea CM) CMG 7 sac as-oe steer LC (Sy) 
GSPEZ. .. so sete ae LOU) HOST. Sees... eee LCC ) 
BCMCS es ss +s cleisisine LOOM) MEST. cilec ooo ose s ce LC CAN) 
CSTE ane oe oes uenee ¢ LCR) FOTCS TEN... 00 ee LC (NN) 


WOMGDs soe oe alelelele « LORMM) 
* only if JOPT = 1 


3.14.5 Nomenclature List 
The definition of all variables is the same as that 


given in the MAIN program nomenclature except for the 


following: 

ADMIN administrative costs, dollars 

CLC construction services labor cost, dollars 
CMC construction services material cost, dollars 
DLC design labor cost, dollars 

DMC design materials cost, dollars 

FACTOR see text 


G1LC weight group 1 labor cost, dollars 





GZ2LC 
G3LC 
G4LC 
G5LC 
G6LC 
G7LC 
GiMC 
G2MC 
G3MC 
G4MC 
G5MC 
G6MC 
G7iC 
GiMH 
G2MH 
G3MH 
G4 MH 
G5MH 
G6MH 
G7MH 


LABCST 


LOA 


MATCST 
ORICST 


RE 
RETRO 
SFCH 


weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 
weight 


total construction labor cost, dollars 


group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 
group 


7 


labor 
labor 
labor 
labor 
labor 


labor 


cost, 
cost, 
cost, 
cost, 
cost, 


cost, 


dollars 
dollars 
dollars 
dollars 
dollars 


dollars 


material 
material 
naterial 
material 
material 
material 
material 
man-hours 
man-hours 
man-hours 
man-hours 
man-hours 
man-hours 


man-hours 


length overall, ft 


cost, 
cost, 
COST, 
cost, 
cost, 
cost, 


cost, 


dollars 
dollars 
dollars 
dollars 
dollars 
dollars 


dollars 
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total construction material cost, dollars 


initial outfit cost, dollars 


same as RGEND in MAIN 


retrofit cost, dollars 


same as SFCHHP in MAIN 


if 
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SFCN same aS SFCMHP in MAIN 

SPARES spare parts cost, dollars 

VE same as VEND in MAIN 

Wi weight group 1 with margin, tons 
W2 weight group 2 with margin, tons 
W3 weight group 3 with margin, tons 
wh weight group 4 with margin, tons 
W5 weight group 5 with margin, tons 
W6 weight group 6 with margin, tons 
W7 weight group 7 with margin, tons. 





3.15 Subroutine SEASPD 
3.15.1 Introduction 

This routine is based on the work of Michael Shen, 
Peter Gale and Kenneth Walker as reported in References 
4% and 5. In these references is given a computer program 
for computing the average sea speed. This computer 
program was altered by dropping the option to use a bow 
sonar dome and by deleting the read and write statements. 
The program was also rewritten in ASA standard FORTRAN. 
Except for these changes, subroutine SEASPD is the computer 
program developed by these three men. The interested 
reader is referred to these references for the dovelopment 
of the routine. 

Briefly, the routine assumes: (1) that the ship 
is operating in the North Atlantic Ocean over a long 
period of time; (2) that it is steaming into head seas 
one half of the time; (3) that the remainder of the time 
(when steaming at other headings) it is not forced to 
reduce speed due to ship motions. Part of the time the 
ship is steaming into head seas, its speed will be limited 
be the sea state. In the routine this is simplified to | 
&@ single sea state limited speed. The remainder of the 
time the ship is assumed to be able to make its maximum 
sustained speed. The percentage of the time the ship 
is limited in speed is estimated based on average wave 
heights in the North Atlantic Ocean and on ship character- 
istics. 





185 
The average speed at sea is calculated assuming 


that the vessel will travel at its maximum sustained 
speed whenever it is not limited in speed by the sea 
state. When the vessel is steaming at headings other 
than into head seas, its speed is not limited. This is 
assumed to be one half the time. When steaming into 
head seas its speed will be limited only part of the 
time. Both the limiting spsed and the percentage of 
time that it applies are calculated in SEASPD. 

This routine has deen included to indicate ths 
effect that changes in the vessels characteristics will 
have on its seakeeping ability. The many assumptions 
made in deriving the average sea speed make the speed 
predicted practically meaningless. However, the speed 
calculated can be used to compare two designs to decide 
which will give the better sea keeping performance. | 

A flow chart for this routine is shown in Figure 46. 


3215.2 Input List for Subroutine SEASPD 


Be wtshnickebetaieinai loi 0GG') DPTRY. 0 0scccee se cl (GG) 
MaIMio1s 320. sieeserera OC (GG) Haute sien since LC GG) 
MED. .isie's's atic LC (GG) OE octenerenicre o's oe c LC CBE) 
GAGs eae 0s eae Le (GG) VM Gihcsrocaiameie ROC’) 


5.15.3 Statement Descriptions 
For a detailed discussion of this routine soe Refer- 


ences 4 and 53 





SEASPD 


FLOW -CHART 


REAL AND 
LABELED 


Lomm or 
STRTEMEWTS 


START SEASPD 


CXPs9 
CX Ha (CX, 314)#,27/93 
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FKYY= .G0S/ 
CPP 2 (CP-,5g)¢2. 76666 


CSP = FKYY 2 CPOP+CXP 
Be im SO Oe 
wLbn «be He (LED + 


TPs 2ePte (Clipz wire 


YP = (Cx- 84) ¥ 8139S 
Cre (Cx~.04 "8.60405 











CM: eee eee 


THe G.288/ Bereta 
RaAnwP) ) 


(a= 
TPs re/(sanr (wen7) 
THUs THAR (CL bY)) 


T= TPL 





exir/4 +(6//C-Z, 
Fke= 599 
P= (CP-,S8)#/.9/666 
CSH sFKRAtCAH+ CXH 
SmM= 22994 OPTLY/G 
AwWwlh=CwPlLew #8 














VLEtVRTHVLE 
AO SPEED= Vie HS aeriien 








Yes 


PELU= 00.- PERI) /2. 
FlR.=1/00-.-FERU 

Avse=(Pcarususe 
PERM SPEEDY / 09, 


SUBROUTINE 
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KY< 24 
G: 32,/5S 
Ral.9Ios 
Pr-= 3.14159 


Ly (KY eH LEW) At 262740, * 
DPRK /G 
CIL= eo /421O Sr CHP ~_ INS 





TV = CIL# Cw LEVDKES 
c-H 

g1= 6/2 

A= LEMA 

Bcc o $(bi/< , Yi) 
YYKic FC 8Cc) 








THRE 








( .9S#@ 


S(2,# PI ))Y/1. ©3878 
W/E He OLS CAESHLEN? 


Wo DPT LY/Coom(srten) nd) 
ww = fiw * "425233. 9 6439 


VLO<(90254-- 


Meee: 


VR =/6+Hh ew) +ecr91 
ALGO ew WK 


ELI =. 
f, (waved) 









Figure 46 , 
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301524; Output for Subroutine SEASPD 


y.10, whe; only; output is the value of average sea speed, 
AVSP. 


3.15.5 Nomenclature List 
P The definitions used are the same as those given in 


References 4 and 5 except for the following: 


Value: used . Value used in 
in SEASPD references 
- °C BEAM 
LEN SL 
w0|UC°6COC CW 
— B 
Ree Nn te, = » 
DPTRY DISPO 
— DKY 
(a YYI 
Megorkh? vt stete 

YJ 
Peeve rt 
LY YI 
VSUS © VVV 
VLR | VPX 
) Ait ‘ 
Op Cor ane a osa 
ores f IHes {<tc . 4 
en oe e 
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3.16 Subroutine OUTPUT 


3.16.1 Introduction 

All program output is printed by subroutine OUTPUT 
except for error messages. A simplified flow chart for 
this routine is shown in Figure 47. Subroutine OUTPUT 
js divided into two main parts. The first prints the 
input data and the secund prints the output generated. 

If an error nessage is printed, only part of the 
output data is available and the data that is available 
is probably not correct. Therefore, when an error message 
occurs, only the input data is printed. Variable I1 has 
avalue of 1 in this case. If the program successfully 
completes all preceeding steps, I1 is assigned a value 
of 2 and output data is also printed. 

Subroutine OUTPUT consists primarily of write and 
format statements. Some calculations are made but only 
to convert previously calculated values to a desired 
output form such as determining ratios between weight 
groups and full load weight. 

Only two decision points are used. One is the decision 
on whether or not to print output. The other checks the ~ = 
option used for inputting vehicle performance data and 
chooses the correct set of write and format statements 


to print the input data. 


3.16.2 Input for Subroutine OUTPUT 
All values which are printed are inputs to this 


routine. 
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OUTPUT SUBROUTINE 
FLOW CHART 


REAL AND 







sT + OuTPUT LABELED 
ART O Common 
STATEMENTS 





VO 
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3.16.3 Statement Descriptions 


Detailed descriptions of the statements of this routine 
will not be given since most of the statements are self- 
explanatory. Use of the sample outputs given in Appendix 
B will be of use in understanding the format statements 


used. 


3.16.4 Output for Subroutine OUTPUT 
No values are generated in this routine for use 


elsewhere. 


3216.5 Nomenclature 
All the variables used in this routine have been 
defined previously except for some dumny variables whose 


definition should be clear from context. 





CHAPTER IV 
EVALUATION OF THE PROGRAM 


A number of test runs were made on an IBM 370, model 
168, computer located in the Information Processing Center 
at M.I.T. The test runs were made using input data for 
four existing designs, the 378° WHEC, the 210° WHEC, the 
150° WPB, and the Navy's Patrol Frigate design. The Cutter 
Model produces feasible designs for each of these ships 
except the 150° WPB. 

The design for the 150° WPB failed because a raised 
deck length greater than the ship length was required. 
The actual ship does not have a raised deck. The probable 
cause of this failure lies with the available area and 
volume relations used in the program. These relations 
were developed for ships with lengths greater than 400 
feet and probably have significant errors when extended 
to ships as short as 150 feet. Recommendations for cor- 
recting this error are given in Chapter V. 

A comparison between the results produced by the 
Cutter Model and the actual ships are shown in Tables 
6a and 6b for the 378° WHEC and in Tables 7a and 7b for 
the 210° WMEC. The differences are reasonably small in 
most areas; however, a few do need discussion. 

The largest difference with both ships occurs in 
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Table 6a 
Comparison of Cutter Hodel Output to the Actual 
Ship for a 378' WHEC 
lpp Beam Draft DO Dio D20 Davg Disp. 
t ft ft ft ft ft ft tons 


Ship 350. 42.0 14.3 34.75 27.25 27.5 28,54 2968.0 
Model 350. 43.32 14.6 38.33 27.83 28.18 29.64% 3061.5 














WEIGHTS VCG 

Ship Hodel Ship Model 
WTG1 917.46 962.8 19.12 19.42 
WTG2 383,6 bp fad: 11.45 11.52 
WTG3 88.12 23049 21076 
WTG4 80.2 80, ; 31.56 S2ee 
WTGS5 268.8 262.3 19.43 19.76 
WTG6 325% 322.3 25.63 25.88 
WTG7 39.0 39.90 35.26 36.30 
WTLSHIP 2102.2 2215.0 19.72 19.98 
LUB OIL T3957 16.0 3.58 3.62 
CREW 13.7 27200 
PER EFFECTS 3 22-03 7.8 326.26 27.00 
PER STORES 107.4 95-1 ie o> 14.36 
CARGO O. 0. 
ALIMO Se 36> 15.99 16.42 
AIRCRAFT 3.4 304 39.01 40.08 
A/C FUEL 19.6 19.6 9.56 9.74 
WTFLLD 2967.6 3061.4 16.61 17.20 
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Table 6b 


Comparison of Cutter Model Output to the Actual 
Ship for a 378* WHEC 


* 460 sq ft included in input area 


** Excess area was used for crew rec, rooms 


ship Model 

Hull Volume 277600 273734 cu ft 
Deckhouse Volume 95611 94705 cu ft 
Total Arrangements Area 30863 29664 sq ft 
Hull Arrangements Area 19615 18522 sq ft 
Deckhouse Arrangements Area 11248 11142 §8sq ft 
Office Space 493 S10. #satt 
essing Facilities 2713 2713 sq ft 
Crew Special 725 787 sq ft 
Officer S.R. 1157 1168 sq ft 
Officer Sanitary 243 192 sq ft 
CPO S.R. 1040 1040 sq ft 
CPO Sanitary 250 208 sq ft 
Crew Berthing 3641 Bieoe sq tt 
Crew Sanitary 716 720 3=6sq ft 
CO S.R., Cabin & Pentry 755 753 eo ft 
Commissary Stores 600 652 sq ft 
Other Stores 3015 3231 sq ft 
Workshops 1161 1254 sq ft 
Repair Lockers 148 100 sq ft 
Chain Lockers 47 40. sq ft 
Uptakes 462 0 sq ft 
Steering Gear & Windlass 736 779 +‘sq ft 
A/C & Fan Spaces 369 396 §=6sq ft 
I.C. Spaces 312 312 sq ft 
Aux. Machinery Spaces 1860 2180 sq ft 
Pilot house, Chartroom & CIC 1125 1111 sq ft 
Passages 2567 2450 sq ft 
Input Area 5283 a8 5900 sq ft 
-Excess Area 1445 0 
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Table 7a 


Comparison of Cutter Model Output to the Actual 
Ship for a 210° WHEC 














Lyp Beam Draft DO Dio D20 Davg Disp. 
ft ft ta ft ft £t rt tons 
Model 200. 35.02 8.98 21.50 19.50 19.70 24.97 959.1 
WEIGHTS VCG 
Ship Model Ship Model 
WIG1 379.5 347.8 16.4 13.92 
WTG2 102.6 102.6 9.5 9.55 
WTG3 43.7 7504 17.5 22228 
WTG4 26.1 26ne 34.2 33.31 
WTG5 110.4 109.6 16.6 19.88 
WTG6 104.4 106.8 22.6 25.93 
¥TG? 8.5. &.5 23.4% 2seeal 
WTLSHIP 77502 776.8 7 1 17.40 
FUEL u52.3 137.3 8.17 8.40 
LUB OIL 3.8 3.9 5 2454 
CREW $8 6 5.5 20.43 18.92 
FER EFFECTS ; 32 20.43 18.92 
PER STORES 47.3 30.4 10.99 11.95 
CARGO GO. QO. 
AMMO Za 2.6 8.4 8.00 
AIRCRAFT 0. 0. 
A/C FUEL 0. 0. 
WTFLLD 990.0 959.5 15.38 15.87 
Ship Model 


Raised Deck Length 164, 120. 





Table 7b 
Comparison of Cutter Model Output to the Actual 


Ship for a 210° WHEC 


Hull Volume 
Deckhouse Yolure 


Total Arrangements Aréa 
Hull Arrangements Area 
Deckhouse Arrangements Area 


Office Space 

Messing Facilities 

Crew Special 

Officer S.R. 

Officer Sanitary 

CPO S.R. 

CPO Sanitary 

Crew Berthing 

Crew Sanitary 

COS.R., Cabin & Pantry 
Commissary Stores 

Other Stores 

Workshops 

Repair Lockers 

Chain Lockers 

Uptakes 

steering Gear & Windlass 
A/C & Fan Spaces 

I.C. Spaces 

Aux. Machinery Spaces 
Pilot house, Chartroom & CIC 
Passages 

Input Area 

Excess Area 


Ship 


124700 
16100 


11242 


oot) 
1897 


130 
1452 
324 
567 
193 
283 
75 
1658 
84 


06 
480 
1294 
Sk4 


110 
2 
0 
270 
165 
180 
820 
422 
(oo 
704 
0 


247 
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weight group 3, electric plant. This difference results 


from the fact that greater generating capacity relative 
to the average load has been built into the model than is 
installed in existing ships. This increase in generating 
capacity reflects current Coast Guard policy for new 
designs and is not an error in the progran. 

The differences in weight group 1 reflect the fact 
that the hull dimensions are not exactly the same for 
the program designs and the actual ships. Weight group 
1 is a function primarily of the hull cubic number. 

The fucl weight is calculated in the program based 
on specific fuel consumption, endurance range and endurance 
spsed. An attempt has been made to approximate the actual 
fuel weight by varying one or more of these quantities 
but an exact match is difficult to obtain. 

The only other significant waight difference occurs 
in personnel stores which includes potable water. The 
mode] results are less than the weights for the actual 
ships due to the fact that the amount of potable water 
per man has been reduced from the amount used on past 
designs. This is in accordance with current Coast Guard 
policy for future designs. 

The vertical centers of gravity show quite a bit 
of variation when compared to the actual ships. However, 
the values in the program are only average values and 
will in general not agree with the values for any one 
Bhip. 
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In the case of the 210° WMEC, a raised deck is 


required. The model predicted a raised deck length 44 

feet shorter than that of the actual ship. This is another 
indication that the available area and volume relations 

are inaccurate for smaller ships. 

Hull areas are given in Tables 6b and 7b. The areas 
used in the cutter model for berthing and sanitary spaces 
are based on current minimum Coast Guard standards. These 
areas in gencral are lower than the areas on current 
ships. | 

Several runs were also made using data for the Navy‘s 
Patrol Frigate design. At the same time, runs were Rade on 
the Navy‘s destroyer model, DDO7, using similar data. 

The results of these runs are shown in fables 8 and 9. 

The Patrol Frigate design is of interest because 
it is at the lower limit of applicability for the des- 
troyer model and at the upper limit of applicability for 
the cutter model. Agreement between the values calculated 
by the two models and the actual ship as shown in Table 
8, is reasonably good. Again, the group 3 weight is highest 
for the cutter model because 3 - 2000 KW generators are 
used instead of the 3 - 1000 KW generators on the actual 
Ship. Also, it appears that the greater emphasis placed 
on habitability by the Coast Guard is reflected in a 
greater group 6 weight. 

Table 9 gives the displacements predicted for greater 
lengths than the actual ship length. As the ship length 





Comparison of Cutter Model Output to the Actual 


Table 8 
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Ship and to DDO7 Estimate for the Patral Frigate Design 


Beam Draft 


z 
ree tt 


Ship 408. 45.0 


Model 405. 
DDO7 405. 46.72 15.3 41.00 33.99 
WEIGHTS 

ship Model DD07 
WTG1 1222.36 1199.6 LSaiee 
WTG2 25601 258. 240.0 
WTG3 155.36 2ilee 156.6 
WTG4 92.58 70-7 TO. 7 
WTG5 351.3 317.5 381.3 
WTG6 268.82 333.8 238.0 
WTG7 95.68 8547 85.7 
WTLSHIP 2442.20 2637.4 2666.3 
FUEL 830.0 eplisl AY, 845.3 
LUB OIL 20.22 20.2 0. 
CREW & PE 20.96 20.9 20.9 
PER STORES 68.05 66.6 84.0 
CARGO 67.0 7520 80.9 
AMMO 58.8 58.8 Sono 
AIRCRAFT LOn0 10.0 10.0 
A/C FUEL. 0, 0. 0. 
WTFLLD 3517. 3707.7 3766.9 


ft 
14.5 








Dio 
ft 


45.95 15.34 41.07 28.92 29.33 31.01 
34.80 35.30 3766.3 








Disp. 
tons 


3517. 
3706-1 





Full Load Displacement vs. Length Estimates 


“op 


Los. 
415. 
425, 
435. 
bbs, 
455. 
465, 


Table 9 
PATROL FRIGATE 


3640.5 


Displacenent 
Cutter 

DDO? Model 
3766.3 3707.7 
3652.2 3679.5 
353761 3672.0 
. 3401.6~- 3696.0 
3458.8 3764.1 
35292 3874.7 


3985.9 
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is increased, the displacements predicted by the cuttor 


model begin to diverge from those predicted by DDO7. 

This indicates that the cutter model begins to lose 
accuracy for ship lengths above 400 feet, the assumed 
upper limit of the cutter model. Table 9 also shows that 
the two models predict different optimums. This is a 
result of the assumptions built into the two models. 


Bx — “ e ; 7 ® ‘ : ? 
Me HO eM yr eee 


CHAPTER V 
, RECOMMENDATIONS AND CONCLUSIONS 


For the 378° WHEC and 210° WREC, the full load 
displacements predicted by the cutter model are within 
four percent of the actual full load displacements of 
these two vessels. This is better accuracy than is 
obtainable with other currently available models for ‘the 
size range of Coast Guard vessels. However, only a limited 
number of desigms were used to develop the cutter model 
and these same designs were used to check it. The program 
should be used with caution until other designs are avail- 
able to confirm the relationships used in the model. 

Although the cutter model appears to give accurate 
results in its current form, there are several areas 
where significant improvements could be made. The two 
most important improvements ares (1) development of avail- 
able. arrangement area and volume relationships which are 
more applicable to ships of the size range used in the 
model and (2) expansion of the C, array used in calculating 
horsepower. requirements. The electrical generator sizing 
criteria should also be given careful consideration to 
determine if-the sizes which result are what is really 
desired. :2 0 ius 4. 
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will require considerable developmental work. As a 


minimum, the variation of area with hull depth, bean, 
Coe Cy, tankage volume, and machinery box length must 
be considered. A computer program should be developed 
to calculate tankage volume and arrangements area for 
systematic changes in these variables. The results could 
then be converted into formulas for use in the cutter 
model. 

The expansion of the C, array could probably best 
be accomplished by fairing data from another model test 
series, such as series 64, into the data from Taylor's 
Standard Series. Where necessary sgome data may have to be 
assumed that is not available. An expansion of the 
array to speed length ratios of 3.0 and C. values up 
to 0.008 would be very desirable. With the current 
array limits, many ships can not be run which may be 
very desirable. Also, if the sustained horsepower is 
given as an input, there are many times when the array 
limits are exceeded during an iteration even though the 
final speed value may be well within the limits of the 
array. 

Finally, the generator capacity predicted by the 
model appears to be excessive. This criteria should 
be given careful consideration. Existing 378° WHEC 
cutters seldom exceed the capacity of one 500 KW generator 
in normal steaming; however, the criteria built into the 
cutter model would result in 1500 KW generator sets being 
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installed. 


Some general recommendations on the use of the model 
are also in order. First, a ship synthesis model is only 
a design tool. Complete reliance should never be placed 
on a model; but, if checked by hand calculations, the 
model will greatly inhance the engineer's ship design 
capability. Also, as a tool, the model must be kept 
current. This implies that periodic reviews must be made 
of the empirical relations used to insure that they are 
in agreement with current policy. Every new design built 
should also be plotted on the curves included in Chapter 
III as additional data. Changes should be made if the 
new data is not in agreement with the old. 

It is expected that the model will most often be 
used for designing new ships. Studied most often will 
probably be changes in length and Cpe These values have 
been arranged so that this can be done by changing only 
one input card. However, changes in any of the input 
variables can be investigated. 

If the model is used to investigate the sensitivity 
to changes in the design assumptions, it will be necessary 
to change statements in the program itself. If many of 
these studies are being performed, some of the fixed 
values in the program could be read in on input cards. 

One constraint on all ship synthesis models is that 
variations in hull form can not be investigated. If 


hull forms much different from those of current designs are 
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desired, a new model has to be developed. It must be 
remembered that a ship synthesis model is really a multi- 
dimensional linearization of a complex problem. If 
large departures are made from the data base designs, 
the errors will increase in magnitude. 

Finally, the cutter model could be programmed as 
part of a ship optimization routine, but it is recommended 
that this not be done. Such a procedure puts the designer 
out of a controlling position while putting total reliance 
on an optinization criteria. Optimization criteria are 
difficult to define if more than just the cost of the 
vessel is in question. It is felt that such questions 
as the balance between payload and cost are best left 
to the designer and not to a rigid formula. 





1. 


Ze 


36 


b, 


5. 


6. 


7. 


8, 


205 


REFERENCES 


Mills, J. Feasibility Studies and Ship Synthesis 
Hodels. lLecvures 7, 3, and 9. Unpublighsd 
lecture notes for Professional Summer at H.I.T. 


Hildebrand, F.B. Introduction to Numerical Anal gis. 
New York: McGraw-Hill Book Company, 1956. 


Comstock, John P., ed. Principles of Naval architecture. 
New Yorks The Society ox Naval Architects and 
Marine Engineers, 1967. 


Shen, Michael and Feter Gale A Precedure for Predicting 
the Spsed of Surface Ships in Rough Head Seas. 
Raval ship Engineering Center Report io. 6114- 
71-4, 1970. 


shen, Michael, Peter A. Gale and Kenneth Walker An 
of Survsce Shins in Revza Head Seas. Naval 
Ship Lngineering Center Keport No. 6114-71-8, 
1970. 


U.S. Coast Guard Cost Estimating Procedure, “Planagan's 
Hethod”, U.S. Coast Guard fechnical Files, 
Washington, D.C. 





Gertler, Morton A Reanalyasis of the Original Test 
Data for the Taylor se niard series. David W. 
faylor Hodel Basin Report 306, washington, D.C.s 
U.S. Government Printing Office, 1954. 


Mills, J. Personal notes on hull subdivision for 
destroyer typo hulls 





APPENDIX A 


DESCRIPTION OF INPUT 





APFENDIX A ae 


DeSeCrI plo of Enpuy 


Input for the cutter model can be divided into two 
Main groups. The first group contains the data for the 
Cy arrays used in estimating propulsion power. The second 
group contains the data for individual ships that are run. 

Data for the C, arrays is read only once regardless 
of how many test cases are run. This data precedes all 
other data and the values never change. The OC, array 
data can be stored on tape or disk if desired. A list 
of the values used is given in Appendix D. 

The second group of data allows the user to specify 
the vayload and performance as well asS various other 
parameters of the ship the computer is to design. In 
order for the computer to produce reliable results, the 
input data must be the best that it is possible to specify. 
Past designs are a good starting point for determining 
reasonable input values. 

The minimum number of data cards for the second 
group of input is 14. These cards are described below. 


Where more than one card may be used this fact is noted. 


Card 1 -- A single number is entered on this card 
in card column 1. The values that can be used 
are? 

QO ~ the program will stop 


1 =- this is the normal value to use 











Card 


Card 


Card 
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2-a fine hull form will be assumed for 


the first iteration, C., = mOoO025.< 


2 -- This card contains a heading that is 
printed above the output. All 80 columns 
can be used for any alphanumeric description 


desired. 


3 -- A Single value is entered in card column 2. 
The total number of armament, aircraft and 
cargo payload data cards should be entered. 


A maximum of 20 can be used. 


4 -- This is really a set of cards containing 
the armament, aircraft and cargo payload data. 
The number of cards used is specified on card 3. 
Bach of these cards has the following format: 
(All data should contain a decimal point) 
columns 1--20 Any alphanumeric title 
describing the remainder of the data 
on the card. 
columns 21-25 Cargo weight, enter the 
weight in tons. 
columns 26-30 Cargo vcg, enter the ratio 
of cargo veg divided by the depth of 
the ship at amidships. 
columns 31-35 Group 7 weight,,enter 


weight in tons. 





209 
columns 36-40 Group 7 vcg, enter ratio 


of kg to depth at amidships. 
columns 41-45 Ammunition weight, enter 
weight in tons.. 
columns 46-50 Ammunition veg, enter ratio 
of kg to depth at amidships. 
columns 51-55 Aircraft weight, enter 
weight in tons. 
columns 56-60 Aircraft veg, enter ratio 
of kg to depth at amidships. 
columns 61-65 Electrical load, enter value 
in KW. 
columns 66-70 Deckhouse Area, enter area 
required to be in deckhouse in sq. ft. 
columns 71-75 Hull Area, enter all remaining 
reguired area in sq. ft. 
columns 76-80 Group 7 cost, enter acqui- 
sition cost in thousands of dollars. 
Except for the moment arms, the values given 
on these cards are summed by category and only this 
sum is used in calculations. Therefore it does not 
matter whether all the data is given on one card or 
on several. The only requirement is that the vcg 
values apply to the weights which immediately precede 
them. All values which are included in light ship 
weignt sroup 7 must be given as data on these cards. 


If there is no payload data, card 3 should 


210 
contain a 1 anda blank card should be included 


for card 4, 

Card 5 -- A single value is entered in card column 2. 
This is a count of the number of data cards 
described under card 6. A maximum of 20 can 


be used. 


Card 6 -- This is really a set of electronics input 
data. cards up to a maximum of 20. The number 
of cards is specified by card 5. Hach of these 
cards has the following formats: (All data should 
contain a decimal point) 
columns 1--20 Any alphanumeric title 
describing the remainder of the 
data on the card. 
columns 21-25 enetn 4 weight, enter 
weight in tons. 
columns 26-30 Group 4 veg, enter ratio 
of kg to depth at amidships. 
columns 31-35 Electrical load in KW. 
columns 36-40 Deckhouse area, enter area 
required to be in deckhouse in sq. ft. 
columns 41-45 Hull area, enter remaining 
required area in sq. ft. 
columns 46-50 Group 4 cost, enter acqui- 
Sition cost in thousands of dollars. 
Again the data may be broken up in any 


way desired as long as the correct sets of 
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weight and vcg are givens. On all of these : 


eards:and :.on the cards .described under card 4, 

ward the-aréa given:under-deckhouse .area should 
included only tthe.area.that is required to be 
in the deckhouse. If an area can be either 
in. the deckhouse or in the hull, it should be 
listed. under hull area. 

All items included in light ship weight 
group 4 must be entered on these cards. fThis 
includes standard navigation equipment. If 
there is no electronics data, card 5 should 
contain a 1 and a blank card should be included 

ari Wor-card 6. 

Card 7.--- The number of officers, CPOs and crew 
members ‘are specified on this card. All numbers 
must be’right-justified integers, 
cal > Columns” 1=--5 -number of officers 
rar! columns 6--10 number. of chief petty officers 
-4°"."eolumns 11-15 number of enlisted men 

Card 8 -- This’card‘contains cost indices input. 

All -values should contain a decimal point. 
“4te fs - columns’ 1---5 labor rate in dollars/hour 
vi t-reolumns 6--10 group 1 index | 
so =. eolumns.11-15- group 2 index 
« i!ecolumns 16-20 group 3 index 
“; ©) Columns 21-25 group 5 index 
* “4. Columns 26-30 group 6 index 








The values of the labor rate and cost aidices a 
may be found under Section 3.14, Figures and 

Card 9 -- The number of endurance days for dry 
provisions and aircraft fuel payload are given 
on this card. All values should contain a 
decimal point. 

columns t---5 endurance days for dry 
provisions 

columns 6--10 weight of aircraft fuel, tons 

columns 11-15 veg of aircraft fuel, 
enter the ratio of kg to depth at 
amidships. 

Card 10 -- A single value is entered in card column 
2. Two options are available. If ai is entered, 
use card iia. This option allows the user 
to input sustained speed and have the computer 
Calculate horsepower. If a 2 is entered use 
card 11 0. Here particulars of the machinery 
plant are input and the sustained speed is 
Calculated. Any other value will have an 
unknown effect. 

Card ila -- Use this card type if a value of 1 is 
entered on card 10. The maximum sustained 
Speed, ,endurance speed and range and propulsive 
coefficients at endurance speed and maximum 
speed are entered on this card. All data should 


contain a decimal point. 


KN Pe TRE Re come meee 
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columns i---5 sustained speed in knots 


columns 6--10 endurance speed in knots 

columns 11-20 endurance range, nautical miles 

columns 21-25 propulsive coef. at endurance 
speed 

columns 26-30 propulsive coef. at maximum 


speed. 


Card 11b -- Use this card type if a value of 2 is 
entered on card 10. The particulars of a specific 
machinery plant are entered on this card. All 
data should contain a decimal point. 

columns 1---5 maximum sustained horsepower, 
enter as HP/1000 
colums 6--10 minimum machinery box 
depth in feet 
columns 11-15 machinery box length in feet 
columns 16-20 total weight of group 2, tons 
columns 21-25 specific fuel consumption 
at max. horsepower in 1lbs/HP-hr 
columns 26-30 specific fuel consumption 
at half power in 1bs/HP-hr 
columns 31-35 endurance speed in knots 
columns 36-45 endurance range, nautical miles 
columns 46-50 propulsion auxiliaries 
electrical load in KW 
columns 51-55 propulsive coef. at endurance 


speed 





columns 56-60 propulsive coef. at eae 
speéd 

columns 61-65 acquisition cost of group 2, 
thousands of dollars 


columns 66-70 lub oil weight in tons. 


Card 12 -- This card contains various indut options 
and dimensions. All values except the first 
should contain a decimal point. 

columns 1---5 MTYFPE, a right-justified 
integer is entered here. 
If MTYPE = 0 the program returns 
to read a new card 1 
If MTYP2 = 1 a diesel plant is assumed 
if the machinery particulars have 
not been given 
If hif¥Pe = 2 a CODOG plant is assumed 
If card 10 contains a value of 2, any 
value of MTYPE except O may be used. 
columns 6--10 free surface correction, 
enter value in feet. This doesnot 
have to be a positive number. A 
negative number can be interpreted 
as a KG below the average built into 
the program. 
.cOlumns 11-15 length between perpendiculars, 


enter the ships length in feet. 





Card 


ZL5 
columns 16-20 Cos enter a value for 


prismatic coefficient 

columns 21-25 Cys enter a value for mid- 
ship section coefficient 

columns 26-30 Gii/Beam, a value of 0.1 
is commonly used 

columns 31-35 design and builders margin, 
enter a value as a decimal, i.e. a 
10% margin is entered as 0.1 

columns 36-45 ACr, this is the frictional 
resistance correction used in horse- 
power calculations. A value of 0.0004 


is commonly used. 


Card 12 may be repeated as many times as 
desired with new data each time. Each new card 


is interpreted as a new data case. 


13 -- If itis Hecired to input an entirely 
new set of payload and performance data or to 
Stop, a value of 0 should be given in column 
5 of card 12. This will cause a new card 1 to 


be read, 


14 -- A new card 1 must be included. If new 
data will be given all of the data cards must 
be repeated. A value of O on card 1 will cause 


the program to stop. 
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A sample listing of the group 2 data cards is given 


below. 
pard 1 
1 
Bard 2 
210 Wisc TEST RUN 
Card 3 
D 
Card 4 
Group 700 
Grouo 701 
Group 702 
Group 750 
Anmuni tion 
columns 
Card 5 
6 
Card 6 
Group 400 
Grous 401 
Group 402 
Grouv 403 iz 


2 
0 
Group 404 9 
Group 450 3 


Card 7 
9 8 5? 
5 610 

Card 8 

pret 1. 


5 10 LS 20 


The data used refers to a 210° 


15 =- columns 


DO eho 15516 39 
25 


30 = colums 


WhIEC. 


385. 76. 


50 75 80 


190. 





Card 9 
60. OF 
e e 
peo 
ward 10 
2 
Card 11ib 
4.5 #19. 
e e 
5 10 
Card i2 
feat St 
e e 
5 10 


Card 13 


O. 
e 


15 = columns 
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APPENDIX E 


HULL GROUP WEIGHT CLASSIFICATION OF 1965 
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109 
110 


111 
ea 2 
113 
114 
115 
halle 
117 
118 
119 
120 
121 
ae 2 
ee 3 
124 
di2 5 
27 
128 
150 
ie 1 


274 
NAVSHIPS HULL GROUP WEIGHT CLASSIFICATION OF 1965 


Group 1 = Hull) Structure 


Shell Plating and Planking 

Longitudinal and Transverse Framing 

Inner Bottom 

Platforms and Plats Below Lowermost Continuous Deck 

Fourth and Lower Continuous Decks 

fhara Beck 

second Deck 

Main Deck and Hangar Deck 

Forecastle and Poop Decks 

Gallery Deck 

Flight Deck, Landing Platforms, Special Purpose 
Decks above Weather Deck 

superstructure 

Foundations for Propulsion Plant Machinery 

Foundations for Auxiliaries and Othsr Equipment 

Structural Bulkheads 

Trunks and Enclosures 

structural Sponsons 

Armor 

Aircraft Fuel Saddle Tank Structure 

Structural Castings, Forgings and Equivalent Weldments 

sea Chests 

Ballast and Buoyancy Units, Fixed or Fluid 

Doors and Closures, Special Purpose 

Doors, Hatches, Manholes and Scuvtles - Nonballistic 

Doors, Hatches, Hanholes and Scuttles - Ballistic 

Kingposts and Support Frames 

sonar Dome 

Masts, Towers, Tetrapods, and Service Platforms 

Welding, Riveting, and Fastenings 

Free Flooding Liquids 





200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
a0 
211 
Ze 
“i 3 
214 
215 
216 
217 
218 
219 
250 
251 


300 
301 
302 
303 


351 


275 
Group 2 = Propulsion 


Boilers and Energy Converters (Non-Nuclear) 
Propulsion Units 

Main Condensers and Air Ejectors 
shafting, Bearings, and Propsllers 
Combustion Air Supply System 

Uptakes (Smoke Pipes 

Propulgion Contr<. Equipment (Non-Nuclear) 
Main Steam System 

Feedwater and Condensate System 
Circulating and Cooling Water Systems 
Fuel O11 Service System 

Lubricating Oil Systen 

Nuclear Steam Generators 

Reactors 

Reactor Coolant System 

Reactor Coolant Service Systens 

Reactor Plant Auxiliary Systens 

Nuclear Power Control and Instrumentation 
Radiation Shielding (Primary) 

Radiation Shielding (Secondary) 
Propulsion Repair Parts 

Propulsion Operating Fluids 


Group 3 - Electric Plant 


Electric Power Generation 

Power Distribution Switchboards 

Power Distribution System (Cable) 

Lighting System (Distribution and Fixtures) 
Electric Plant Repair Parts 

Electric Power Generator Fluids 





L400 
LOL 
402 
403 
404 
O05 
406 
407 
4O8 
OQ 
410 
411 
442 
413 
450 
W51 
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Group 4 - Communication and Control 


Navigation Equipment (Non-Electronic) 
Interior Communication Systems and Equipment 
Gun Fire Control Systen 

Countermeasure Systems (Non-Electronic ) 
Electronic Countermeasure Systems (ECM) 
Missile Fira Cont~ol Systens 

ASW, Torpedo Fire Control Systems (Surface Ships) 
Torpedo Fire Control System (Submarines ) 
Radar Systems 

Radio Communication Systems 

Electronic Navigation Systems 

space Vehicle Electronic Tracking Systems 
sonar Systens 

Electronic Tactical Data System 
Communication and Control Repair Parts 
Communication and Control Operating Fluids 





500 
501 
502 
503 
504 


505 
506 


50? 
508 


599 
510 
511 


mae 


at? 
Group 5 - Auxiliary Systens 


Heating Systen 

Ventilation System 

Air Conditicning Systen 

Refrigerating Spaces, Plant, and Equipment 

Gasoline, HEAP, All Liauid Cargo Piping, Oxygen, 
Nitrogen, and Avieticn Lubricating Oil Systens 

Plumbing Installations 

Firemain, Flushinz, Sprinkler, Washdown, and Salt 
Water Service Systens 


- Fire Extinguishing Systvens 


Drainage, Eallast, frimming, Heeling, and Stabilizer 
Tank Systens 

Fresh Water Systen 

scuppsrs and Dsck Drsins 

Fuel and Diesel O41 Filling, Venting, Stowage and 
Transfer Systems 

Tank Heating Systens 

Compressed Air System 

Auxiliary Steam, Exhaust Steem, and Steam Drains 

Buoyancy Control Systen 

Hiscellaneous Piping Systens 

Distilling Plant 

steering Systems 

Rudders 

Mooring, Towing, Anchor and Aircraft Handling Systems 
and Deck Machinery 

Elevators, lioving Stairways, Stores Strikedown, and 
Stores Handling Equirent 

Operating Gear for Retracting and Elevating Units 

Aircraft Elevators 

Aircraft Arresting Goar, Barriers, and Barricades 

Catapults and Jet Blast Deflectors 

Diving Planes and Stabilizing Fins 

Replenishment at Sea and Cargo Handling 

Auxiliary Systems Repair Parts 

Auxiliary Systems Operating Fluids, Gases, and 
Stabilizer Fluids 
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Group 6 - Outfit end Furnishings 


Hull Fittings 

Boats, Boat Stowage, and Handling 

Rigging and Canvas 

Ladders and Gratings 

Nonstructural Bulkheads end Nonstructural Doors 
Painting 

Deck Covering 

Hull Insulation 

Storerooms, Stowages and Lockers 

Equipment for Utility Spaces 

Equipment for Workshops, Laboratories, and Test Areas 
Equipment for Galley, Fantry, Scullery and Commissary 
Furnishings for Living Spaces 

Furnishings: Offices, Control Centers, Machinery Spaces 
Furnishings for tledical, Dsntal, and Pharmaceutical 
Radiation Shielding for Nuclear Support Facilities 
Outfit and Furnishings Repair Parts 

Outfit and Furnishings Operating Fluids 


Group 7 - Armament 


Guns and Gun Mounts 

Ammunition Handling Systems 

Ammunition Stowage 

Special Weapons, Handling and Stowage 

Rocket and ijissile Launching Devices (Surface to Air, 
Surface to Surface, end Sub-Surface ) 

Rocket and Missile Launching Devices (Antisubmarine 
Warfare ) 

Rocket, Missile, and Components Handling Systems 

Rocket, Missile, and Components Stowage 

Torpedo Tubes 

Torpedo Handling and Stowage 

Mine Handling Systems and Stowage 

small Arms and Pyrotechnic Stowage 

Air-Launched Weapons Handling Systems 

Air-Launched Weapons Stowage 

Cargo Munition Stowage 

General Arrengement - Armament Drawings 

Armament Repair Parts 

Armament Operating Fluids 
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